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CHAPTER 1 

THE MOTIVATION FOR ANALOG COMPUTERS 


Analog Computation, based on the modern electronic 
analog computer, is of fairly recent date. The first 
commercially available general purpose electronic 
analog computers appeared on the market in the 1940's. 
These early machines were an outgrowth of an emerging 
electronics technology and a critical need for automatic 
computing machines that could solve complex dynamic 
problems. Slide Rule and manual equation solving could 
no longer be relied upon as a practical approach to 
seeking engineering solutions. Therefore, analog com¬ 
puters became important tools in the design of air¬ 
craft, jet engines, atomic reactors, oil refineries, 
chemical plants, etc. 

Many types of analog computers have evolved over the 
years. The family has included the mechanical differ¬ 
ential analyzer, electromechanical differential ana¬ 
lyzer, and most recently, the iterative differential 
analyzer. 

The analog computer has always had several advantages 
compared to a digital computer. These are primarily 
speed, more simulation capability per dollar, an ability 
to integrate, and an excellent man-machine interface. 


The main feature of the analog computer is that it can 
integrate time-varying voltages. There is no easy 
way to differentiate. Consequently a mathematical 
model of a physical system which is expressed in terms 
of differential equations cannot be solved with the 
machine directly. It is necessary to reformulate the 
mathematical model in terms of integral equations, 
either implicitly or explicitly. 

The analog computer can integrate only with respect 
to time. Thus, a mathematical model which contains 
partial integrals (corresponding to partial derivatives) 
with respect to several variables must be approximated 
by a set of ordinary integral equations with respect to 
time. Of course, computer time need not correspond 
to time in the physical world, although it usually does. 

Before 1959 the analog computer was a synchronous 
machine; all its integrators operated in unison. In 
1959 the DYSTAC 1 was introduced. The name was an 
acronym for 1 dynamic storage analog computer 1 . This 
machine was the forerunner of the iterative differential 
1 

DYSTAC is a registered trademark of CSI. 



Two large SD 80 computers are used in this engineering laboratory to simulate the behavior of 
drone helicopters with different load configurations and various automatic flight control sys¬ 
tems under a wide range of operating conditions . Simulation saved time and money, reduced the 
exposure of personnel and materiel to possible damage when actual flight tests were made, and 
it also eliminated the necessity of waiting for special environmental (weather) conditions . 

(Photo courtesy, Gyrodyne Company of America, Inc.) 
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An SD 80 computer mounted inside a Boeing 
experimental jet transport, selected for 
the NASA sponsored Supersonic Transport 
(SST) Program. In this actual in-flight 
application, the SD 80 is inserted between 
the pilot f s controls and the aerodynamic 
control surfaces of the jet plane. This 
permits the total control system to assume 
the dynamics of any of a wide variety of 
SST types. 

(Photo courtesy, the Boeing Company) 


Analog computers have long served in the 
fields of chemistry and process control as 
a convenient, low cost means to observe, 
analyze, control, and predict the effect 
of varying parameters in a dynamic problem. 
Some of the problems may pertain to enzyme 
reaction, chemical kinetics, continuous 
distillation, heat transfer or transport 
delay -- just to mention a few basic appli¬ 
cations of analog computers. The university 
student shown here is performing a research 
problem involving the effect of potential 
barriers on kinetic energy levels. 

(Photo courtesy, University of California) 


analyzer which appeared in 1960. The iterative differ¬ 
ential analyzer is an asynchronous computer; the 
integrators need not be controlled in unison. They 
can operate independently either in groups or singly. 

In 1962 the analog computer was augmented with digital 
logic. This innovation first appeared in a machine 
called the HYDAC. 1 The HYDAC had a very large 
quantity of synchronous digital logic. Since that time 
the inclusion of a smaller complement of asynchronous 
digital logic has become accepted practice. 

Digital logic can be used for the implementation of 
logical decisions. These are based on results obtained 
from the analog portion of the computer during the 
solution of the problem. Digital logic also can be used 
for mode control of the analog computer. All of the 
Systron-Donner 10/20 and 40/80 series analog com¬ 
puters can be operated in the iterative mode and can 
be augmented with digital logic. 

Analog computers have found widespread acceptance 
in virtually every area of scientific investigation. This 
growing interest in analog computers has created a 
need for complete software, specially designed for the 
beginner and less experienced user. Also, the recent 
addition of digital logic control has greatly improved 
the problem-solving capability of analog computers. 
How this new feature can be used in analog computation 
is thoroughly illustrated in this publication. 

It is the purpose of this handbook to provide students 
as well as experienced computer users with compre¬ 
hensive and up - to - date analog computer software. 
Chapters 2 to 5 develop the basic fundamentals of com¬ 
puter operation and illustrate the solution of elementary 
problems. Chapter 6 provides a detailed description 
of the operating controls and computer logic of the 
SD 10/20 and 40/80 series computers. This informa¬ 
tion serves as useful reference material to problem- 
solutions illustrated in the more advanced discussions 

1 HYDAC is a registered trademark of EAI. 



A student in mechanical engineering is shown 
how to simulate a mass-spring-damper system 
on an SD 3300 analog computer. Stating 
equations is unnecessary. Using block 
programming techniques, illustrated in 
Chapter 4, the student need only understand 
the basic relationship of physical variables 
and constants. Following the program block 
schematic for a given system, the student can 
easily program the problem on the computer, 
observe results on the oscilloscope, arid make 
further parameter adjustments to seek an 
optimum solution. 


which relate iterative programming techniques to the 
Systron-Donner computers. Chapters 7 to 22 develop 
more sophisticated programming techniques and appli¬ 
cations on a progressive basis. 
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On the Apollo Program, Douglas Aircraft 
coupled an SD 40 computer to a large centri¬ 
fuge, The computer calculated and integrated 
errors in human performance during Apollo 
lifting body reentry simulation studies, The 
Systron-Donner computer is seen next to the 
programmer, in upper left portion of picture, 


(Photo courtesy, Douglas Aircraft Co.) 
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Analog computers are now standard computing equipment for classroom teaching and research work 
in colleges and universities, Students in the departments of Electrical Engineering, Mechan¬ 
ical Engineering, Chemistry and Biosciences receive instruction in the use of desk top analog 
computers as basic electronic model builders of dynamic problems. The ease and swiftness of 
presenting a solution on a readout (oscilloscope, XY recorder), and the ability to vary problem 
parameters and observe immediately their corresponding effects, have made the analog computer 
an important teaching aid, 

(Photo courtesy, University of Santa Clara) 
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CHAPTER 2 

BASIC ANALOG COMPUTING ELEMENTS 


This chapter shows how electronic equipment and 
circuits are used to implement mathematical rela¬ 
tionships in an analog computer. 

OHM'S LAW 

Ohm r s law describes the relationship between the 
current through, and the voltage across a passive 
impedance. A passive impedance is a collection of 
passive elements such as resistors and capacitors 
connected together in an arbitrary way. Such an impe¬ 
dance is generally considered (in analog computing) 
to be a two-terminal network which can be denoted by 




For a capacitor 



so that Ohm's law is 



or E = ZI 


where 


Z 


_ 1 _ 

Cs 


where Z is the dynamic impedance of the element. 
Ohm's law states 


E = ZI 


for 



+ E 


C 



SERIES and PARALLEL IMPEDANCE 
Impedances are additive in series: 


where I is the time-varying current (in the direction 
indicated) through the passive element generated by 
the impressed time-varying voltage, E (with the 
polarity indicated). 

In order to simplify what follows, transform notation 
will be used where 


•7 

I 

-7 

Z 1 


Z 2 


z t 1 z 1 + z ; 


+ E, 


+ e 2 


s f(t) = -jjf f(t) 

and 

if = f f(t)dt. 


This can be proved in a simple way. The total voltage 
from point a to point b is E T = Ei + E 2 and by 
Ohm's law 

e t = z x I + z 2 I = (z 1 + z 2 ) I = z T i. 


For a resistor Z = R and Ohm's law is 
E = RI 


When impedances are connected in parallel the total 
impedance can be found as the reciprocal of the sum 
of the reciprocals: 
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—+— 
Zt z, z 2 z 3 


This again can be proved by simple application of 
Ohm’s law: 


= h + h 


I 


T 



1 





The total impedance of two parallel impedances has 
the simple formula: 


Z 1 + Z 2 


Z 1 Z 2 
Z 1 + Z 2 


DIODE 

A diode is a non-linear resistor. Its resistance or 
impedance depends on the direction of the current flow¬ 
ing through it. It is denoted by 

I 

— 

+ E- 


where 


E - ZI 
Z = Rp E > 0 

= R 2 , E < 0 


and generally 


r 2 » R i • 


Typically, Ri ranges between 1 ohm and 100 ohms; 
and R2 between 100, 000 ohms (100K&) and 1, 000, 000 
ohms (IMft). Thus the diode is an approximation to a 
switch for which 


OPERATIONAL AMPLIFIER 

An operational amplifier is of the type that is called a 
d-c amplifier; it amplifies not only time-varying volt¬ 
ages, but also d-c or constant voltages. It is charac¬ 
terized by its excellent stability and extremely high 
low-frequency gain (amplification). It is denoted by 



where 


e = -A e. . 
o 1 

These voltages are measured with respect to ground 
(zero reference). The gain, A, is usually frequency 
dependent and will decrease with increasing frequency. 
As will be seen presently, this becomes a limiting 
factor in the use of a computer at high speeds. An 
equivalent circuit for the amplifier is 


e; o 




Io 


where, generally, Rj is greater than 10° ohms and 
Ro less than 1CT 4 ohms (closed-loop) at zero fre¬ 
quency. However, the amplifier is current limited. 
That is, it will only perform satisfactorily if the out¬ 
put current, Iq, is less than some value. The cur¬ 
rent limit for the S-D* amplifier is ±25 milliamperes 
at ±100 volt output. The amplifier is also voltage 
limited; it will not function satisfactorily if the output 
voltage, e 0 , is greater, in absolute value, than some 
upper limit. The limit for the S-D solid-state 100 volt 
amplifier is 105 volts. Since the output of these amp¬ 
lifiers contains unwanted noise whose magnitude typi¬ 
cally can be 10“2 volts, their effective useful range is 
about three and one-half decades (5x10*2 to 10 2 volts). 


GENERATION OF TRANSFER FUNCTIONS 


In what follows, it will be assumed that the input 
impedance of the operational amplifier is infinite, the 
output impedance is zero, and the gain is infinites 


1 Abbreviation for Systron-Donner. 
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These assumptions introduce negligible error at zero 
frequency and are a good approximation at mid-fre¬ 
quencies. At high frequency, the assumptions cannot 
be made (particularly for gain). 

Consider an amplifier with input and feedback (from 
output to input) impedances and applied input voltages 
as shown. 


and 



The simplest case is 



From Ohm's law 



Now if the input impedance, R., of the amplifier is 
assumed to be infinite 1 


which can be expressed 



By a suitable choice of impedances many desired trans¬ 
fer functions can be generated. (Henceforth, the units 
of megohm and microfarad will be used for resistance 
and capacitance respectively.) 

SUMMER (ADDER) 

Consider the configuration shown below, which is 
called a summer. 



for finite e. . Thus 
J 



Further, if the gain, A, of the amplifier is assumed 
to be infinite, then 



This is a reasonable assumption at low frequencies 
since as noted above e 0 , at most, will be in the neigh¬ 
borhood of 100 volts in absolute value. Thus, equation 
(1) becomes 
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A set of typical values (in megohms) for these resistors 
in a summer in a computer is 


Substitution in equation (3) of 


1 

Cs 


R 1 = R 2 = R f = 1.0 
R 3 = 0.1 

so that 

e o = - e l -e 2 ” 10e 3 • 



In a computer, access is usually provided to the input 
(summing) junction so that additional input resistors 
or feedback components can be added to the summer 
circuit externally. Thus, the summer circuit is 



so that the output voltage, e , of the amplifier will be 
the integral with respect toTime of the time-varying 
input voltage, e^ 

If a switch is included as shown 



where the gains (multiplying factors) are indicated at 
each input, and the summing junction terminal by J. 
The program symbol for the summer is 




where the input gains are omitted if they are unity. If 
it is desired to indicate the junction or summing junc¬ 
tion (high-gain input), the symbol becomes 



INTEGRATOR 

The high-gain input is labeled only if pertinent. Con¬ 
sider the configuration shown below, which is called 
an integrator 


C 



and if 

1 ) there is an initial charge stored on the capacitor 
which results in the voltage, e Q (0) 

2 ) the switch closes at t=0, 

then 

t 

e o = e o (0) ' / e i (t)dt - 

o 

It is common practice to call this initial voltage the 
initial condition' for the integrator. The derivation 
of this terminology is obvious from mathematics. 

It is necessary to find a practical way to establish the 
initial condition or initial voltage for the integrator. 
To accomplish this, the circuit below is used. 



C 


2-4 


Here 


1 _ r , . i_ = HCs + 1 

~ CS R R 

Z - R 
f RCs + 1 

Z. = R 
i 

and substituting in equation (3) 

e o = ^ RCs + 1^ ^' e IC^ = RCs + 1 e lC 


or 


RC 


^0 

dt 


+ 


e 


o 


e 


IC* 


The solution of this differential equation is 
e 0 = A exp ^ + e IC’ 

where A is a constant depending on the initial voltage 
stored on the capacitor before e IC was applied. Then 


Lim p 
t —00 e O 


e 


IC’ 


For practical cases, it is only necessary that t>10RC, 
since e“^<. 0001, which is compatible with the accu¬ 
racy of the circuit. The result in simple terms is: 
to guarantee the establishment of the initial condition 
it is necessary to wait 10RC seconds after the voltage, 
-e^, has been applied to the circuit. 

The two previous circuits, combined with appropriate 
switches as shown below, constitute a practical inte¬ 
grator 


0.1 0.1 



The operation is as follows: 

With S 2 closed and Si open, e 0 = ejc after 
1.0 second (or 0.1 second if C = 0.1, etc.). 


With S 2 open, e 0 still equals ejc, since the 
current through the capacitor is zero. (This 
must be true since both Si andS2 are open and the 
input impedance of the amplifier is assumed to be 
infinite.) 


With S^ closed at, say, t=0, then 

t 

e o (t) = e IC " / e i (t)dt - 

o 

If at t - T, Si is again opened, e 0 will stop chan¬ 
ging and remain at the last value before Sj 
opened, namely 


T 

e o < T > = e IC - / e i (t > dt 
o 

If S 2 is closed again while S^ is open, the out¬ 
put of the integrator will return to 

e = e in . 
o IC 

Each integrator in a computer has these switches as 
part of its circuit. They are open or closed depending 
on what the programmer wants the integrator to do. 
The state of these switches is called the integrator 
mode. If all the integrators are controlled in unison, 
the switch states are determined by the main computer 
mode. 

The modes have simple names with obvious interpre¬ 
tations. During the R, or reset (IC) mode, initial 
voltages are impressed on the integrator capacitors. 
During the compute or C mode, the integrators inte¬ 
grate input voltages. During the hold or H mode, the 
integrator outputs remain constant at the last value 
achieved before entering the hold mode. Thus, the 
mode permits the programmer to stop the computation 
at any time, enabling him to evaluate what has happened 
thus far in the calculation. The table below shows the 
states of Si andS 2 in the various modes. The numeral 1 
indicates the switch is closed (logical 1) while 0 
indicates the switch is open (logical 0). 

MODE S 2 S 2 

R(IC) 0 1 

H 0 0 

C 1 0 

Finally, the general integrator circuit is shown below 
together with a definition of its transfer function for 
various modes. 
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MODE 


TRANSFER FUNCTION 


R 

C 


H 


e o = e IC 


e o (t) = e IC 


- y[e l( t) + e 2 (t) + 10e 3 (t)jdt 


e o( T > = e IC 


- yjej(t) + e 2 (t) + 10e 3 (t)J dt 


Note that the gain (multiplying factor) of the computing 
inputs is noted next to the terminals and is determined 
by the value of the input resistor. The SJ terminal is 
provided to allow the other external input resistors to 
be added. The U terminal provides the ability to gen¬ 
erate an initial condition which is the sum of several 
voltages. The J terminal permits the connection of 
additional feedback elements around the amplifier (to 
be discussed in chapter 10). Also, these terminals 
provide an external connection for the use of solid- 
state switches in place of the mechanical switches 
(relays) Sj and S 2 . In this case Sj and S 2 are con¬ 
strained to be open regardless of integrator mode and 
the necessary switching is done by the external solid- 
state switches. 



Then 

X 2 = r i + l L 

so that 


and 


Since 




e 

o 




+ 


1 




+ 



R 1 + “ R t 


The program symbol for the integrator is 


where R^. 


total potentiometer resistance, 


R 2 R t " R 1 


IC 



The input gains are usually omitted if they are 
unity. 

POTENTIOMETERS 

It is necessary to have a device for entering constant 
parameters in the computer program. This is accomp¬ 
lished with a potentiometer. The circuit for this ele¬ 
ment is shown below. 



and 


Rt 


R, 


R, 


e 

o 


R 1 R L 


R t (r l + R i } " R 1 2 


e = a e.. 
o 1 


Now, obviously, it would be time-consuming to deter¬ 
mine Rj, with the knowledge of Rt, Rl for each new a. 
Consequently, in practice, the pot is set with the load 
connected, by reading the output voltage, e , with a 
meter, foraknown input voltage, which is usually 100 
volts. 

Thus the meter reading is equal to 100 a and the pot is 
changed until the desired value for a is obtained. The 
symbol for the potentiometer is 


If the output voltage, e 0 , of the potentiometer (pot) is 
applied to an input of another element of the computer, 
then the input impedance, Rl, of the other element is 
connected from the pot output to ground. 
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where H is the notation for the high end, or input. The 
His usually omitted since it is obvious from the com¬ 
puter program which side is the input to the pot. The 
element discussed above is called a two-terminal pot 
since it has two available terminals (or connections) 
on the computer program board (patchboard). Some¬ 
times it is desirable to connect the bottom (low) end 
of the pot to some computing-element instead of ground. 
In this case the low end is made available at the patch¬ 
board. The circuit is 



and the program symbol is 




where notation for the arm, A, is omitted when it is 
obvious from the computer program. 

ARBITRARY FUNCTION GENERATOR 

Many mathematical problems to be solved with a com¬ 
puter require the generation of an arbitrary function. 
This is accomplished with a device called an arbitrary 
function generator (or frequently, diode function gen¬ 
erator, since the internal circuitry uses diodes). This 
device allows the programmer to approximate the de¬ 
sired function with straight line segments. An example 
is shown below 


The principle of operation of the function generator 
depends on the use of a non-linear input impedance 
for an operational amplifier. That is, the impedance 
generates an input current proportional to the function 
which in turn constrains the output voltage to have this 
functional relationship to the input voltage. A simpli¬ 
fied circuit is shown below. 


1.0 



The program symbol for an arbitrary function gener¬ 
ator is 



FIXED FUNCTION GENERATOR 

Fixed function generators are used to generate often- 
used functions such as sine, cosine, log. etc. These 
are similar to arbitrary function generators in opera¬ 
tion, but do not permit the programmer to change the 
parameters within the device. Generally, a fixed 
function generator is more accurate and has better 
frequency and noise specifications than an arbitrary 
function generator. The program symbol for a sine 
generator, for example, is 




Sin x (t) 


MULTIPLIER 

The multiplier used with the S-D computer is called a 
quarter- square multiplier. The name derives from the 
equation 


Here the function f(e^)is approximated by f'(ej) with 
three line segments. A In general, each function gen¬ 
erator, depending on how it is used, will provide either 
10 or 11 line-segments with maximum slope changes 
of 2 or 2. 5:1. Several of these devices can be used to¬ 
gether if more segments are required. The location 
of the slope discontinuity is called the breakpoint and 
is adjustable. Detailed instructions for the setup of 
this element appear in Chapter 6. 


XY = j [(X + Y) 2 - (X-Y) 2 ] 

The actual equation to be implemented with hardware 
in order to provide multiplication is 



The multiplier module contains two fixed function gen¬ 
erators. A block diagram for one of the generators is 
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where K is a constant that determines the feedback If module is connected to an output amplifier 
resistor of the output amplifier (to be shown below). 

The block diagram for the other is 



e = ~ e l e 2 

The multiplier module block diagram is o 100 



The program symbol for the multiplier module (with¬ 
out the output amplifier) is 



With the designations 

M: The module is used as a multiplier 
I: The output is a current 
X, Y: Bipolar inputs 
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Note that the "minus sign" associated with the "M" indi¬ 
cates a negative current is generated for a positive 
product. 


If the bipolar inputs are not naturally available from the 
program (e.g. other variables) then summers can be 
used to generate them. This can be done in two ways: 

Case I: 




PROGRAM SYMBOL 


Note that the "plus sign" associated with the "M" indic¬ 
ates a positive current is generated for a positive pro¬ 
duct, e x ey. The sign is generally omitted since this 
is the normal condition. 


The program symbol for the multiplier module together 
with an output amplifier is 



e Q = -■ 


e x e* 


100 


As before, summers can be used to generate the bi¬ 
polar inputs if they are not otherwise available. 

Case I: 



Case II: 



CIRCUIT 



PROGRAM SYMBOL 



PROGRAM SYMBOL 

(Note that here the sign inversion is due to the output 
amplifier.) 

Case II: 
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A multiplier module can be used to make a divider with 
the same output amplifier. (This connection can be 
made conveniently at the patchboard.) 



Here 


I 


e e 

k ° y = 
* 100 



so that for R 


1 _ 

K 


e 


x 


or 


e 

o 


Vy 

100 


100e 


The programmer’s symbol is 


Rf < 0. Then, If will be in the direction shown, for 
e 0 > 0, and If will produce an ej <0. Since the gain 
of the amplifier is negative and large, the circuit will 
be unstable. An additional restriction, |e x | < |e y |, is 
necessary to ensure that the output amplifier not ex¬ 
ceed the voltage limit. 

The divider can be used to generate square root in a 
simple way. If e y - e Q , 



then 


or 



e 2 = -100 e 
o x 



From the circuit it can be seen that the condition ey>0 
must hold in order to avoid instability. This is easy 
to show. The multiplier module can be considered to 
be a non-linear feedback resister, R^. The sign of 
R f is the same as the polarity of e y . Suppose in the 
circuit 


so that 


e o = 10 /sT - v°- 

The diode insures that the system will not saturate in 
the wrong direction if e x goes positive inadvertently. 

This circuit has the program symbol 


e x 



e o 


IO-/e x 
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FUNCTION SWITCH 


COMPARATOR 


A function switch is a manual switch which can be oper¬ 
ated from the control console. It is used to change the 
computer program during execution. The switch has 
three positions: up, down, and center (off). Its circuit 
is shown below. 



D 


This is the program symbol as well. The terminals 
of the switch are located at the patchboard. 

FUNCTION RELAY 

A function relay is used to make program changes auto¬ 
matically during execution. It is energized (logical 1) 
for an input equal to or greater than +28 volts and is 
de-energized (logical 0) for an input equal to or less 
than 0 volts. The circuit is shown below 



This is also the program symbol where N is the relay 
number. 

ELECTRONIC SWITCH 


An electronic switch is used in much the same way as 
afunction relay; the only restriction being that it must 
be in series with the junction of an amplifier. Its speed 
of operation is much greater (10"5 sec) than afunction 
relay (10 _ 3 sec). It requires the same input voltages 
(logic levels) as the function relay (i.e. logical 1: 
28<e. <100, logical 0: -100<e. <0). 

The program symbol is 



where N is the switch number. 


Frequently, it is required to determine the sign of the 
sum of two variables. 

A circuit which accomplishes this is 



If (ei + e2) > 0, then e 0 is limited to approximately 
0 volts by Di because of its low resistance when it con¬ 
ducts (logical 1). If (ei + e2) < 0, then e 0 is limited 
to approximately +28 volts by Dg (logical 0). The 
latter is true since 

1) I F mustbein the direction shown (e 1 + e 2 >0) 
and e^ > ej. 

2) e Q = 28 volts since e^ must be approximately 
0 volts due to the high gain of the amplifier. 

The program symbol for the comparator is 


a 

b 




U 


U = I if (a+b)<0 
U = 0 if (a + b) > 0 


If both logical comparator outputs are required, then 
e Q is complemented with a biased analog inverter. The 
program symbol is 


a 

b 



JJ 

U 


and the analog circuit is 
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CHAPTER 3 

ELEMENTARY ANALOG PROGRAMMING 


All dependent variables in an analog computer are vol¬ 
tages. Consequently, it is necessary to equate these 
voltages to physical variables. Since the elements of 
the computer are voltage limited, scaling will be re¬ 
quired. The art of scaling is discussed in Chapter 5. 
In this chapter scaling is ignored since the necessary 
concepts can be developed without reference to scaling. 

As pointed out in Chapter 1, it is necessary to remem¬ 
ber that an analog computer cannot differentiate easily. 
It is a machine whose main feature is its ability to 
integrate. Consequently, problems which are defined 
by differential equations must be reformulated in terms 
of integral equations either explicitly or implicitly. 
Consider the problem 

y = y(x) 

I * b » - ° 

y(0) = A. 

This can be rewritten as the integral equation 
y + B J ydx = 0 

y(o) = a. 

For this problem computer time will represent x, and 

y(t) + &Jydt = 0, t = x 
y(0) = A 

An integrator, as shown in Chapter 2, integrates from 
t = 0 onward. It also requires an initial condition: 
namely the value of the integral at t = 0, which in 
this case is A. Thus, the problem is formulated for 
the computer as 

t 

y(t) = -f By(t)dt, (1) 

O 

y(o) = a 

The right half of equation (1) can be generated by 


where the constant B and the initial condition A are 
inserted by means of pots. Remember that the inte¬ 
grator inherently has a sign change. Equation (1) states 
that y(t) is equal to the negative integral. All that is 
needed, to complete the computer program for this 
problem, is to connect the output of the integrator to 
the input of the pot set to B. Hence, the computer 
program is 


A/100 



The independent variable, y, can be recorded from the 
output of the integrator. 

In this example the differential equation was reformu¬ 
lated explicitly in terms of an integral equation. It is 
frequently possible to do this in an implicit way, as 
the next example will show. 

Consider the problem 

y + Ay + By = f(t), y = y (t) (2) 

y(0) = C 

y(0) = D 

This is equivalent to 

y = f(t) - Ay - By 

y(0) = C 

y(0) = d 

First, y can be generated from y by 


A/100 
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(Note that again the pots are used to generate initial where equation (2) is satisfied by applying the correct 
conditions.) Using this circuit, the complete program inputs to integrator 1. 
can be generated by 

D/100 C/100 Next, the use of non-linear elements is demonstrated. 
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As another example, consider 


y + f(t)y = o 
y(o) = a 

y(0) = B 

The computer program is 



The next example shows how the DCU f s (Chapter 8) can 
be used to determine the solution of a differential equa¬ 
tion for a specified value of the independent variable. 
Suppose it is required to find y (. 17) for 

y + y = 0 

y(0) = 100. 

From above, the analog program which generates y is 



In order to determine y(. 17) the computer is allowed 
to compute for 0.17 seconds and then put into hold. 
The latter can be accomplished by applying +28 volts 
to the’'problem hold” terminal at the patchboard after 
0.17 seconds. This will halt the entire computer and 
in particular the computation carried out by the above 
program. The digital program which will do this is 
shown below: 



The operation is as follows: The clock pulses are not 
applied to DCU 1 until the machine is in compute be¬ 
cause the output of the RH f bus is a logical 0 in I C 
and this is an input to gate 1. The FF is initially in 
the R state when the machine is in compute having been 
reset by the FRT’ Logic output, divider 1 reduces the 
clock frequency of 1000cycles/sec to 100 cycles/sec. 
The output of the 0 count of counter 4 is 10 cycles/sec 
and is connected to the input of counter 5. The 7 count 
output of counter 4 and the 1 count output of counter 5 
are connected to gate 2. Thus, whenever t =* 0.17 sec 
the output of gate 2 will change from Oto 1, which will 
change the FF to the S state. At this time the computer 
will go to the hold mode due to the S output of the FF 
being connected to the problem hold terminal at the 
patchboard. Also, the clock and counters will be reset 
because of the logical 1 applied to the 0V terminal. 
As soon as the machine goes to the reset mode the 
FRT’ logic output will cause the FF to return to the 
R state. Thus, the logic circuit is ready to beusad 


again as the computer is manually put first in the reset 
mode and then the compute mode. 

ALGEBRAIC EQUATIONS 

The first examples were intended to illustrate the basic 
approach to programming. The following illustrates 
a more serious application of the analog computer. 
Suppose it is desired to solve a set of simultaneous 
algebraic equations which are expressed in matrix form 
by 

AX = C (3) 

where the n x n matrix A and the column matrix C 
are known. For the purpose of illustration it will be 
assumed that A is of rank 2, although the derivation 
which follows is perfectly general. Consider 

X + AX = C 


(3a) 








When the system has reached steady-state (assuming 
there is one) 

X = 0 
AX = C. 

Thus the solution to equation (3) is obtained from the 
steady-state solution of equation (3a). The latter can 
be written 

X 1 + a ll x l + a 12 x 2 = C 1 

x 2 + a 21 x l + a 22 X 2 = C 2 (4) 

and the analog program is 



Generally, integrator capacitors are chosen equal to 
0. 001 fjifd to decrease the solution time. However, 
there is no guarantee that the system of differential 


equations, (4), is stable. A different approach will 
provide an unconditionally stable system of differential 
equations. If 

AX = C 

it follows that 

A'AX = A T C (5) 

where A T is the transpose of A. The solution of equa¬ 
tion (5) is the same as that for equation (3). Again at 
steady-state the solution of 

X + A T AX = A f C (6) 


will be the solution of equation (5) and therefore equa¬ 
tion (3). The stability of equation (6) is determined by 

X + A T AX = 0. (7) 

If equation (7) is multiplied by the row matrix X, then 
the result is 

XX + XA T AX = 0. 

It is well known that XA'AX is positive definite (i. e. 
non-negative for all values of X). Unless X = 0, 

XX< 0 


and stability is guaranteed. The computer program for 
equation (6) is 
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In this program it is easy to change individual coeffi¬ 
cients because each pot is associated with only one co¬ 
efficient. The program can be simplified so that it 
requires no more equipment than the previous pro¬ 
gram if the following coefficient changes are made. 


OLD 

NEW 

a n 

2 2 
a ll a 21 

a 12 

a ll a 12 + a 21 a 22 

a 21 

a ll a 12 + a 21 a 22 

a 22 

a 22 12 

C 1 

a ll C l + a 21 C l 

o 

to 

a 12 C 2 + a 22 C 2 


This simplified program saves equipment but is in¬ 
convenient to use when the solution is required for 
several different values of the coefficients. 


f(t) = A(t + a) n ,* t > 0: 

df _ n f 
dt t + a 

f(0) = Aa 11 

Generating Equation 

Aa n 


100 



Program 


ARBITRARY CLOSED FUNCTIONS 

Many analytic or closed functions can be generated by 
representing them in terms of their generating dif¬ 
ferential equations. The following examples illustrate 
the technique. 


f(t) = In (t + oj, t > 0: 

df = _1_ 

dt t + a 
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f(t) = A cosh(t), t>0: 


This program is 


f = f 
f(0) = A 

Generating Equation 


A 

100 



Mean value : 

The mean value, <p (t), of a function, f (t), over the 
interval T ^ + tj is defined by 

T +t 

*<t) = 4/ 1 f( ^ )d ^ i - Q - (8) 

T i 

The generating equation for <f> can be obtained by dif¬ 
ferentiation: 

3T * t{^^ 1 + ^ ” 0 ft)} (9) 

The program is 



Suppose the,integrators 1, 2 are in the initial condition 
mode until t # = T^, at which time their mode is changed 
to compute. Then, theoretically, the program will 
generate the solution of equation (8). However, the 
program will not work because the right side of equa¬ 
tion (9) is indeterminate for t = 0. Thus the divider, 
in the program will have an unstable output. Practi¬ 
cal limitations require the Y-input to be at least 3 
volts for useful operation. Consequently, it is neces¬ 
sary to choose a different initial condition: 

i Tl+e 

0 («) * — jf f(Ud£, 



As before, when t 1 = T^ the mode of integrators 1, 
2 is changed from initial condition to compute. Initially 
SI is closed and S2 is open. The output of integrator 
1 will not be <f> (t) until t = « . When t = €, the cor¬ 
rect initial condition, namely 0(c), will have been es¬ 
tablished for integrator 1. At this time the comparator 
will be activated, which will open SI and close S2. 
Also the output of integrator 2 will be c. Thus, the 
generating equation is implemented by the program 
starting when t = e 1 . The output of integrator 1 will 
be 0 (t), t>€. 

COORDINATE TRANSFORMATION 

Suppose two coordinate systems have the same origin 
and are displaced by an angular rotation, 6 : 



An X-coordinate point, sayx, will have the coordinates 
x' = x cos e 
y • = -x sin Q 

in the x T , yVcoordinate system. Similarly, for a 
y-coordinate, say y, 

x' = y sin 6 
y' = y cos e. 
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Thus, a vector, R, in the X, Y-system with components 
x, y, will have the components 


positive rotations 0, <p,\ respectively about the X, Y, 
Z axes correspond to the transformations 


x' = x cos 0 + y sin 0 
y/ = -x sin 0 + y cos 0 

in the X', Y' - system. In matrix notation 
R' = TR 

where R, R' are the representations of the vector in 
the unprimed and primed coordinate systems respec¬ 
tively, and where 

T = T cos 6 sin ^1 

[7 sin 0 cos 0 ). 

In the three-dimensional coordinate system 


Z 



1 0 0 


0 cos 

0 

sin 0 

0 -sin 

0 

cos 0j 

COS 0 

0 

-sin <j> 

0 

1 

0 

sin (p 

0 

COS <t> 


R 


z 

A 


cos A sin A 0 
-sin A cos A 0 
0 0 1 


Any arbitrary three-dimensional rotation can be repre¬ 
sented as a product of these three matrices. The super¬ 
scripts refer to the axis about which the coordinate 
system is rotated and the subscripts to the name of 
the angle of rotation. 

x 

The transformation of a vector, V, by R is 


V' = RV 


V' 

x 


= V 


V' = V cos e + V sin e 

y y z 


V' = -V sin 0 + V cos 0 . 
z y z 


The program is 



Note that 0 is limited by the range of the sine and cosine function generators. 
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The program symbol used frequently is 


Note that AB ^ BA. The unit matrix, I, is 



The transformations R y , R Z have similar programs. 
MATRICES 

A matrix, A, is a collection of elements: 

A 

(in what follows it is assumed that n = m). Thus, if 
n = 3 


a ll a 12 a 13 

A = a 21 a 22 a 23 

a 31 a 32 a 33 


I, has the property 


IA = AI = A 


Differentiation: 


d_ A = ^_( a ij) = da 
dx A dx J dx 


or for, say, n = 2 


da.. da.^ 
11 12 

d a ll a 12 dx dx 


^ * a 21 a 22 " da 21 da 22 

dx dx 

Transpose: The transpose, A 1 , of A is 


Multiplication by a constant, a, follows the rule 


A' = a.. 
Ji 


11 

a 12 

a 13l 

-1 

& 

i—*■ 

H- 1 

" a l2 

0^13 

'21 

a 22 

1! 

CO 

CM 

oS 

“ a 21 

aa -22 

“ a 23 

31 

a 32 

a 33j | 

.^31 

aa 32 

" a 33 


which also can be denoted 


aA = o?a.. aA = = raa..] 

iJ L x Jj [ *jJ 


when it is understood that the a .. are the elements of A. 

Addition follows the rule 

A + B = a.. + b.. = (a + b).. 

ij iJ 'i J 


or for, say, n = 2 

A*B-P U *»!♦ 


a ll a 12 b ll b 12 a ll +b ll a 12 +b 12 


21 b 22 b 21 b 22 a 21 +b 21 a 22 +b 22 . 


Multiplication: 


AB = y a. b 

2 , IQ! ( 


or for, say, n = 2 

A' - 


a ll a 21 


a 12 a 22 


Submatrices: The submatrix, Sy, corresponding to 
the element ay in A is obtained by 
deleting the i-th row and j-th column 
of A. Thus for, say, n = 3 


a ll a 12 a 13 

S 13 = a 21 a 22 a 23 

a 31 a 32 a 33 


a 21 a 22 


a 31 a 32 . 


Determinant: The determinant, |A|, of A is a scalar 
quantity and can be found by induction. 
By definition 

N 

| A | = £ (-l) 1+] a. i |s..|, for any i. (10) 
j=l J J 


Thus for, say, n = 2 


a ll b ll + a 12 b 21 a ll b 12 + a 12 b 22 


a 21 b ll + a 22 b 21 a 21 b 12 + a 22 b 22 . 


For n = 1 


(Note that some particular choice is 
made for i). 


I A I - a n- 
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Thus, for n = 2, 

M - - a n s n - “ 12 s i2 

= a ll a 22 " a 12 a 21 


or 

2 

l A l = ^ a 2jl S 2jl = " a 21 S 21 +a 22 S 22 

= " a 21 a 12 +a 22 a ll 

and either calculation leads to the same result. For 
n = 3, 


Cofactor: The cofactor, C.., 

ment a.., in A is 
ij 


corresponding to the ele- 


or for, say, n 


and 


C.. = (-l) 1 +j |s..|, 

1 ] ' ' 1 J 1 ’ 


21 

a 22 

31 

a 32 

a ll 

a 12 

a 31 

a 32 



£<-*>Piil 


and in particular, any of the equations 

| A | = j | 1 (-D 1 +j a lj |S lj | 

IA | .^ 1 (- 1 ) 2+3a 2j| S 2j[ 

l A l = | 1 ( - 1)3 %l S 3 j l 

will lead to the same (A I. If i = 1 


l A| = a ll S ir a 12 S 12 +a 13 S 13 


Inverse: The inverse, A’ 1 , of matrix, A, is 



For example, if n = 2 


A 


-1 


1 _ 

a ll a 22" a 12 a 21 


22 



" a 21 

a ll 


Orthogonality: A matrix, A, is said to be orthogonal 
if 


A' 


A 


-1 


a ll^ a 22 a 33' a 23 a 32^ " a 12^ a 21 a 33 _a 23 a 31^ 

+ a 13 (a 21 a 32' a 22 a 31 ) 

a ll a 22 a 33" a ll a 23 a 32' a 12 a 21 a 33 + a 12 a 23 a 31 

+ a 13 a 21 a 32" a 13 a 22 a 31 
Alternatively, | A| can be determined from 

l A l = I (-l) 1+i a i1 | S |, for any j. (11) 

i=l J J 

A determinant can be multiplied by a constant, a, in 
a simple way: 

|qa| = q?|a|. 

From equations (10) and (11) it can be seen that if a 
matrix B can be generated by m adjacent row inter¬ 
changes and p adjacent column interchanges starting 
with matrix A, then 

|b| = (-i) m+ P |a|. 


Note that R , R , R are orthogonal 
matrices. The geometric interpreta¬ 
tion of orthogonality is rotation. 

Logical Matrix Operations: 

(AB)' = B'A' 

(AB)~ 1 = B~ 1 A -1 

DECOMPOSITION OF ORTHOGONAL MATRICES 

Frequently the programmer is presented with a com¬ 
posite orthogonal (rotation) matrix which is a function 
of several angular variables. Usually these angular 
variables represent successive rotations about coordi¬ 
nate axes. The simplest way to generate the computer 
program is to decompose the matrix into a product of 
the RX, rY ? rZ matrices and then generate subpro¬ 
grams for each of these matrices. This can be accom¬ 
plished by following a simple set of rules. These are 
stated without proof and examples are given. R x , R Y , 
R z are said to be type 1, 2, 3 matrices respectively. 
Suppose the composite matrix, M, can be represented 

M = 
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The outer matrices are A-, A-., and the inner matrix 
Mj, is defined by w 

M i = A 2 A 3 ‘ ‘ ‘ ■'Sr-i . 


The simplest element of Mis Sin 0 and is in row 2 and 
column 3. The products of column 3 contain the vari¬ 
able * and those of row 2 the variable 0 . Therefore 
A l> a 3 are of types 2 , 3 with variables *, 0 , respec¬ 
tively. 


Now suppose the simplest element (lowest degree homo¬ 
geneity) of M is located in row i and column j. Then A. 
is of type i and A^ is of type j. The angular variables 
corresponding to Aj[, A^ will appear in the products 
of columns j, i respectively. The sense of the rotation, 
Ai, (positive or negative) can be found by setting all 
the variables corresponding to A 2 ..., A^ equal to 
zero. The sense of the rotation, Ajf, is similarly 
found by setting the variables corresponding to A, ..., 
Ajsf-1 equal to zero. If the simplest element of M is 
zero, then M has the representation 

M = A 1 a 2 

and there is no inner matrix. If the simplest element 
of M is not zero, then there is an inner matrix and Mj 
can be found by setting the variables corresponding to 
A l> a n ec l ua l to zero. The further decomposition of 
Mj is carried out with the same procedure as for M. 


Example: 


M 


Cos * 

-Cos 0 Sin * 
Sin 0 Sin * 


Sin * 

Cos 0 Cos * 
-Sin 0 Cos * 


0 

Sin 0 
Cos 9 


The simplest element is zero and is in row 1 and col¬ 
umn 3. Therefore there is no inner matrix and Aj, 
a 2 , are of types 1, 3 respectively. The products of 
column 3 contain the variable 9 and the products of 
row 1 contain the variable for *. Therefore, the var¬ 
iable for Ai is 9 and the variable for A 2 is *. If * = 0 


M 


If 9 = 0 


10 0 
0 Cos 9 -Sin 9 
0 -Sin 0 Cos 9 


M 


Cos * Sin * 0 

-Sin * Cos * 0 
0 0 1 


Finally, 

M = R X R^ 
9 * 


Example: 


M 


Cos 0 Cos* + Sin 0 Sin 0 Sin* 
-Cos 0 Sin 0 

Cos 0 Sin* + Sin 0 Sin 0 Cos* 


If 0 = 0 = 0 , then 


A i = 


Cos * 0 -Sin * 

0 1 0 

Sin * 0 Cos * 


If * = 0 = 0, then 


*3 


If 0 = * = 0, then 


A 2 = M l - 


Cos 0 Sin 0 0 

-Sin 0 Cos 0 0 

0 0 1 


10 0 
0 Cos 0 Sin 0 

0 -Sin 0 Cos 0 


Finally, 




As can be seen from the previous example, it is not 
necessary to find all the elements of A j to determine 
the sense of *. Since it is known that Ai is of type 2 
it is only required to evaluate the element in row 1 and 
column 3 with 0 = 0 = 0. Similarly, if the element 
in row 1 and column 2 is found for *= 0 = 0 , then 
the sense of 0 is determined. 


Example: (Here we use the abbreviations: C = Cosa, 
S = Sinof) a 


M 


V.VVAVV. 0 * 

-SWA 


W*VAV C A C > SVW* 
W.'VWVWx ww. 


■s A as-cs 

* x * * x 


C . C _ 
X * 


The simplest element of .M is C X C which is in row 3 
and column 3. The products of column 3 contain 0 and 
those of row 3 contain 0 . Thus, Ai is type 3 and A N 
is also type 3. 0,0 are the variables for Aj, A N 

respectively. Consider the element in row 1 and 
column 2. For X = * = 0 = 0 its value is Sin 0 and 
consequently 


A 


1 



Sin 0 Cos* - Sin 0 Cos 0 Sin* 
Cos 0 Cos 0 

Sin 0 Sin* - Sin 0 Cos <p Cos* 


-Cos 0 Sin* 

, Sin 0 

Cos 0 Cos* 


For 0 = X = * = 0 its value is Sin 0 and consequently 
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To determine the inner matrix, set 0 = <f> = 0 so that 


Finally, 


-e x * 0 




-S,S T 
X * 

c s 

A * 


C x S x C * 
■ s x c x c * 


The simplest element of Mjis a zero in row 1 and col¬ 
umn 2. Also for * = 0 the element in row 2 and col¬ 
umn 3 is SinX. The products of column 2 and row 1 
contain the variables X, *, respectively. Thus 

v y 

m t = Rf r: . 

I X * 


COORDINATE TRANSFORMATIONS WITH UNRE¬ 
STRICTED ANGLES 

The coordinate transformation program is the same as 
before, except for the use of Sine and Cosine function 
generators. For unrestricted range we use the vari¬ 
able 0 and the equations 

Sin0 = J<9Cos0 dt 

Cos0 =-/0Sin0 dt 


The program is, 


-IOO 



The lower two multipliers and summer generate 

2 2 

E = -(100 - Sin 0 - Cos 0). 

This error is used to change the damping of the inte¬ 
grators to ensure orthonormality of Sin0, Cos0. 
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CHAPTER 4 


BLOCK PROGRAMMING FOR PHYSICAL SYSTEMS 

A Valuable Aid to Simulation of Lumped-Parameter Physical Systems 


The analog computer, or differential analyzer, is a 
powerful tool for obtaining dynamic solutions to differ¬ 
ential equations. Generally the equations to be sol¬ 
ved relate to a physical system of particular interest 
to the computer programmer, in which case the pro¬ 
grammer wishes to simulate his physical system by 
means of an electrical network whose defining equa¬ 
tions are analogous (hence the name analog computer) 
to the subject physical system. It is desirable then, 
to arrange the computer mechanization to have a one- 
to-one correspondence between the problem board and 
the physical system, so that the programmer "sees” 
his system rather than an abstract electrical network. 

The following discussion will illustrate a method of 
programming in which the physical elements are 
treated as blocks rather than as sets of differential 
equations. Certain useful basic building blocks will 
be developed, and methods for combining the basic 
blocks to form any desired physical system (simula¬ 
tion) will be discussed. Simulations generated by the 
methods to be described will have the advantage that 
any or all physical parameters appear as single con¬ 
trols in the computer so that any parameter may be 
changed without affecting any other setting. In addi¬ 
tion, any or all variables appear as computer outputs 
and may be individually monitored. 

A lumped-parameter system, by definition, consists 
of a number of interconnected discrete elements whose 
individual or collective responses to impressedforces 
or stimuli are of analytical interest. The forces and 
responses are related by a mathematical operator. 
Block programming starts with the selection of basic 
definition and rules. 

MATHEMATICAL OPERATOR (O) 

The operator, or combination of operators, describes 
the functional relationship between quantities. The 
most common operators are the following: 

Summer 

Constant Multiplier (potentiometer, commonly 
referred to as "pot") 

Inverter - -1 Sign Changer 

Integrator - ^ (La Place Transform Notation) 

Differentiator - s (This operator is never used in 
an analog simulation if it can 
be avoided, as it usually can). 

Variable Multiplier 

Arbitrary Function - f(x) 

Analytical Function - Trigonometric, Log, 

Hyperbolic, etc. 

Inverse Operator - (O) 1 s, —- are inverses 


ELEMENT 

The element is the basic unit of the block representation 
of a system. It has a pair of terminals or nodes, with 
a value associated with each node. It also has a trans¬ 
ference quantity between nodes, or "through" the ele¬ 
ment which is functionally related by some operator to 
the node-pair value. The transference quantity may 
also be referred to as the branch transference or trans¬ 
mission. 

The elements required to depict the passive lumped- 
parameters of most physical systems are very few in 
number. In particular, five basic elements will suffice 
for many linear electrical, mechanical, and thermal 
systems. These are the Summer, Constant Multiplier, 
Sign Changer, Integrator, and Differentiator. 

BLOCK 

A block is the computer mechanization related to a 
physical element. The input/output variables of a com¬ 
puter block are voltages and currents, but the operator 
relationship is the same as the simulated element. 

Table I gives examples of single elements and their 
corresponding blocks. Table II gives a few commonly 
occurring two-element combinations. Combinations 
of more than two elements can be generated as shown 
in the following examples. Useful tables of more com¬ 
plicated blocks or transfer-function simulations appear 
in many publications. * 

DRIVING SOURCE 

A driving source is a source of energy or power for a 
network of elements. 

SYSTEM STABILITY 

A stable system is one whose responses are bounded 
(finite amplitude limit) for any finite input. An un¬ 
stable system is one whose responses are not bounded. 

NETWORK 

A network is an interconnected set of elements and 
sources. Elements may be joined together at their 
nodes, providing the connected nodes have common 
dimensions and value, and provided the transference 
quantities have common dimensions. 

This discussion concerns elements and networks in 
which the following two rules apply: (Kirchoff T s laws) 

1. The algebraic sum of the node values around 
any closed path (loop) in a network is zero. 

2. The algebraic sum of the transmission quan¬ 
tities at any node is zero. 


* Korn and Huskey, "Computer Handbook", Chapter 2, 
McGraw-Hill, New York, 1962. 
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TABLE I. Single-element Computing Blocks 


PHYSICAL 

ELEMENT 


SCHEMATIC 

SYMBOL 


MATHEMATICAL 

RELATIONSHIP 


PROGRAM BLOCK 



Torsional 

Spring 


Heat 

Conductor 


Viscous 

Damper 






















TABLE I (Continued) 


PHYSICAL 

ELEMENT 


Electrical 

Capacitor 


Viscous 

Damper 


Mechanical 

Mass 


Thermal 

Mass 


Electrical 

Inductor 


Mechanical 

Mass 


SCHEMATIC 

SYMBOL 



MATHEMATICAL 

RELATIONSHIP 


( E 1- E 2>=IC 


< X l- X 2> = iB 


(V 1 ' v 2> = IM 


( T l- T 2>=iM- T 


(Ei - E 2 ) 


(v l ‘ v 2> 



COMPUTER PROGRAM MECHANIZATION 
COMMENTS 


1. Operational amplifiers operate with the junc¬ 
tion at virtually zero potential; the currents 
into the junction through all the input and 
feedback paths must sum to zero; and there 
is a sign inversion between the amplifier in¬ 
put and output quantities. 

2. Any closed-loop in a problem mechanization 
simulating a stable physical system wiU gen- 
eraUy have an odd number of sign inversions. 

3. In mechanizing an equation, set the highest 
order derivative only on the left of the equals 
sign to obviate the necessity for differen¬ 
tiation. 

4. The mechanization is accomplished by: 


a) assigning aposition in the mechanization 
diagram to all necessary variables and 
operators. 

b) accomplishing the necessary intercon¬ 
nections so that each output has its pro¬ 
per inputs. 

c) checking to see that all input/output rela¬ 
tionships in the mechanization are satis¬ 
fied. 

EXAMPLES OF NETWORK PROGRAMMING 

The following examples illustrate some of the tech¬ 
niques used in applying the block programming meth¬ 
ods. Examples are given for linear elements only. 
Non-linearities can be included simply by replacing 
the potentiometer corresponding to the non-linear ele¬ 
ment with a multiplier, function generator, or other 
suitable computing element. Backlash, hysteresis, 
stiction, etc can be included as required. 
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MASS - SPRING - DAMPER SYSTEM 



Element - force equations 

Spring Ky F gl = Kj (X 1 - 0) 

Spring K 2 : F g2 = K 2 (X 2 - Xj) 

Damper B: F^ = B (v 2 - v ) 

Forces acting on masses 


Ml 

-F - 
si 

+ F s2 


-F 

- F 

M2 

s 2 

D 


Steps to mechanize 

1. Assume X r -X 2 , -v^ v 2 and F are available 

2. Draw blocks for Mp M 2 

3. Draw blocks for Kp Kg, B 

4. Make input connections to mass blocks to produce Xp -Xg, -Vp v 2 


Initial values at t = 0. 
Xj (0), X 2 (0), Vj (0), 








The above mechanization makes each parameter avail¬ 
able as an individual adjustment. The programming 
is easily accomplished by generating the forces im¬ 
posed on the elements and then operating on the forces 
to obtain velocity and displacement. Any or all of the 
initial conditions, Xj(0), X2(0), vi(0) and V2(0), may 
be zero, or may also be variable. 

Note that the block representing the damper contains 


two summers. The first summer produces (vj - V 2 ), 
and the second summer acts only as a sign changer. 
The two-amplifier, one-pot combination can be reduced 
to a more simple two-pot combination by connecting 
v^ and V 2 through B pots to the mass block inputs. 
The less complex alternate is shown below. This con¬ 
figuration requires fewer amplifiers, but variations 
in Mj, M 2 , or B necessitate two or three adjustments. 



* Do not connect to -v- to get pot setting of B instead 
B ^ 

of . A coefficient pot output cannot be connected 

to another coefficient pot input. This is a practical 
limitation of the computer hardware. 


Rigid Bar Supported on Springs 



Rigid bar of mass M, moment of inertia about C.G. of 
I. Assume 9 small such that 9 «sin 9. Vertical motion 
only. Yq is displacement of C.G. 


. F- + F« + F 

Y = — - --- 

1 0 sM 


F 1 = ’ K 1 Y 1 


F 0 = -K 0 Y. 


K 2 Y 2 


y i = Y o-V 


Y 2 = Y 0 +< 2 e 
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T 


CM 

•O 

J<2 

1 


M 3 C 3 

1 


K 


Mg C 2 


K = 0 


K- THERMAL CONDUCTIVITY 
M-MASS 

C-SPECIFIC HEAT 











E 2 E 3 /R 3 



fIN-T NETWORK 



E 1 E 3 *3 ” I 2 


*2 R x ’ E 3 sC 3 


_ E 2 - E 3 


3 R 2 

Is 

, -«s*y 

1—1 

1 

<v 

«4 - E 2 ■ SC 2 

<; r l 



, , ■[ - h 

i 

a 

(E 1 ' E 4 ) sC- 

-±- 

1 


E 4 = ®i “ l 2 + *3 + V R 3 


E^ = (E^ - E^) + E^; Eg = E^ 










CHAPTER 5 
SCALING 


The analog computer is voltage limited to a practical 
range of operation of ±100 volts as a maximum and ±1 
volt or zero as a minimum. It is usually not possible 
to numerically equate analog voltages to physical vari¬ 
ables. That is, the magnitude of physical variables 
must be scaled to fall within the useful range of the 
machine. This procedure is called amplitude scaling. 

The analog computer is also limited as to the speed with 
which it will solve a problem. Practical bounds on 
solution time are 100 seconds as a maximum and 100 
milliseconds as a minimum. Events in the physical 
world usually occur in time intervals which fall out¬ 
side these limits. Thus, the simulation of real world 
phenomena with the analog computer is ordinarily 
faster or slower. The procedure for relating computer 
time to physical time is called time scaling. 


The general procedure for scaling is: 

1. Generate an unsealed program having a math¬ 
ematical structure which agrees with the pro¬ 
blem (model) to be solved. 

2. Identify the location of all voltages corres¬ 
ponding to physical variables. 

3. Associate a scale factor with each of these 
voltages (e.g. V x = a x x )- 

4. Label the program with the scaled physical 
variables (i.e. a x X) rather than the voltages. 

5. Estimate maximum values of the physical 
variables. 


Correct scaling is an important factor in reducing sim¬ 
ulation errors. An error analysis of any particular 
program is about as complex as the problem which is 
being solved by the program. Consequently, except 
for unusual cases, an analysis of error is not made. 
An estimate of the accuracy of the program can be 
found from check cases for which the answer is already 
known. Error is reduced by scaling the program so 
that all potentiometer values and amplifier gains are 
reasonable. Frequently, in the process of scaling a 
problem, it will be found that some parts of the model 
are not significant and can be eliminated. Scaling is 
a good check on reasonability. 

Scaling is an art, not a science. A good deal of exper¬ 
ience is required to become proficient. Scaling can be 
done in several ways. The techniques discussed in this 
chapter are those in most prevalent use among analog 
programmers. 

If X is a physical variable corresponding to computer 
voltage, V, then a scale factor, a, is chosen so that 


6 . Choose the amplitude scale factors and the 
time scale factor, N, so that all pot settings 
and gains are reasonable. 

In order to accomplish step (4) it is necessary to know 
how input and output scale factors are related for 
various computer elements. 


POTENTIOMETER 


Physical equation: Y = bX 
Scale factors 


Scaled equation 
Program 


V = a X, V = a Y 
x x ’ y y 


v = b-± V 
y ax 

J X 


a X 
x 


o 


a Y 

y 


b a 


V = a • X 

will fall within the practical range of the computer. It 
is good practice to make V as large as possible. If 
tp is physical (real) time corresponding tot c , the com¬ 
puting time, then these two can be related by a scale 
factor, N: 

, ,,m m m 

t = Nt , dt = N dt 
c p’ c p 

It follows that 

/(F)dt c = N/(F)dt p . (1) 

N is chosen so that the problem solution time on the 
analog computer falls within reasonable limits. Prac¬ 
tically, this is accomplished by a choice of N which 
will not result in prohibitive amplifier gains. 


Note that a scale change can be made with a pot. Assum¬ 
ing either or a^ has been established previously, 
the other is chosen so that ba^/ay is a reasonable pot 
setting. The pot can be used solely for a scale change. 
In that case b = 1. 

SUMMER 


Physical 

equation: 

Z = X + Y 

Scale 
factors : 

V = a Z, V = a X, 
z z ’ x x 9 

Scaled 

a a 

equation: 

> 

+ 

>* 

n 

> 
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Now, if only one-gain inputs for the summer are used, 
then it must be true that a z = a x = a anc * the P ro " 
gram is z x y 



Thus, it is seen that input scale factors must be the 
same for a one-gain summer and that a change of 
scale cannot be made through a summer. However, 
a different input scale factor can be used for a 10-gain 
input. The program for this case is. 


-a x X- 

-a y Y_!0 


where it must be true that a x = lOa^. 



INTEGRATOR 

As with the summer, the input scale factors (except 
for the 10-gain input) must be the same. However, a 
change of scale can be made with an integrator. Both 
amplitude and time scaling can be done. Scaling for 
an integrator is derived below, where C is the value 


of the capacitor in 

fjLfd, and R is in megohms. 

Physical 

equation: 

Y -SX dt p 

Scale 

factors : 

V = a Y, V = a X, t = 
y y x x 9 c 

Scaled 

< 

'C 

ii 

equation: 


Thus, the integrator gain, 


1 

RC’ 


must be 


_ 1 _ 

RC “ 


and the program is 


( 2 ) 



MULTIPLIER 


Physical 

equation: 

XY = 

z 

Scale 
factors: 

V = 

a X, V 


X 

x y 

Scaled 

equation: 

V V 
x y 

100 a 

z 


100 * a a z 

x y 

Program: 



It must be true that 


100 a 

z 


From this requirement it is seen that the relationship 
among the scale factors is the same as that for input and 
output multiplier voltages: 


DIVIDER 


a 

z 


a a 

x y 

100 


Physical x 
equation: — = Z 

Scale 

factors : V x = a^X, V y = a y Y, V z = a g Z 


Scaled V a a 

equation: 100^- . = V z 

y x 

Scale 

factor a x 

constraint: a = 100“ 
z a 

y 


Program: 
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SQUARE ROOT 
Physical 

equation: Y = \fx 

Scale 

factors : V y = a y Y, V x = a x X 

Scaled a _ 

equation: V = Y- . 10 \/v 

y “>AT 

Scale 

factor 

constraint: a =10 Va 
y v x 


Program: 


°x 



FUNCTION GENERATOR 


Physical 

equation: 

Y = f(x) 

Scale 

factors : 

V = a Y, V 

y y ’ x 

Scaled 

V 

equation: 

y 

4—1 

II 

>> 
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EXAMPLE: 

Mass-spring-damper system 



Physical 

equations: Spring Force, F g = KX 

dJC 

Damper Force, F Q = B 

P 


Mass ^2^ 

Acceleration, - 5 

dt 2 
P 

M = 1; B = 4 x 10"' 



K = 6.4 x 10 


^ (0) = 0, X (0) = . 95 
P 


Estimates: ~ 8 x 10 ; Xmax. - 1 

\ p/max 

The unsealed program is: 



The scaled program is obtained by assigning literal 
scale factors, determining their numerical values, and 
calculating the coefficient potentiometer settings. 



The \ and— factors appear 
N 

due to the inherent gain of N 
through each integrator, from 
Eq. (1). 


In general, the determination of the numerical values 
involves a trial and error approach. Any convenient 
parameter can be picked for the starting point. If a 
selected value leads to unreasonable potentiometer 
settings or unreasonable amplifier gains, then new 
values may be required. Reasonable pot gains are 
between . 1 and 1. Reasonable amplifier gains are 
between 1 and 20. Amplifier gains should be integer 
values. Non-integer values can be obtained from a 
potentiometer, in combination with gain-of -10 inputs, 
if necessary. 

1. Determination of a. In order to utilize the 
full dynamic range of integrator 2 , set 
aX max = 100 (volts). 

a = 100 
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Potentiometer 2 setting; 


2. Potentiometer 1 setting. 


4. 


PI 


a X (0) 
100 


.95 


(3) 


3. Determination of ~ At this time, a trial 

value for N maybe selected. Analog solution 
times of . 1 to 100 seconds, or computer nat¬ 
ural frequencies from . 1 to 100 radians per 
second are the preferable operating ranges. 
An estimate of the physical time or frequency 
ranges can be used to select a value for N. 

For this example, a value of N = .01 appears 
reasonable. 

_b = 100 (volts) 

N (dX/d^)max. * ' 

Looking ahead, note that integrator 2 will 
aN 

require a gain of —. The gain of N is in¬ 
herent in the integrator from Eq. 1. There¬ 
fore, the value for the integrator will be 
p K is highly desirable that this factor be 

an integer, preferably either 1 or 10. There¬ 
fore, b = 100 will be used. The effect of this 
is to force integrator 1 to utilize something 
less than its full dynamic range, but the effect 
is not serious in this example. As a general 
rule, it is not possible to scale so that all 
amplifiers work over their full range. 


P2 = — = (6.4 x 10 -3 )c (5) 

NaM 


The value of c is chosen to assure that a rea¬ 
sonable value for P2 is obtained and to assure 

that^ is an integer. 

c = 100 will satisfy both conditions. 
P2 = .64 

5. Potentiometer 3 setting: 


P3 


cB 

NbM 


.04 


This setting for potentiometer 3 is below the 
desired minimum potentiometer setting value 
of . 1. In this example, however, it is not pos¬ 
sible to have a large setting for potentiometer 3. 
This follows from a consideration of the phy¬ 
sical problem. The damper, B, was chosen 
to have a very slight effect on the mass-spring 
system, so it necessarily follows that the 
portion of the analog program related to the 
damper will have only a slight effect. It should 
be noted that abnormally low potentiometer 
settings will frequently arise in simulations of 
systems containing negligibly small elements. 
It should be recognized that the low settings 
stem from the actual system characteristics. 
The programmer, therefore, need not spend 
time in futile attempts to improve the scaling. 


The final program with numerical values is: 


-100 



X s 


100 




d 2 X 


3-f = IO-8(V p2 +V p3 ) 


dt 
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GENERAL COMMENTS ON TIME SCALING: 

The following points concerning time scaling are of 
particular importance: 

1. If all integrator gains are simultaneously 
changed by a factor k, the time scale changes 
N 

from N to-^, and all computer frequencies 

increase by k. The amplitude scaling, how¬ 
ever, does not change. All numerical scale 
factors, all potentiometer settings, and all 
computer problem voltages remain the same. 
This property of amplitude invariance under 
time scale change allows the programmer to 
easily speed up or slow down the computer 
solution. 


2. The most convenient way of changing all inte¬ 
grator gains simultaneously is by changing 
the capacitors by a fixed factor. Integrator 
capacitors can be simultaneously changed by 
means of the computer Time Scale control*. 

3. An event which occurs in a fixed real-time 
interval can occur in several different com¬ 
puter time intervals depending on the choice 
of integrator capacitors, that is, depending 
on the time scale factor. 

4. Events of identical character which occur in 
different real-time intervals can all be made 
to occur in the same computer time interval 
by appropriate time scaling. 


On the SD 10/20 and SD 40/80 analog computers, a master Time Scale switch is a standard feature. This per¬ 
mits fast and convenient switching of capacitors that are patched — up to 1000:1 time scale (xlO, xlOO, xlOOO, 
depending on patch panel connections). 
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CHAPTER 6 

COMPUTER OPERATION 


This chapter is based on materials from the Systron-Donner 10/20 and 
40/80 Computer Operating Manuals. The descriptions of the computer 
mode controls, computer logic, and various computing modules permit 
the reader to relate all applications and programming examples given in 
this handbook to an existing computer. For example, Chapter 8 offers 
detailed information on the use of digital logic. The modules needed to 
implement the applications described in Chapter 8 are shown in this 
chapter. 

Not all available S-D computing components are described here, only 
those which illustrate and serve to explain a basic function. For the 
latest, up-to-date availability of all computing elements and modules in 
the SD 10/20 and 40/80 series, please consult Systron-Donner or your 
nearest SD engineering sales representative. 



SD 10/20 Computer 


The Model SD 10/20 is an all solid state analog com¬ 
puter with an operating range of ±100V and expansion 
capabilities for up to 24 operational amplifiers. The 
principle features of the computer are full ±100V com¬ 
puting range, all solid-state design, compact modular 
construction of all computing components, removable 
problem board which mates directly with computing 
modules to eliminate all problem board cabling, a 
control center with readout selectors for amplifiers 
and potentiometers, high-speed reed relays and solid- 
state switches, complete logic capability for repetitive 
and iterative operations, 0.01% accuracy of computing 
resistor and capacitor networks, a panel voltmeter 
having a null mode of measurement with 0.02% (of full 
scale) resolution. The amplifiers are protected against 
accidental grounding of their outputs as are all power 
supplies. 


All computing components in the Model SD 10/20 are 
modular plug-in units. They are identical to those 
used in the larger SD 80. This high degree of standard¬ 
ization between SD 10/20 and 40/80 computers maxi¬ 
mizes usefulness of computing modules. Since all 
modules are standard plug-in units, the computer user 
is offered a wide choice of computing components. 
As new modules become avai lable, the SD 10/20—40/80 
series computers can thus be constantly up-graded to 
include the new modules. Furthermore, a user may 
start with an SD 10/20 and later use the same 10/20 
modules in an SD 80. This common design feature 
results in important cost savings. As the need for 
more computer modules increases, the expansion cost 
from an SD 10/20 to an SD 80 is minimized. 



m Panel Meter Ranges: 1, 3, 10, 30, 100, and 300 volts, and 
± null. Full scale accuracy: 3%. Null position provides 
0.02% F.S. resolution with reference potentiometer having a 
±0.05% linearity at 25°C. 

iyi Compute Time Selector — compute time continuously vari- 
able from 5 msec to 10 sec. Reset time varies from 5 msec 
to 5 sec, depending upon coarse steps of compute range. 

m Meter Range Selector — with positions for 300 v, 100 v, 30 v, 
10 v, 3 v, and 1 v. Serves also as sensitivity adjustment for 
± null. 

m Function Selector — for rapid choice of: -f null, — null, 
Meter, External (connects selected bus to external jack). 

[5] Address Selector — address of all amplifiers and potentiO' 
meters. 

nn Null Reference Potentiometer — provides high accuracy 
readout using null method with 0.02% F.S. resolution. Lin¬ 
earity is ±0.05% at 25°C. 

[7] Function Switches — provide manual switching flexibility 
in problem solutions. 


nn Mode Selection (lighted pushbuttons): 

Hold — places problem solution on all integrators into 
hold position. 

Compute — applies problem voltages to all integrators. 

Reset — applies initial condition voltages to integrators. 

Rep-Op — places integrators into a repetitive operation 
cycle. Compute time variable from 5 msec to 10 sec. 

Bal/Pot Set — disconnects junction and grounds the 
input resistor summing junctions of all amplifiers. 
Each amplifier is converted to a gain of 2500 for 
precision monitoring of junction offset. 

Pwer On/Off — energizes and de-energizes computer. 

Oven — indicates + 28-volt oven power is on to maintain 
constant temperature of computing capacitors. 

OL Hold, OL Reset — lights up when any amplifier is 
overloaded. When depressed, computer goes into 
Hold; when released, normal operation is resumed. 

Time Scale — activates relays in each integrator module 
to change computing capacitor, (x 10, x 100, x 1000, 
depending on patchpanel connections.) 

Slave — permits operation of computer control circuitry 
from a second console. 

[9] Coefficient Potentiometers — up to 24, available in groups 
of 6, featuring 10-turn wire-wound with lockable counting 
dials. 
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SD 40/80 Computer 


The same compact design and high performance com¬ 
puting capability found in the small SD 10/20 computer 
is also available in the large SD 80 computer. With an 
expansion capability to 84 amplifiers (126 special 
order), 125 potentiometers, and a full set of digital 
logic control modules, the SD 80 provides a problem¬ 
solving power capable of handling large, complex 
problems normally associated with the simulation of 
dynamic electrical, mechanical, thermal and similar 
systems. 

Yet, the large-capacity SD 80 fully retains the 
versatility of a desk-top computer. Mounted on a desk 
or table equipped with lockable casters, it can be 
easily moved to new locations. In addition to creating 
a highly compact and rugged computer, the all solid 
state design increases reliability, simplifies mainten¬ 
ance, and eliminates the need to operate the computer 
in an air conditioned laboratory environment. 

Programming a problem on a fully expanded SD 80 
is virtually as simple as operating a small SD 10/20. 
Visual computer circuits in symbols everyone under¬ 
stands show at a glance how a circuit is patched. 
This feature saves time, reduces patching errors, and 
makes patching a task that the non-expert can perform. 

To allow maximum operator control over a computer 
with a capacity as large as that of the SD 80, all com¬ 
puter controls are centralized in the left wing which 


is movable to any convenient position to suit the 
operator. Major SD 80 controls that speed and simplify 
problem checking and solving include: 

4 A 199-point digital address pushbutton selector 
system for rapid access to amplifier outputs and 
potentiometer arms. 

4 Static and Dynamic checks by simple pushbutton 
control to verify patching and test time scaling 
rapidly. 

4 Instant alarm controls, both visual and audible, 
with automatic transfer to HOLD at the occurrence 
of an overload. 

4 Individual control of compute and reset times for 

precise adjustment in Iterative and REP-OP modes. 

4 Time scaling the problem with a single switch on 
the Control Wing. 

4 One pushbutton to check any amplifier's balance 
without repatching and to set any potentiometer. 

The problem board is removable, just as it is with 
the SD 10/20. There is no restriction as to grouping 
of the modular patch panels in the overall problem 
board layout. In both the SD 10/20 and 40/80 compu¬ 
ters, digital logic control is provided by modular plug¬ 
in units that slip into the central patchbay. They thus 
form an integral part of the computer and the main 
problem board, not a separate added-on element. 



1. Digital Voltmeter — 5 digit readout, automatic ranging and 
polarity. 

2. Address Selector Panel — with an amplifier and potentio¬ 
meter address capability of 199 points. 

3. Control Panel — for selection of: 

Op. Modes: Hold, Compute, Reset and Rep-Op. 

Master Time Scale Switch. 

Compute Time: 5 msec to 5 sec. 

Reset Time: 5 msec to 5 sec. 

Slave: for remote control. 

Oven and Power pushbuttons. 

Check Mode switches for Static and Dynamic checks. 
Overload Controls for: 1) audio alarm, 2) transfer to hold 
mode, 3) removal from Hold mode. 

Amplifier Balance and Pot Set: changes all amplifiers to 
gain of 2500 for simplified monitoring of junction off¬ 
set; and for proper loading of potentiometers. 


4. Function Switches 

5. Coefficient Potentiometers, No. 1-5 

6. Null Meter Reference Potentiometer 

7. Removable Patchboard. — An SD 80 board can hold up to 

42 individual computing module panels. All modules 
are mounted in a universal patchbay directly behind 
the removable problem board. Nomenclature on patch 
panels matches textbook circuit diagrams. Operator 
can thus "see” the circuit while patching a problem. 

8. Built-in Digital Logic Modules, consisting of Flip-flop, 

Gate, and Time/Event Control Modules. 

9. Potentiometer Wing can hold up to 6 panels for a total of 

125 potentiometers (first 5 pots are on Control Wing.) 
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Computer Logic 


Each operating mode switch energizes a logic buss 
which is then distributed to the logic relays in the com¬ 
puting modules (primarily to the integrators). The logic 
levels in the Analog sections of the computer are: 

+28V = True = logical 1 
OV = False = logical 0 

A logical one energizes a buss and places the com¬ 
puter in that mode of operation. All contradictory modes 
have their switches interlocked to prevent two modes of 
operation from being energized at the same time. 

The table that follows outlines the logic busses and 
operating modes. 

In most cases the nomenclature for each logic buss 
is self-explanatory. The following explanations should 
clarify the remainder of the logic nomenclature. 


Logic Buss 

Description 

F. R. (R. 0.) 

Forward Reset (Repetitive Operation). 
This is the normal Reset buss and also 
is the Reset buss energized in the reset 
portion of Repetitive operation. 

R. R. (R. 0.) 

Reverse Reset (Repetitive uperation). 
This buss is the logical complement of 
forward reset. It is used for complemen¬ 
tary integrators and track/store circuits 
in iterative computation. Energized in 
the Compute mode of operation and in 
the Compute portion of Rep-Op. 

F. R. (R. T.) 

Forward Reset (Real Time). This buss 
is used to reset non-repetitive operation 
integrators. It is energized only when 
the reset switch is engaged. 

F. H. 

Forward Hold is the normal Hold buss. 

It is energized in the Hold mode of 
operation or if an overload occurs when 
the Overload Hold switch is engaged. 

R. H. 

Reverse Hold. Reverse Hold is the logi- 
cal complement of F.R.R.T. It is used for 
complementary integrators and Track/ 
Store circuits in iterative computation. 
Energized in the compute mode of opera¬ 
tion and in Rep-Op. 


Table 1. CONTROL LOGIC TABLE 


LOGIC 

BUS 

COMPUTER OPERATING MODE 




POT-SET 

BALANCE 

RESET 

COMPUTE 

REP-OP 

MANUAL 

HOLD 

PROBLEM 

HOLD 

O.L. 

HOLD 

STATIC 

CHECK 

DYNAMIC 

CHECK 

TIME SCALE 
CHANGE 

Bal 

+28 

0 

0 

0 

0 







+28 

0 

0 

0 

0 


+28 




F R (R 0 ) 

0 

+28 

0 

0, +28 

0 

0 

0 




R.R/R.O.) 

0 

0 

+28 

+28, 0 

0 

0 

0 




F.R.(R.T.) 

0 

+28 

0 ^ 

0 

0 






F.H. Bus 

0 

0 

0 

0 

+28 

+28 





R.H. Bus 

0 

0 

+28 

+28 

0 






Static Ck. 








+28 

0 


Dyn. Ck. 








+28 

+28 


Time Scale 










+28 


( 1 ) 

( 2 ) 

(3) 

(4) 


Blank indicates either 1 or +28 condition can exist, two values indicate cyclic operation. 

Manual Hold is actuated by pushbutton switch on Logic Con ^ ro ‘ Pan * L , ... p , . o , 

Problem Hold is actuated by application of +28V or greater Problem Hold Trunkline on Problem Board... 

sri. ox. —<..*« - >—> ~ -> - 

the computer. 


I 


Use of the Mode Controls 

The mode controls are of the push button type with 
interlocks to prevent two contradictory modes from being 
energized at the same time. 

All control in the computer is accomplished in the 
integrators. Thus, when a mode is energized it opens 
or closes a relay in each integrator. Reference to the 
illustration, Integrator Block Diagram, on following page 
(6-6), while reading the following explanations will 
assist in understanding the control system. 

Reset 

The Reset or Initial Condition (1C) mode enables 
the integrators to be set to their initial values before 
starting computation. When the Reset relay is energized, 
an input and a feedback resistor are associated with the 
integrator amplifier making it, in effect, an inverter. 
Applying a voltage to the input resistor (in this case the 
1C terminal) will cause that input to appear at the output 
and to simultaneously charge the integrator capacitor 
to the same voltage. 

Compute 

In the compute mode there are no energized relays. 
The logic rs such that if the computer rs not in Hold and 
not in Reset, it rs in compute. 

Hold 

In the Hold mode a normally closed relay contact is 
operfed which disconnects the input, signal from the 
amplifier. As there rs no discharge path for the inte¬ 
grator capacitors the output will remain at its last com¬ 
puted value. In this manner outputs may be read or re¬ 
corded and the problem resumed at will hy returning to 
Compute. 

Rep-Op 

In the Rep-Op (Repetitive Operation) mode the com¬ 
puter is automatically switched between the Reset and 
Compute modes of operation. This mode enables the user 
to obtain repetitive solutions which are easily monitored 
on an oscilloscope. 

The amount of time the computer stays in the compute 
mode is determined by the setting of the Compute Time 
switch located at the lower left of the control center. 
The outer dial of the switch may be positioned in any of 
three ranges: .005 to .1 sec, .05 to 1 sec, and .5 to 10 
sec or in the External (EXT) position. The inner dial 
adjusts the compute time over the range selected. 

The compute and reset times may be adjusted arbi¬ 
trarily by observing the solution to the problem and ad¬ 
justing the times until the desired display is obtained. 
If more precise compute and reset times are required the 
Compute-Reset cycle may be observed on an oscillo¬ 
scope. Connect the oscilloscope to any of the FOR¬ 
WARD RESET terminals on the integrators or to Test 
Point 6 at the rear of the computer. A rectangular wave¬ 
form will be observed. The positive (+28V) portion is 
Reset, the negative (0V) portion Compute. The Reset 


and Compute time controls may then be adjusted to give 
the desired times. 

If the Compute Time switch is placed in the External 
(EXT) position external circuits may be used to generate 
the Rep-Op cycle. The drive signals must be compatible 
with the logic used in the computer (+28V = True, 0V = 
False) and should be connected to the.R (reset) and 
C (Compute) terminals located on the problem .board. 

OL Hold - OL Reset 

When engaged, the Overload Hold portion of this switch 
will cause the computer to automatically switch to the 
Hold mode at the occurrence of an overload in any ampli¬ 
fier. When the overload is found and eliminated, the 
switch is re-engaged to reset the overload hold circuitry. 

If an overload occurs, the light within the OL Hold 
/OL Reset will begin to blink on and off and will con¬ 
tinue to do so until the overload is removed. 

Time Scale 

The Time Scale switch determines which of the 
integrator capacitors connected to the A and B terminals 
on the integrator modules is in the circuit. With the Time 
Scale switch not engaged, the capacitor connected to 
the A terminal is in the circuit. When the Time Scale 
switch is engaged the capacitor connected to the B ter¬ 
minal is in the circuit. Any two of the three available 
capacitors may be patched to the A and B terminals. If 
time scaling is not required the capacitors may be patch¬ 
ed directly to the output. 

Special Trunklines 

There are six special trunklines located on the pro¬ 
blem board of the Computer. These are located on the 
top left row of modules. The first module contains the 
M, X and Y special trunk lines and the third module 
contains the H, C and R special trunklines. The special 
trunklines are used as follows: 

M: This line is connected to the panel meter or the EXT 
jack as determined by the position of the Meter Select 
Switch. If the Meter Select Switch is in the METER 
function the M line is connected to the EXT jack; if 
the Meter Select Switch is in the EXT function,the 
M line is connected to the panel meter. 

X & Y: The X and Y lines are connected directly to 
the X and Y terminals on the left side panel of the 
computer. These lines are intended for use with an 
X-Y plotter but may be used to connect signals to 
any external device or they may be used to connect 
external inputs into the computer. 

H, C & R: The H (Hold), C (Compute), and R (Reset) 

lines are used as external input lines to the Control 
Logic circuitry. A logical 1 applied to any of the 
lines will place'the computer in that mode of opera¬ 
tion. To use these lines the computer must be in 
REP-OP and the Compute Time switch must be in 
EXT. 
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Integrator Modules 



50K U 50K 



Integrator Block diagram illustrating flexibility 
of circuit and patching convenience. 


Model 3320 

Dual integrator with two uncom¬ 
mitted operational amplifiers. 


Model 3329 

Quad integrator, dual multiplier, 
operational relay. Contains four 
integrating networks, two four 
quadrant multipliers and one 
operational relay. Al so avai lable, 
less multipliers, Model 3329A. 


# • 



Moc/e/ 3320, Dual Integrator 


Model 3329 - Quad Integrator, 
Dual Multiplier, Function Relay 


The Model 3320 is divided into three sections: the 
top and bottom integrators and a logic section common 
to both. Basic patching of the integrator requires only 
a patch from the end of the integrating capacitor to the 
A terminal of the time scale relay. By patching the end 
of another capacitor to the B terminal a second time 
scale is set up and is controllable by the TIME SCALE 
switch on the control panel. 

Initial conditions (1C) are applied at the 1C terminal. 
The source of the 1C voltage may be a pot, amplifier 
output, the reference voltage, or from an external source. 

The logic section consists of a 7-terminal matrix 
near the center of the patching block. The nomenclature 
is the same as that explained in the Logic Control sec¬ 
tion of this manual. The two center terminal s of the matrix 
are connected to the coils of the Reset and Hold relays. 
For normal and repetitive operation, patch from the F 
(forward) terminals to the relay coils. For reverse and 
iterative operation, patch from the R (reverse) terminals 
to the relay coils. For nonrepetitive operation, patch 
from the F(RT) terminal to the reset relay coil. (In this 
mode the integrators will reset ONLY when the reset 
button on the control panel is pushed.) Note that the 
logic patching controls both integrators in this module. 

To use the integrator as a summer, patch from the 
input side of one of the resistors to the output of the 
integrator. This will connect a resistive feedback around 
the amplifier; the other resistors may be used as summer 
inputs. 


The Model 3329 has four integrating networks without 
operational amplifiers. The operational amplifiers neces¬ 
sary for integration are normally patched into the Model 
3329 from an adjacent Model 3325. The input resistor 
network of the Model 3325 then becomes the input to 
the integrator. Relay logic connections are in the lower 
left hand section of the Model 3329 patchboard. Nomen¬ 
clature is the same as that used on the Model 3320, Dual 
Integrator. Each Integrator Network has a Hold (H) and 
Reset (R) terminal associated with it, these terminals 
must be individually patched to the Relay Logic termi¬ 
nals. The Squid (multiple) patch cords should be used 
when more than one integrating network is to be used 
at the same time. 

The Dual Multiplier section of the Model 3329 may 
be patched for division, squaring and square root func- 
tions as well as for multiplication . For multiplication 
or division operation, three operational amplifiers are 
required. In the division and square root functions arbi¬ 
trary limitations as to polarity and relative amplitude of 
the inputs are imposed. For example, in division the 
Y input must always be positive and equal to or greater 
than the X input. Note that in square root an external 
diode must be patched into the circuit. 

Patching and operation of the Function Relay in the 
Model 3329 is identical to that of the Function Relays 
in the Model 3322A. The relay is energized by a nominal 
+28 volt signal and will not be damaged by a signal of 
up to ±100 volts. 
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Summer Modules 



Model 3321 — Dual Summer 


Model 3325 — Quad Summer 


The Model 3321 contains two operational amplifiers and 
the components necessary for summing operations. Basic 
patching consists of connecting one end of the feedback 
resistor to the output. The top summer has three input 
resistors, the bottom summer two input resistors. There 
are four spare summing resistors located in the center 
section of the module which may be used as either input 
or feedback resistors with either of the summers. 

The module also has two 20 kohm resistors and 
two uncommitted diodes associated with each summer. 
These may be used to patch the summer as a limiter- 
comparator or for any of the common diode circuits used 
in computation. 

A function switch termination is provided on each 
summer module. The function switch is a single pole- 
double throw type switch and is physically located near 
the control center. Normally only the first three summing 
modules will have switches terminated in them. 

The center section of the summing module has ter¬ 
mination for a diode function generator. Use of these 
terminations and of the diode function generator is ex¬ 
plained in that section of this manual. 


The Model 3325, Quad Summer is identical in per¬ 
formance to the Model 3321, Dual Summer. As illustrated, 
patching of the Model 3325 is, however, considerably 
different than that of the Model 3321. Basic patching for 
any of the four summers in the module consists of patch¬ 
ing the SJ’ terminal to the J terminal and patching the 
right end of the feedback resistor to any of the output 
terminals. Any of the resistor input networks may be 
used with any amplifier. Thus summers with six, nine 
or twelve inputs may be patched. 

The Model 3325 also contains four spare resistors, 
four diodes and four trunk lines. In addition, it has 
terminations for a Function Switch and for a Diode 
Function Generator. 

The Model 3325 is also used as a source of amplifiers 
for the Model 3329. When used in this manner, the Model 
3325 furnishes the amplifier and input resistors for the 
integrator circuit. 
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Model 3322A - Dual Inverter, Dual Function Relay 



Model 3322A: Dual inverter, dual operational relay. Two 
uncommitted operational amplifiers with 
two operational relays. 


The Dual Inverter contains two operational amplifiers 
for use as inverters or high gain amplifiers in comparator 
and limiter circuits. Operation as an inverter requires 
two patches: (1) From the summing junction to the junc¬ 
tion. (2) From the input side of one of the resistors to 
the amplifier output. The inverter may be used as a 
summer by patching additional resistors into the sum¬ 
ming junction. (Note: Potentiometers used as inputs 
to the Model 3322A should be set in the Reset mode of 
operation.) 

The Model 3322A also contains two function relays. 
Each is a two form C (2FC) relay corresponding to a 
double pole-double throw (DPDT) switch. The relays may 
be energized by any voltage from +28 volts to+ 100 volts. 
The energizing source is typically the output of an amp¬ 
lifier or a logic signal, but any source within the +28 
volt to +100 volt range will energize the relays. Any 
voltage in the range plus and minus 100 may be applied 
without damage. 

As in the Dual Summer, the Dual Inverter contains 
uncommitted diodes, limiter resistors, terminations for 
a Function Switch and Diode Function Generator ter¬ 
minations. The use and operation of these elements is 
the same as that described in the section on the Model 

3321. 
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Mo del 3323 - Dual Inverter, Dual Multiplier 



Model 3323: Dual multiplier, two four quadrant quarter 
square type multipliers with two uncom¬ 
mitted operational amplifiers. 


Multiplication Operation. The multiplier/divider is 
based upon the quarter-square principle, in which the 
following relationship holds: 

1/4 [ (X + Y) 2 - (X-Y) 2 ] = XY 

In ±100 volt computers, the multiplier output voltage is 
always scaled to -0.01 XY. 

The circuitry consists mainly of two diode squaring 
networks, one based upon positive polarities and the 
other upon negative. For four-quadrant operation, volt¬ 
ages representing X, -X, Y, and -Y must be present at 
the inputs to the squaring networks regardless of the 
polarity of input voltages. 

The center section of the module contains the termi¬ 
nations for the multipliers. The section is divided into 
two identical halves, one for each multiplier. The multi¬ 
plier will also perform division, squaring, and square 
root operations 


0.1 



Division Operation. Division is accomplished by the 
terminal connections shown in the figure below. The 
concept is illustrated by the circuit diagram in the same 
illustration. If the output of amplifier 1 is designated Z, 
the output of the multiplier must be+.OlYZ. By reference 
to the basic discussion on operational amplifiers (para¬ 
graph 2.8.)/ it follows that X = -.01YZ, or Z = -100X/Y. 
Since +.01YZ and Z must always be opposite in polarity 
to X, it also follows that Y must always be positive in 
order to maintain this relationship. 

If Y occurs in the problem only as a negative voltage, 
apply it to the -Y multiplier terminal and patch the input 
and output terminals of the inverting amplifier to the -Y 
and Y terminals, respectively. Z still equals -100X/Y. 
The absolute value of Y must always be larger than X; 
otherwise, the output would tend to be larger than 100 
volts, causing an overload. 



CIRCUIT FOR DIVISION OPERATION WITH INSERT SHOWING 
HOW TO CONVERT FROM MULTIPLIER PATCHING. 


Squaring. Squaring is performed as a special case 
of multiplication where X = Y. Only one input inverter 
amplifier is required. Patch together the R and Y input 
terminals and the -R and -Y terminals. Leave the other 
connections as for multiplication. The output voltage 
taken at the amplifier output terminal represents-0.01X2 

Squareroot Operation. The squareroot operation is a 
special case of division where Z = Y. The relationship 
Z = -100X/Y of division becomes Z 2 = -100X for -100 
<X<0, Z= 10VTT, for 0 <X<+ 100 , z = -ioV5r 

If X is always negative, patch together the R and Y 
terminals and the -R and -Y. If X is always positive, 
patch R to -Y and -R to Y. Leave all other connections 
as for division. 

The Inverter sections of the Model 3323 are similar to 
those described for the Model 3322A. The two basic 
patches required for inverter operation are: From the 
summing junction to the junction and from the input 
side of one of the resistors to the output of the amplifier. 


CONNECTIONS POP MULTIPLICATION 
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Model 3324 — Dual Inverter, Quad Electronic Switch 



Model 3324: Electronic switch, four high speed elec¬ 
tronic switches with two uncommitted 
operational amplifiers. 



The center section of the Model 3324 contains four 
electronic switches. There are two normally open (NO) 
and two normally closed (NC) switches. The normal volt¬ 
age required to excite the switch is a nominal +28 volts; 
however, voltages of up to ±100 volts will not damage 
the switch. The current required to excite the switch is 
less than three milliamperes at 28 volts. 

The Electronic switches were designed and are nor¬ 
mally used to replace the reed type relays in the inte¬ 
grator modules. The Electronic switches may be used as 
normal switches as long as the input voltage and current 
through the switch limitations are observed 

Patching to the Model 3320, Dual Integrator is shown 
in the accompanying illustration. Patching to the Model 
3329, Dual Integrator is similar. Also included in the 
Model 3324 are two .001 mfd. capacitors and four .01%, 
10 K ohm resistors for use in high speed Sample-Hold 
circuits. 

The inverter sections of the Model 3324 are identical in 
operation and performance to those of the Model 3322A. 
Patching is identical to that of the Model 3322A and the 
instructions for that Model may be used. 



Integrator with Electronic Switches 


Electronic Switch Block Diagram 
(Connections shown for comparator amplifier) 
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Digital Logic Modules 




Model 3327 

Twelve three-input AND 
Gates with four relay drivers. 


Model 3326 

Twelve RST Flip-flops with 
four relay drivers. 



Model 3326A 

Six RST Flip-flops, six AND 
Gates and four relay drivers. 



"““TloSel 3328 . 

Time-event Control. 1 kHz 
clock with three decade 
dividers, three decade counters 
and one gated relay driver. 


Digital control, the new way of multiplying the efficiency 
of an analog computer, can be included in the SD 40/80 as 
well as the small 10/20 computers. The SD hybrid computer 
expansion system is comprised of three types of compact plug¬ 
in modules. 

The advantages gained by digital control in an analog com¬ 
puter are of far-reaching significance. Here are some important 
new advantages made possible by SD’s Digital Logic Control: 

1. Track and hold operation by individual integrators. 

2. Sub-routines can be flexibly programmed at different speeds de¬ 
pending on decisions made by logical equations. 

3. Program statements can be arranged into a flow chart quite sim¬ 
ilar to those used in digital computation. 

The flexibility gained through this interplay of analog/ 
digital equipment results in: 

1. Better and greater problem-solving capacity. 

2. Ability to solve a wide range of problems that before could not 
easily be handled by an analog computer. 

3. Speed. Problem solving time is greatly reduced. 

Through digital logic control, sub-routines start and termi¬ 
nate when the corresponding binary control variables change 
state as logical functions of: 

1. External control (switches, relays controlled by external devices). 

2. The states of timers or sub-routine counters. 

3. Analog-comparator decisions. 

The interplay of binary control variables and analog com¬ 
putation results in a special hybrid analog-digital structure. 
Relays or electronic switches implement analog sub-routine 
changes under control of digital (binary) control variables 
and constitute the digital-to-analog interface of the computer. 
Analog solutions, in turn, can modify digital control. 


Additionally, in the SD system it is possible to insert a 
delay of the Hold-mode command to the Complementary 
Integrators. This delay is important because it enables the 
C Integrators to store the final values of their inputs in a 
problem solution where the C Integrators are tracking rapidly 
changing problem variables. 

The Reset-Compute-Hold modes are, of course, settable to 
any desired intervals, and it is this new degree of freedom 
imparted to integrators that illustrates the importance of 
digital logic control. 

The combination of the three SD digital logic modules 
(Flip-Flops, Gates, Time/Event Control) results in a most 
flexible hybrid analog-digital structure. The operator can 
easily set all Reset-Compute-Hold intervals of integrators with 
the Time/Event Control Module. 

A typical example that illustrates the use of SD’s digital 
logic is mode control of iterative integrators. In iterative 
operation (IO), results obtained during or at the end of one 
solution of the problem are used to change parameters or the 
circuit configuration (switching) for the next solution. The 
following figure shows how in IO, integrators are paired into 
normal and complementary (opposite) logic to implement 
iterative solutions: 


Normal 


Complementary 



Mode Duty Cycle for Integrators in Iterative Operation 


* When programming involves the use of digital control, logic levels for the digital logic elements are the 
opposite of those for the mode control and function relays. 

* Refer to Chapter 8, “Basic Operation of Digital Logic Elements,” for complete operating descriptions 
of the Systron-Donner digital logic control modules. 















































Function Generator Operation 


Principle of Operation: The Function 

Generator approximates the curve of an arbitrary function 
with a series of connected straight-line segments as 
shown in Figure 1A. This graph represents the sum of 
the individual line segments shown in Figure IB: Each 
line segment corresponds to the output voltage of a diode 
network in the function generator. The point along the X 
axis where each line begins corresponds to the selected 
bias or “breakpoint” voltage. The diode networks are 
arranged so that six conduct when the X input is positive 
(0 to +100 volts) and six when the X input is negative 
(0 to -100 volts). For better accuracy, the input voltage 
may be biased to use all 12 networks above or below 
zero. Two function generator units can be operated as 
a single-channel of 24 segments for more accurate 
simu lation. 

The slope of each line segment corresponds to the 
ratio of output voltage per volt-change of input and can 
be varied between the values of+2 and -2 by adjusting 
the current flow through the diode network. For steeper 
changes of line slope at any point, two or more networks 
can be operated in parallel by shifting to the same 
breakpoint. 

When the individual outputs of all diode networks are 
added together by an operational summing amplifier, the 
total output over the entire range of X-input will repro¬ 
duce the original-straight-line approximation. No output 
is produced when X = 0; therefore, a constant voltage, 
Yo/ must be added at the output summing amplifier, when 
Y is not equal to 0 at X = 0. 

Set Up Procedure: A function may be set up on the 

function generator either by the preparation of a table of 
values to which the individual slope and breakpoint 
potentiometers are adjusted, or by adjusting the potentio¬ 


meters so that a trace of the output voltage on an XY 
plotter will match a prepared graph as closely as pos¬ 
sible (visual method). The second method may be faster 
but is limited in accuracy by the read-out instrument. 




TABLE 2. DATA FOR SETTING UP FUNCTION Y = COSX 


X (DEGREES) 

Y (COSINE) 

X (VOLTS) 

Y (VOLTS) 

CONTROL 

BREAK 

★ 

SLOPE 

POINT 

-90 

0 

-90 

0 


-6 

-6 

-60 

+0.500 

-60 

+50.0 

-6 

-5 

-5 

-45 

+0.707 

-45 

+70.7 

-5 

-4 

-4 

-30 

+0.866 

-30 

+86.6 

-4 

-3 

-3 

-20 

+0.940 

-20 

+94.0 

-3 __ 

» -2 

-2 

-10 

+0.985 

-10 

+98.5 

-- 


-1 

0 

+ 1.00 

0 

+ 100 

±i 



+ 10 

+0.985 

+ 10 

+98.5 

+ 2 ^ 7 ^ 

* + l 

+ 1 

+ 20 

+0.940 

+20 

+94.0 

+3 

► +2 

+2 

+30 

+0.866 

+ 30 

+86.6 

+4 

+ 3 

+ 3 

+45 

+0.707 

+45 

+70.7 

+5 

+4 

+4 

+60 

+0.500 

+60 

+50.0 

+6 

+5 

+5 

+90 

0 

+90 

0 


+6 

+6 


* Arrows indicate sequence of adjustments. 
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SCALING: X: 1 VOLT = 1* (INPUT) 

Y: 100 VOLTS = 1.00 (OUTPUT COSINE VALUE) 

Figure 2 Typical Graph prepared for setting up the Function 
Y = COSX. 

Tabular Method (for Model 3352) 

General: The maximum slope change for all segments, 
except the first, Kx, segment, is ±2.5 : 1. The first, 
Kx, segment has a maximum slope change of ±5:1. 
Breakpoints may be stacked to achieve greater slope 
changes. 

Step 1: Set all Break-point/slope change (polarity) 

switches to the correct combination of polarities as de¬ 
termined by the function to be simulated. (See previous 
discussion on tabulating a function.) 

Step 2: Set all Slope potentiometers to approximately 
center position. Set all Break-point potentiometers to 
one end as determined by the arrows associated with 
the polarity switches. (E.g., if the polarity switch is in 
the + BP, +AS or in the -BP, -AS position the break¬ 
point potentiometers should be set to the counter-clock¬ 
wise end.) During this procedure the output amplifier 
may go into overload. Although no harm is being done, 
the computer may be placed in the Pot Set/Bal mode to 
prevent overload. The computer must be returned to 
Reset before proceeding with the adjustments. 

Step 3: Set the Kx/Yo switch as determined by the 
function and by the following instructions: 

a. The Yo potentiometer sets the value of Y when X 
equals zero. This value may be set anywhere within 
the range of plus to minus 100 volts. 

b. The Kx potentiometer sets the slope of the first line 
segment. The Kx slope is adjustable within the range 
of plus to minus five (5). 

c. If the Kx/Yo switch is placed in the Off position 
both Yo and Kx will be zero. 


d. If the Kx/Yo switch is placed in the Kx position Yo 
will be zero, Kx will be adjustable over its range. 

e. If the Kx/Yo switch is placed in the Yo position, 
Kx will be zero, Yo will be adjustable over its range. 

f. If the Kx/Yo switch is placed in the Kx & Yo 
position both Kx and Yo are adjustable over their 
ranges. 

When the switch has been set, the Yo adjustment 
should be made. With X = 0, the potentiometer may be 
adjusted for the proper value of Y q at the output of the 
Model 3352 (0 2 ). The Null voltmeter or a precision ex¬ 
ternal meter should be used for this and all succeeding 
voltage adjustments in this procedure. 


Step 4: The Kx adjustment should be made next. 

Apply the value of X at the first break-point (BP) to the 
input of the Model 3352. (If a potentiometer is used to 
obtain the X value its output should be applied through 
an inverter to the Model 3352. This will prevent errors 
due to the diode function generator loading the potentio¬ 
meter.) Adjust the Kx potentiometer until the proper value 
of Y is obtained at the output of the Model 3352. 

Step 5: The first break-point (BP) is now adjusted. 
Adjust the first break-point potentiometer until the Kx 
value previously set is offset by approximately 150 mi 11 i- 
volts. (0.15 volts) This offset compensates for the diode 
characteristics and will help produce a more accurate 
function. 

Step 6: Set the input to the value of X at the second 

break point. Adjust the slope potentiometer associated 
with the first break-point until the output is equal to the 
value of Y at the second break-point. 

Step 7: Adjust the B.P. potentiometer associated with 

the next segment until the value of Y set in the preced¬ 
ing specification is offset approximately 150 millivolts. 
(See discussion in Step 5). 

The remaining slope and break-point potentiometers 
are set alternately as indicated by the large arrows on 
the panel, by repeating Steps 6 and 7 until all segments 
have been adjusted. 



CHAPTER 7 
LOGICAL ALGEBRA 


It is necessary to have a rudimentary understanding of 
logical algebra in order to use logic elements. The 
fundamentals of logical algebra are presented for this 
reason. 

LOGICAL VARIABLES 
A logical variable, L, has two values: 

1 (true) 

0 (false) 

ADDITION (OR) 

The logical sum, L^ = L^ + Lg, is defined by the fol¬ 
lowing table: 


L 1 

L 2 

L 3 

0 

0 

0 

0 

1 

1 

1 

0 

1 

1 

1 

1 


This is commonly called the OR function and has the 
alternative symbology 

L 3 l i u L 2 • 

MULTIPLICATION (AND) 

The logical product, L^ = L^ • is defined by the 
following table: 


L 1 

L 2 

L 3 

0 

0 

0 

0 

1 

0 

1 

0 

0 

1 

1 

1 


This is commonly called the AND function and has the 
alternative symbology 

L 3 - L 1 n L 2 . 


COMPLEMENT 

The complement of L is denoted by L. Loosely speak¬ 
ing, L is the inverse of L. That is, whenever L = 1 
or 0, then L = 0 or 1, respectively. In manipulating 
logical formulas, the two obvious relations below are 
often useful. 

L + L - 1 

L • L - 0 

COMMUTATIVE LAWS 

L 1 + L 2 = L 2 + L 1 
L 1 ' L 2 = L 2 L 1 

ASSOCIATIVE LAWS 

L 1 + < L 2 + V = < L 1 + L 2 } + L 3 

Lj • (L 2 • Lg) = (Lj ' L 2 ) • L 3 

DISTRIBUTIVE LAW 

Lj • (Lg + Lg) = L^ • Lg + L^ Lg 

COMPLEMENTED SUM 

T-TT^ - Lj Z 2 


This relationship can be seen from the following table: 


L 1 

L 2 

L l +L 2 

L l+L 2 

L 1 

L 2 

V L 2 

0 

0 

0 

1 

1 

1 

1 

0 

1 

1 

0 

1 

0 

0 

1 

0 

1 

0 

0 

1 

0 

1 

1 

1 

0 

0 

0 

0 


COMPLEMENTED PRODUCT 

L 1 L 2 = L 1 + L 2 
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DE MORGAN FORMULAS 



These can be obtained by complementing the comple¬ 
mented sum and product formulas. 

EXCLUSIVE SUM 

The exclusive sum, ©, of logical variables is defined 
by 

L 1 ® L 2 = ^ L 1 + L 2 ^ ' ( L 1 ' L 2 ^‘ 

This is also called the exclusive OR. is the 

shaded area of the "set" diagram below. 
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CHAPTER 8 

BASIC DIGITAL LOGIC ELEMENTS 


Digital logic can be included in an analog computer. 
It is used for making program decisions based on events 
occurring during the analog simulation. It is also used 
for program mode control. The application of these 
elements is discussed in other chapters of this hand¬ 
book. This chapter describes the basic operation of 
the logic elements. 

OPERATION OF DIGITAL LOGIC IN AN ANALOG 
COMPUTER. 

Logic levels for the elements are 

+28 volts = false (logical zero) 

0 volts = true (logical one) 

Most of the elements have a logically complemented 
output in addition to the normal output. All of the ele¬ 
ments except the drivers and the dividers have status 
lights. When the light is on, the element is in the 
"true" state. 

All of the elements, with the exception of the drivers, 
are designed so that their outputs can be connected 
together without damage. This feature is very impor¬ 
tant because it permits an OR operation by connecting 
two or more outputs to the same point on the patch¬ 
board. (See driver exception below.) 

However, when two outputs are connected together in 
this way, they cannot be used as separate logical var¬ 
iables. To avoid this difficulty diodes can be used. 
Thus, if Li, L 2 are two logical outputs and Li + L 2 , 
Li, L 2 are all required for logical inputs elsewhere 
in the program, the circuit below can be used 



Each of the logic elements will drive any reasonable 
number of others. However, the logic elements are 
not capable of activating analog control relays or func¬ 
tion relays. For this purpose drivers are supplied. 
The driver acts as both a power buffer and logical in¬ 
verter. It should be noted that the logic levels for 


mode control and function relays are the opposite or 
complement of those for digital logic. That is 

0 volts = false (relay de-energized) 

+28 volts = true (relay energized) 

Thus, the driver not only delivers the current neces¬ 
sary to operate several relays but also provides the 
logic level inversion required. To use a driver, the 
output of the logic element is connected to the driver 
and the output of the driver to the relay input. These 
connections are made readily at the patchboard. It 
should be noted that the driver is a special element 
and its output cannot be used as an input to the flip- 
flops or DCU T s. It can be used as a gate input. 

Important: The driver outputs must never be connected 
to any other digital logic output. Con¬ 
necting a dr iver output to ground (logical 1) 
will destroy the output transistor, that is, 
if connected to a gate, it should not be 
OR-ed with any other element. 

The remainder of this section describes the basic oper¬ 
ation of the digital logic elements. 

DRIVER 

The driver is used as a power buffer and logical in¬ 
verter from digital logic elements to analog mode con¬ 
trol, function relays, and electronic switches. Its 
logic levels are 

INPUT OUTPUT 
0 (false) 0 volts 

1 (true) +28 volts 

The program symbol is 


GATE 

The logical gates are of the coincidence type. With the 
exception of input number 1, they operate as AND gates 
only on the inputs connected . That is, if no logical 
variable is connected to one of the inputs at the patch¬ 
board, then this input is essentially de-activated. Note 
that this is the same as connecting a logical 1 to the 
input. 
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The gate is shown on the patchboard as 


Input number 1 is different because it must be connec¬ 
ted; otherwise the gate will not operate. However, if 
only input 1 is connected, the gate will act as an in¬ 
verter if its complemented output is used. 

The program symbol for the gate is 




(NORMAL) 

(INVERTED) 


where the input numbers are shown. 


The truth tables below will clarify the operation of the 
gate. 


1 

INPUTS 

2 3 

NORMAL 

OUTPUT 

INVERTED 

OUTPUT 

none any any 

INPUTS 

0 

NORMAL 

1 

INVERTED 

_ 1 _ 

2 3 

OUTPUT 

OUTPUT 

0 

none none 

0 

1 

1 

none none 

1 

0 



INPUTS 

NORMAL 

INVERTED 

1_ 

2 3 

OUTPUT 

OUTPUT 

0 

0 none 

0 

1 

0 

1 

0 

1 

1 

0 

0 

1 

1 

1 

1 

0 



INPUTS 


NORMAL 

INVERTED 

i_ 

2 

_ 3 _ 

OUTPUT 

OUTPUT 

0 

none 

0 

0 

1 

0 


1 

0 

1 

1 


0 

0 

1 

1 


1 

1 

0 



INPUTS 


NORMAL 

INVERTED 

1 

2 

3 

OUTPUT 

OUTPUT 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

1 

0 

0 

0 

1 

1 

0 

1 

0 

1 

1 

1 

0 

0 

1 

1 

1 

1 

1 

0 



The true state response to inputs 1 and 2 is delayed 
by 50 usee, while input 3 is not. As will be seen in 
Chapter 9 , covering logic applications, this permits 
the gates to be used in connection with flip-flops to 
form a shift register. Also, this feature frequently 
obviates the necessity of using one-shots for delays 
in a logic program. 

Each gate has three inputs on the patchboard. If the 
logical product of more than three variables is desired, 
the additional elements can be connected in through 
externally patched diodes, as shown below. 



The additional elements can be connected to any of the 
three patchboard inputs, except number 1. 

The gate can also be used for pulse-shaping. It is ON 
(logical 1 at the output) when no input voltage is greater 
than 4 volts and OFF when any input is greater than 
6.5 volts. When the net input is between 4 and 6.5 
volts, the output corresponds to the last state of the gate. 

GATE-DRIVER 

This is a special element found in the TIME/EVENT 
CONTROL MODULE. It is a power buffered, logically 
inverted AND gate which performs the combined func¬ 
tion of an AND gate and driver. The program symbol 
is 



where Li, L2 and L3 and e 0 are respectively logical 
input variables and the output voltage. This gate is 
the same in logical operation as the AND gate except 
there is no delay introduced at any input. When L]/ 
L2* L3 = 1, the output of the gate is 28 volts; other¬ 
wise it is zero. 
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The primary use of this gate is for implementing ana¬ 
log mode control changes based on logic outputs. For 
example, in this way a multiple decade counter can be 
used either to stop the analog computer after a pre¬ 
scribed number of runs or in conjunction with a clock 
to specify precisely the length of time the analog com¬ 
puter is in one of several modes. This is accomplished 
by patching the output of a gate-driver to either the C, 
R or H inputs on the patchboard. 

DECADE COUNTING UNIT (DCU) 

The DCU is either a decade divider or a decade counter. 
It contains ten internal binary states which correspond 
to integers from 0 to 9 . Note that this is different than 
a binary counter which counts in the binary number 
system. There are two types ofDCU’sin the comput¬ 
er. Both of these are in the TIME/EVENT CONTROL 
MODULE (also referred to as Counter Module). One 
has all the integer binary states available as a M ten 
line” output. This element is called a decade counter. 
A status light is associated with each of the integer 
binary states. The other type has only the decade out¬ 
put available and is called a decade divider. 

There are two master resets for each group of 3 decade 
counters together with 3 decade dividers. These are 
the OV and +V terminals. The OV input is directly 
coupled to the clock. As long as a logical one is applied 
to this terminal the clock will be stopped. A state 
change from zero to one will cause all the counters 
and dividers to be reset. As long as a logical one is 
maintained at this terminal, all other reset signals are 
inhibited . 

A state change from logical 1 to 0 at the +V terminal 
will cause all the counters and dividers to be reset. 
(This is called a trailing edge trigger.) The +V input 
does not inhibit other reset signals, and it does not 
turn off the clock. 

All counter or divider inputs (including reset inputs) 
require a change from a logical 0 to a logical 1 for 
operation. Thus, the application of a logical 1 to an 
input will cause it to operate once. It will not operate 
again until the logical one is removed and reapplied. 
This is true also for the master resets for each count¬ 
er/divider group. 

The program symbol for the decade counter is 



The counter outputs are the ten integer counts from 0 
to 9 . Whenever a logical 1 is applied to the reset ter¬ 
minal, the counter will return to count 0 . Whenever 
a logical 1 is applied to the input terminal, the counter 
will advance to the next higher count. The counter is 
cyclic: when it is in count 9 and receives an input 
pulse, it returns to count 0. 


Note that the counters can be used in series. This is 
accomplished by connecting any output from the first 
counter to the input of the second counter. Then, when¬ 
ever the transition from alogical 0 to a logical 1 occurs 
at the output of the first counter, this transition from 
0 to 1 is applied to the input of the second counter. 
This causes the second counter to advance one count. 
Normally, the second counter receives its input from 
the ” 0 ” digit of the first. Then, when a transition from 
count 0 to count 1 occurs in the first counter, the out¬ 
put from count 0 goes from a logical 1 to a logical 0 
and the second counter will not be advanced in count: 
its input is activated only with alogical change of state 
from zero to one. Note that when the first counter 
goes from count 9 to count 0 , then the count 0 output 
goes from 0 to 1, and the second counter is advanced. 


The program symbol for the decade divider is 


IN - 
R- 


DIVIDER 


• OUT 


The divider has ten internal integer counts in the same 
way as the counter. However, these are not available 
as logical outputs. An input to the reset terminal will 
cause the divider to return to count 0 internally. The 
output will only return to state 0 if it is in state 1 . If the 
logical variable applied to the input terminal changes 
state from 0 to 1, then the divider will count inter¬ 
nally. The divider is cyclic: count 0 follows count 9 . 

With a sequence of input pulses, the output is alogical 
1 except when the internal counter is in count 8 or count 
9 . Thus, when the internal transition from count 9 to 
count 0 occurs, there is a transition from a logical 0 
to 1 at the output which can be used as input to other 
logic elements. However, it should be remembered 
that if the divider is reset during the time when inter¬ 
nally it is in count 8 or count 9 , the transition from 0 
to 1 will occur at the output then. 


CLOCK 

Each TIME/EVENT CONTROL MODULE provides a 
1000 cycle clock output at the patchboard. 

Maintaining a logical 1 at the OV master reset termi¬ 
nal will stop the clock. 

The program symbol for the clock is 



The logical cycle of the clock is repetitively a logical 
1 for 800 /i sec and logical 0 for 200 n sec. 
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FLIP-FLOP 

A flip-flop (FF) is a one-bit memory. It has a set (S) 
and a reset (R) state. The set state corresponds 
to true (logical 1 ) and the reset state to false. The 
program symbol is 







There is a status light connected to the S output so 
that when the FF is in the S state the light will be on. 
The FFcan be placed in the S state by applying a logi¬ 
cal 1 to the S input. It can be placed in the R state by 
applying a logical 1 to the R input. Applying a logical 
state change from 0 to 1 to the T input will cause the 
FFto change to the opposite state. The flip-flop R and 
S inputs are direct-coupled. Each can be a-c coupled 
by inserting a 200 pF capacitor in series with the input 
so that these inputs will be activated by only logic state 
changes from 0 to 1 and not by logic levels. 

Internally, the application of a logical 1 to the S input 
turns off the R output, and a logical 1 at the R input 
turns off the S output. The R and S inputs are level 
sensitive in contrast to the T input which is activated 
by state changes. Consequently, if a logical 1 is ap¬ 
plied to both the R and S inputs simultaneously, both 
the R and S outputs will be a logical 0 , and a state 
change at the T input will have no effect on the state of 
the FF. If both R and S inputs are a logical 1 and then 
both are changed to a logical 0 simultaneously, the 
final state of the FF will be either R or S with nearly 
equal probability. If the change to 0 does not occur 
simultaneously, the input which sees the logical 1 last 
will control the state of the FF. 


If either the R or S input is energized by a logical 1 and 
if a pulse or logical variable state change is applied to 
the T input during this time, the FF will act like a one- 
shot mono-stable multivibrator. That is, it will make 
a transition to the opposite state for about 5 M sec. 

Two master reset terminals are located at the patch¬ 
board for each FF module. One will put all the FFs 
in the R state for a transition from 0 to 1 . The other 
reset is level sensitive. As long as a logical one is 
applied all the FFs are constrained to be in the R state. 

COMPUTER CONTROL LOGIC OUTPUTS 

The analog control logic outputs in the integrator mod¬ 
ule must not be used as an input to any logic element 
except a gate. Therefore, special F-RT and R-H out¬ 
puts are located in the FF and gate modules. These 
can be used as an input to all of the logical elements. 
They are identified as F-RT T and R-H T . The logic 
levels for both F-RT T and R-H f are consistent with the 
digital logic elements. There is a logical inversion 
between F-RT and F-RT T as well as R-H and R-H f . 
For example, when F-RT f is connected to a flip-flop 
R input, the flip-flop will be in the reset state when 
The F-RT bus is energized. When the F-RT bus is 
de-energized the FF is free to respond to signals at 
its T and S inputs. 


The F-RT and R-H busses are logical complements, 
which is the reason for their selection in logic control. 

LOGIC SWITCH 

There is a manual logic switch in the counter module 
which is convenient for single-stepping programs. The 
output of this logic switch (ungrounded terminal) can¬ 
not be connected to any other logic output. The reason 
for this is because the output of the logic switch is 
shunted internally by a large capacitor. Consequently, 
if this element is connected to another logical output 
it will cause degradation of the rise-time of that ele¬ 
ment. 
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CHAPTER 9 

CIRCUITS FOR SIMPLE LOGICAL FUNCTIONS 


Several fundamental circuits which are used frequently 
for logical programming are presented below. In all 
cases the trigger is a negative-going voltage. 

ONE-SHOT 


This is due to the 50 /zsec delay at input 1 of the gate. 
Its behavior for multiple triggering is the same as for 
the 5 nsec one-shot. 

Variable: 


A one-shot circuit will produce a logical pulse of pre¬ 
determined length whenever its input is triggered. 
That is, before the trigger, its output will be a logical 
zero. When the trigger occurs its output changes to 
a logical one for a fixed time interval. A one-shot can 
be constructed with logic elements in several ways. 

5 jisec: 

The circuit below will provide a 5 jusec logical pulse 
when triggered. 






If a second trigger occurs during the pulse, then the 
pulse will last for 5 pisec after the second trigger. In 
other words, after the one-shot is turned on by the 
trigger, it will stay on until 5 nsec afterthelast pulse. 
The sensitivity to multiple triggering can be eliminated 
with a gate: 



The logical pulse from this circuit will be delayed 50 
nsec due to the delayed input of the gate used for L^. 

50 /isec: 

The circuit below will provide a 50 ^sec pulse when 
triggered. 



R 



This one-shot will produce a logical pulse of fixed 
duration. The pulse duration depends on the choice 
of R, C. Its behavior for multiple triggering is some¬ 
times erratic. The sensitivity to multiple triggering 
can be eliminated by gating. 



This one-shot will produce a logical pulse whose dur¬ 
ation is determined by both the clock frequency and the 
output line number from the counter. The pulse dura¬ 
tion for the circuit shown is n times the clock period. 
Obviously several dividers may be inserted in series 
with the clock if desired and several counters may be 
used. This circuit is insensitive to multiple triggering. 


HALF-ADDER 


A half-adder is a circuit for implementing an exclusive 
sum, © : 
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Although digital inverters are shown in this circuit, 
they are not usually necessary because all logical ele¬ 
ments except for the counters have both complemented 
and non-complemented outputs. 

SHIFT REGISTER 

A shift register is a binary register which is shifted 
one bit each time a trigger is applied. The triggered 
half-adder is needed for this circuit shown at right. 

If Li0 L2 = 1, then whenever a trigger, T, is ap¬ 
plied, L.Q will be a logical pulse whose length is equal 
to the duration of T. Note that T is applied to input 3. 

If the duration is less than 50 jusec then changes of Lj, 
L.2 after the application of T will have no influence on 
Lq due to the built-in delay for gate inputs 1 and 2. 



The circuit for a shift register is 



Notice that the inverters for the triggered half-adder 
are not required since complemented inputs are avail¬ 
able. An initial setting for this register can be entered 
by resetting the flip-flops individually, or with the com¬ 
mon reset, and subsequently gating a logical one into 
each S input as required. This circuitry is not shown 
because it requires little explanation. Whenever a 
trigger, T, is applied which is less than 50 /xsec (a 5 
Atsec one-shot is a convenient source) the following 
state transfers will be made 


Thus, when each trigger occurs, the current value of 
L i is entered into the first flip-flop. The last flip-flop 
can be connected to the first flip-flop, instead of L^, 
to obtain a register with an end-around shift. 

UP-COUNTER (BINARY) 


The circuit for the binary up-counter is (shown for 
three bits) 


L x —FF1 
FF 1 — FF2 


FF . -*FFn 
n-1 


FFI 


L|- 




FF2 


vd/ 


FF3 
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The counter is initialized with the common reset. Al¬ 
ternatively, an initial value can be entered in the same 
manner as for the shift register. This circuit is acti¬ 
vated by a logically positive change of state. It will 
count the number of positive changes of state in the 
binary system. Each time FF n-1 returns to the R 
state, the state of FF n is changed. The table below 
shows the states of the FF T s for successive logically 
positive changes of state of L. 

Number of 

FF NUMBER State Changes 

3 2 1 for Lj 


FF NUMBER 
3 2 1 


0 1 1 
1 0 0 
1 0 1 
110 
1 1 1 
0 0 0 


Number of 
State Changes 
for 

3 

4 

5 

6 

7 

8 


0 0 0 0 
0 0 1 1 
0 10 2 


DOWN-COUNTER (BINARY) 


The circuit for the binary down-counter is (shown for 
three bits) 


FF I FF2 FF3 



This circuit is activated in the same way as the up- 
counter. The table below shows how it counts. (It is 
assumed that the initial value is the binary number 111.) 


FF NUMBER INPUT 

3 2 1 STATE CHANGES 


1 1 1 
1 1 0 
10 1 
10 0 
0 1 1 
0 10 
0 0 1 
0 0 0 
1 1 1 


0 

1 

2 

3 

4 

5 

6 

7 

8 
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TIME INTERVAL MEASUREMENTS USING DIGITAL 
LOGIC 


The digital logic modules, which are part of the Systron- 
Donner Analog Computer complement, can be used to 
obtain measurements of time intervals with a resolution 
of up to one millisecond. The overall timing accuracy 
depends on the number of digits displayed. It can be 
as good as . 01%, which is a significant improvement 
over the typical 3% accuracy obtained through the use 
of an oscilloscope. Also, this particular capability 
of using digital logic modules to make accurate time 
interval measurements can eliminate the need for a 
frequency counter, which is generally an extra cost 
accessory. 

The basic time interval measuring program is given 
in Figure 1A, and the timing relationships are shown 
in Figure 2. This program measures the duration of 
a logic signal, L. Initially, all flip-flops and DCUs 
(Decimal Counting Units) are in the Reset state. When, 
and only when L changes from 0 to 1, Flip-Flop 1 is 
toggled into the Set condition. Gate 3 is opened and 
the clock timing signals are sent to the Time Interval 
Indicator. This is a cascade of DCUs of sufficient 
quantity to measure the total time interval. At the end 
of the time_interval being measured, L returns to the 
0 state. L from the inverted output of Gate 1 resets 
F-Fl. F-Fl, in turn, toggles F-F2 into the Set state 
and turns off Gate 1. The Time Interval DCUs stop 
counting and hold the elapsed time display. The Set 
output of F- F2 inhibits the output of Gate 2 through 
Diode 1, preventing any further changes in L from 
toggling F-Fl. F-F2 f s Set output also opens Gate 4, 
and clock signals are sent into the Display Time DCUs. 

The Display Time Generator is also a cascade of DCUs. 
Outputs from the Display Time Generator are logically 
summed in Gate 5. When the count reaches the value 
corresponding to the desired display time, Gate 5 is 
changed to the 1 state and toggles F-F3. F-F3 is 
connected as a 5 microsecond one-shot. That is, F-F3 
goes into the Set state for 5 microseconds before 
returning to the Reset state. The Set output of F-F3 
resets all flip-flops and DCUs, and the circuit is ready 
to measure another time interval. The automatic 
Display Time Generator can be dispensed with, and 
the system can be manually reset by using the logic 
switch on the clock module (Model 3328). 

If L changes state before the end of the display time, 
nothing happens due to the Inhibiting action of F-F2 and 
Dl. If the reset action at the end of the display time 
occurs while L is in the 1 state, the system remains 
in Reset. Interval measurement will not commence 
again until L changes from 0 to 1. This is illustrated 
in the Timing Sequence diagram, 3rd interval time. 

To summarize, this program takes no action until L 
changes from 0 to 1. It then measures the current 
time interval, holds the display until automatically or 
manually reset, and then waits for the next state change 
from 0 to 1 in L. 

An alternate program is shown in Figure IB for meas¬ 
uring the time interval between two independent logic 


signals. L a initiates the timing action. Gate 2 acts 
as a buffer so that L a is not affected by the inhibiting 
action through Dl. terminates the time interval 
when it goes into the 1 state. The interval terminates 
at that time independently of the state of L a . If 
causes trouble by being in the 1 state during part of the 
interval to be timed, it can be connected to a one-shot 
similar to F-F3 and then into the Reset input of F-Fl. 

The program of Figure 3 allows measuring of the Nth 
time interval after reset. It is identical to the program 
of Figure 1 except that a pre-set event counter is in¬ 
serted between L and the toggle input of F-Fl. The 
event counter can be made of DCUs, or a cascade of 
flip-flops. 

The period of N intervals can also be_measured, by 
inserting the pre-set circuit between L and the Reset 
input of F-Fl. This is useful for measuring equal- 
period intervals with an effective resolution of 1/N 
milliseconds, or for averaging readings over several 
cycles. 

Figure 4 shows how to convert analog signals to logic 
levels. The logic gates can be operated from analog 
inputs; the output will conform to logic states. How¬ 
ever, the gates have hysteresis. That is, the gate is 
off for inputs greater than about + 8 volts. Once opened, 
the gate will not turn on until the input becomes less 
than about+4 volts. Negative inputs to the gates cause 
no action. Also, the gates are not damaged by input 
voltages within the ±100 volt computing range. 

Gate B is turned on at the start of the time interval 
period, thereby toggling F-F A into the Set state. 
Gate B turns off during the interval, but this has no 
effect on F- F A. When Gate B turns on again, F- F A 
resumes its Reset state and the timing is terminated. 

Figure 5 indicates the method used to convert inte¬ 
grator mode logic bus signals to logic levels. Gate 
input 1 must be grounded or connected to a logical 1. 
When the mode bus is energized, it turns the gate off. 
L then is taken from the inverted output of the gate 
since this is a logical 1 when the mode logic bus is 
energized. 

Display of the less significant digits of a time interval 
(when the more significant digits are known), can 
appear on the Model 3328 which allows only a 3-digit 
display. A 3-digit display of the less significant digits 
will suffice when the more significant digits are known. 
For example, to read an 11 second interval to a reso¬ 
lution of 1 millisecond requires a 5-digit display. If 
the first two digits are already known, the last 3 can 
be indicated on the clock. The counters are cyclic, 
so overflows of the more significant-digit counters do 
not affect the count in the less significant digits. 

As another example, the period of a 1 radian per second 
sine wave can be measured to 1 millisecond resolution, 
and. 015% accuracy, without displaying the most signif¬ 
icant digit since the period is known to be 6. + seconds. 
The display would indicate the . 283 part of the 6.283 
second period. 
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Figure 1A. Time Interval Measurement Program 
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Figure IB. Alternate for Interval Between Two Logical Events 
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Figure 2. Timing Sequence 



Figure 3. Alternate To Measure Nth Interval After Reset 
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Figure 4. Program to Convert Analog Signal to Logic Signal 


(LOGICAL I) 


INTEGRATOR 
MODE LOGIC 
BUS ( + 28 V) 



L 


Figure 5. Program to Convert Integrator Mode Logic (+28V) Bus to Logic Signal 
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CHAPTER 10 

CIRCUITS FOR SIMPLE LINEAR AND NON-LINEAR FUNCTIONS 


Several simple operations are used frequently for pro¬ 
gramming. They include division by a constant, dif¬ 
ferentiation, time delay, and a variety of sectionally 
linear transfer functions. These are described below. 
(It should be remembered that all electronic switches 
in this section must be used as J inputs in series with 
the input resistor.) 

DIVISION BY A CONSTANT 

The circuit 


DIFFERENTIATION 

Programs are usually written to avoid differentiation. 
This is done because differentiation inherently stresses 
noise and because the circuits used may lead to insta¬ 
bility. However, in some instances it is difficult to 
avoid forming derivatives. This is especially true for 
on-line data processing or editing applications. For 
this reason, it is helpful to be familiar with available 
differentiation techniques. 

A single amplifier program which can be used is 


IM 



The symbol 




which has the transfer function 


e 

— = RCs 
e l 


1 

rCs + 1 


r is usually about 10K producing reduced gain at fre¬ 
quencies above (2TrrC)"^ cps. Since C is generally 
large, the circuit tends to load unduly the preceding 
computer element. Also, a large C may lead to sta¬ 
bility problems in the amplifier driving the differenti¬ 
ator. The small resistor, r, helps prevent drive cir¬ 
cuit instability. 


The more sophisticated program 


satisfies the condition 

cte Q + e 1 + 10e 2 = 0. 

(The input gain labels are usually omitted if they are 
unity.) The transfer function is 

1 10 

e 0 a e l a e 2 • 


This circuit is useful for dividing a sum of variables 
by the same constant, a . Sometimes it is more con¬ 
venient to set a pot to a, if a is a natural parameter 
in the program. Otherwise the pot would be put between 
the amplifier output and the following input and would 
require the calculation of aT^ for each new setting. 

This configuration is also useful for obtaining gains of 
greater than 10 without requiring more input resistors 
or more amplifiers. 



has the transfer function 

fo _ S 1 + Ts 
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where T = CRin, Rin> the input resistance at ej, is 
usually 1 Mfi. Thus, expressing C in microfarads 
T (sec.) = C(/if). 


a is large but not close to unity. This is due to the 
stability of the loop composed of amplifiers 1 and 2, 
and pot 3. A good choice for C is 


Ideally, 


a= 1 - K^C 


but in practice j a is adjusted to a slightly different value 
due to amplifier and potentiometer phase shift. The 
error of the circuit is strongly dependent on a causing 
its setting to be critical. C must be selected so that 


C 


0.3 

* 1*2 


The circuit will give good performance for 10 < K1K2 
< 1000. The lower limit is necessary in order to limit 
the size of C and thus avoid stability problems. It fol¬ 
lows that the setting for a is roughly 0. 7. 

The test circuit for adjusting a is (the block denotes 
the differentiator circuit), 


DIFFERENTIATOR 



ERROR 


The oscillator is set to the highest frequency of interest 
and is adjusted until the error is minimum, using a 
scope. In practice a>0. 7. 

An excellent low-frequency differentiation circuit is 


This circuit works well for 


1 



< 100. If desired, 



pot 1 can be adjusted precisely with the previous test 
circuit, in which the integrator input is set to 


*2 . 



TIME DELAY (Continuous) 

Sometimes a program for time delay is needed for 
solving a problem. The time-domain definition of time 
delay is 


It has the transfer function 


e Q (t) = e x (t - T). 


— = -K.s 
e l 1 


. Is + 1 


001s^ + .01s + 


and at low frequency 


This has the transfer function 


e 


-sT 


K i 


*2 


s . 


There is no convenient way of generating a continuous 
time delay with only analog computer elements. How- 
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ever, various approximations can be made; most of which are due to Pade. The simplest (first order) is 


1 


1 


+ 



which works reasonably well for a delay of <M0. 6 where 4> is expressed in radians. Since 4> = «T, where w is the 
frequency, 

T * ^6 . 

0) 


That is, at very low frequencies, the approximation will provide relatively large time delays with reasonable 
accuracy. The program is 



In practice, the highest significant frequency in the simulation is estimated in order to determine an upper 
bound for T. 

A second-order approximation which is reasonably accurate for 4> -1.8 is 


1 - 


Ts 

2 


1 + 


Ts 

2 


tV 

12 

T2s 2 

12 


The program is 



A fourth-order approximation which is reasonably accurate for <t> - 7. 5 is 

% _ 1 - jW + jsfrs) 2 - M< Ts ) 3 2 l*S8o< T °» 4 

'l l* 5 <Ts) ♦ ^(T«) 2 . jijlTs) 3 * 


10-3 


The program is 



DEAD-ZONE 1 

The dead-zone operation relates e^, e Q by the sectionally linear function 



The analog program is 




+100 



In the following circuits note that mode relays and electronic switches are operated by logical zeroes. The 
circuits have been shown with electronic switches for mode control. Function relays can be used instead of 
electronic switches, if high-speed operation is not required. 
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Pot 1 is set so that its diode will not conduct until e 1 > a. Pot 2 is set so that its diode will not conduct until 
e^ < -b. Thus, for -b < e^ < a, neither diode will conduct and e Q = 0. Also if pot 3 is set so the total gain is m 
then 

e o = ( e ! - a), e x - a 

= 0 , -b< e^ a 

= -m (e^ + b), e^ 2 = -b. 

Note that pot 3 cannot be set to m exactly, since the source impedance due to the pot-diode networks is not zero. 
If digital logic is used, the program is 


ELECTRONIC 

SWITCHES 



me^ 

-mb < me^ < ma 

+, > ma 
< -mb 



U 1 U 2 U 3 

1 0 1 

0 1 1 

0 0 0 


ES x ES 2 ES 3 

o o o 

c c o 

c o c 


e 


o 


0 

-m(e 1 - a) 
-m(e^ +b) 


Pots 2 and 3 are used to set the breakpoints for the 
functions. Their outputs are used both to set the com¬ 
parator references and to provide a biased input to the 
amplifier as required. The amplifier output will be 
zero when all three switches are open. Logical vari¬ 
ables Lp L 2 are defined by 

= sgn (e^ - a) (output of Cl) 

L 2 = sgn (e^ + b) (output of C2) 


ES i is open when L 1 = L 2 = 1, which is equivalent to 


-b - e^ ^ a. 

Since pot 1 is set to m, e 0 will have the desired func¬ 
tional relationship to e*. Normally closed and nor¬ 
mally open contacts of the electronic switches are used 
to provide the logic level inversion required. 
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SATURATION 

The saturation function is 



The analog program is 



Pots 1 and 2 are set to limit at -a and b respectively. This circuit gives a "soft" limit. That is, e n does not 

remain constant at the saturation level, but continues to increase with increasing e x due to the non-zero feedback 
resistance when the diodes are conduct ing 1 «iuieeaDacK 


The program with digital logic is 





The operation of the electronic switches is governed by 
the logical equations 


II 

CO 

W 

W h 2 

(=* 

m 

“ e l 

es 2 = 

r 2 

(e i 

a-) 

m 

ES 3 = 

L 1 

(e i 

■?|s 

VI 


The normally closed or open contacts of the electronic 
switches are used to provide the logic level inversions 
where required. The logic program will provide hard- 
limiting; i. e., the slope after saturation will be zero. 
This is not true for the purely analog program. 

BINARY (bang-bang) 

The binary function is 


«o 



The analog program is 



-100 


Pots 1 and 2 are set to limit at -a and b respectively. 
As in the "soft" saturation function, the output in the 
limiting region has a non-zero slope. 

The program with logic is 
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ABSOLUTE VALUE (variable slope) 


The absolute value function with variable slope is 



The analog program is 



The program with logic is 



BACKLASH 


The backlash function is 
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The analog program is 



Pot 3 is set to m where m is the desired slope. Pots 
1 and 2 are set to limit at -b and a respectively. The 
limiting network and amplifier 3 comprise a dead-zone 
circuit whose output is an integrator rather than a sum¬ 
mer. Integrator drift will be about 50jj.v/sec/p.fd, so 
it is necessary to consider two factors when selecting 
the integrator capacitor: 

1 - Sufficiently small IC time constant to permit 
adequate tracking 


2 - Sufficiently large C to ensure that the inte¬ 
grator drift is permissible. 

The logic program for this function is obtained by re¬ 
placing the analog dead-zone circuit with the equivalent 
logic program. 

BISTABLE FLIP-FLOP 

The bistable flip-flop function is 




The analog program is 



Pot 1 adjusts the value of e^ which results in a change 
of state, namely C. 


DITHER 

The dither function is 
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In both the analog and the logic programs, the gain of integrator 1 determines the dither frequency. 
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CHAPTER 11 

SIMULATION OF CONSTANT COEFFICIENT TRANSFER FUNCTIONS 


The simulation of control systems requires programs 
for constant coefficient transfer functions of the type 


T(s) 


P(s) 


where P, Q, have polynomial form: 


P(s) 


= £ a i sl 

i = 0 


m 


Q(s) £ 


b.s 

0 


where G(s) is the forward path transfer function and 
H(s) is the feedback path transfer function. 


E 2 


Ei 



Eo 


If n > m differentiation is implied. Differentiators are 
treated inChapter 10. Here, it is assumed that n < m. 
For physically realizable systems the aj<, bf are all 
real numbers. It is well known that any polynomial 
with real coefficients can be represented as a product 
of quadratic factors. Further, in many practical 
problems the physical model is composed of dominant 
second or first order transfer functions. For these 
reasons only cases for which m < 2 are considered. 
Programs for specific higher order transfer functions 
can be derived easily with the approach outlined in 
Chapter 3. The materialhere is intended to comprise 
a dictionary for second or first order programs. With 
this information, generating a program for a complex 
system diagram with many blocks is easy. The pro¬ 
gram is written for each block and the system is rep¬ 
resented by making appropriate interconnections be¬ 
tween blocks. 


The transfer functions of this system are 



G 1 °2 


HGi G 2 


E, 


HGi G 2 



That is, the response to inputs other than Ej is equal 
to the normal response divided by the transfer function 
preceding the actual input. 

Specific transfer function programs are given below. 


In developing transfer functions using only pots, am¬ 
plifiers, and integrators, the following fundamental 
characteristics of feedback or closed-loop systems 
are of interest. 

General Closed-Loop System 



SUMMING POINT SYMBOL 



The transfer function of this system is 



G 

1 - HG 




a 
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For a, p > 0, a 2 > 4 P, another program is preferable. 
Let 6=f, \ 2 = S 2 -/?. Then 

y 

= (s + 6 + a) (s + <5 -A) 
and the program is 



Note that the condition a,/3 > 0 ensures (<$+A), (8-A)>0. 


T(s) - 2 + *■- 

s + as + p 


The program for this function uses types (1), (2), & 

y 

(3) of T(s) =—g-, shown above, with the 

S + aS + /3 

addition of two pots and one or two summers. In each 
of (1), (2), (3), pot 1 is set with 7=1 
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CHAPTER 12 

CONTROL SYSTEM SIMULATION 


A model of a control system is usually represented by 
a block diagram The block diagram contains constant 
coefficient transfer functions, non-linearities, and 
frequently, special differential equations. The easiest 
way to generate the complete computer program for 
the model is to 


1 - derive a program for each block 

2 - interconnect the blocks. 

This method is called block programming and is simple 
using information developed in Chapter 4. Some illus¬ 
trative examples are shown below. 


MODEL: 



s + a 


1 . 
s + 3 * 



<D=j 
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MODEL: 







CIRCUIT 


PROGRAM: 


R 



Note that no pot is required for the parameter, B, be¬ 
cause this parameter is used to translate the dead- 
zone function. 
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CHAPTER 13 

FUNDAMENTALS OF VECTOR ANALYSIS 


The development of vector operators given below is 
within the framework of Euclidean (x, y, or x, y, z) 
space. Due to the invariance property of vectors and 
vector operators under transformation to other coor¬ 
dinate systems of engineering interest, it is reasonable 
to proceed in this manner. Sometimes we will show 
the representation of the vector operator under discus¬ 
sion in other coordinate systems for the purpose of 
illustration. Indeed, the main motivation for using 
vector notation is its non-dependence on a particular 
coordinate system. 


For a vector V = Vx, Vy, Vz; 


a 3V 

V x V = —- 

9 y 


d V 
~d 


dV 


aV dV 
_z _y 


aV 


z z xx 

is a vector-valued function known as curl V. 


_ _ v 2 ai £_ 

= 9 x 2 9 y 2 9 Z 2 


is known as the 


DEFINITIONS 

(A • B) is the scalar (dot) product of two vectors. 

A • B = |A|. |B|Cos (A, B) 

(A x B) is the vector (cross) product of two vectors. 


If 


C = A x B 


then 

|c| = |a|- |b| Sin (A, B). 

Note that 

(C x (A x B)) = 0 

Vis called the vector operator del, 

A _a_ d_ d_ 

dx ay’ az 

(Here ordered n-tuples are used instead of the cus¬ 
tomary i i k unit vector notation, i. e. = ■ ~ i 

7 ax , ay ax 

a .v 

If 0 = <p (x, y, z) is a scalar-valued function, then 

v . _ d< P d< P ; 90 

^ ax' ay az 

is a vector-valued function also known as grad <f>. 

Vhas different forms in other coordinate systems. 
Example: In cylindrical coordinates 

ar ? rae’ az 


Laplacian. 

SUMMARY OF OPERATOR ALGEBRA (greek letters 
are scalar-valued functions, capital letters are vec¬ 
tor-valued functions) 

V* (aV) = a (V • V) + V * (Va) 

V x (aV) = a(vxv) + (v a ) xV 

V • (VxU) - U * (VxV) - V • (VxU) 

V(V • U) = V x (VxU) + (V • V) U + Ux (WV) + 
(U • V) V 

Vx (V x U) = (u • V) V - (V • V) U - U (V • V) + 

v(v- u) 

V X (va) = o 

V . (v x V) = o 

Vx (Vx A) = V (v • A) - V 2 A (this defines V 2 
which is not a true vector operation). 

AxB =- BxA 

A • (B x C) = (A x B) • C (triple scalar product) 
DIFFERENTIAL OPERATIONS 

The remainder of what follows will be derived for two 
dimensions but may be readily extended to three. 

Divergence: 

If f (x, y) is a scalar-valued function then the family 
f = c of curves in the plane are level lines for f. If 
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f is thought of as a potential then these curves are 
equipotential lines. 


outfiux v x (x+h , y +£ ) - V x > y + ? ’> 

unit area h 


df = 


af_ 

ax 


dx + 


af 

ay 


dy 


= o 


Vy(x+£, y + k) - V y (x + | T , y) 


Vf • (dx, dy) = 0 

where (dx, dy) is a vector element of arc along the 
curve f = c. (This is so because we constrain the 
values of x, y so that f = c). Thus we see that 

Vf _L (dx, dy) , 


In the limit as h, k—►O we have 


Outfiux _ 9V x + dV y 
unit area ax + ay 


V • V = div V 


and the family of curves obtained from the differential 
tangent field, Vf, are orthogonal to the family of level 

A df 

lines and represent a flux (or flow), vf • N = ^ is 

called the directional derivative of f in the N direction 
(N is a unit vector) at the point (x, y). The meaning 


is obvious. V • V 


av av 

= —- + ^ is a scalar-valued 

a a 
x y 

function known as the divergence of the vector-valued 
function,V. It is commonly written as div V. 


Consider a small rectangle in the x, y-plane 


j’ -1 | (x + h,y+k) 

(x.y) 


Thus if <4 is a conservative potential field 

V • V <p = 0 or V 2 <*> = 0 

2 source • 

At a point x, y; V <p ^ 0 implies a local sin ^ in the 

case of incompressible flow. Note that with compres¬ 
sible steady flow V 2 $ * 0, in general, without the 
presence of sources or sinks. In transient physical 
processes v 2 <p is a measure of the change of energy 
from one form to another. For example, in heat trans¬ 
fer processes heat is stored locally by the distributed 
heat capacity of the medium in which the process occurs. 
Thus, if T is the temperature at (x, y), the local con¬ 
ductive heat transfer is k VT where k is the thermal 
conductivity, and the local heat stored by the heat ca¬ 
pacity of the medium is cAT where c is distributed 
heat capacity. The loss of thermal flux (heat flow) 
from the conductive process equals the heat flux re¬ 
quired for heat storage and this time-varying relation¬ 
ship is given by 


We think of V as a flux, and determine the net emer¬ 
gence of flux from the rectangle, 


Vy(x+£,y +k) 


-Vx(x, 



We have 

Outflow = 


V x (x+h, y + f) - V x (x, y + f 


+ R 


1 


V y (x+£, y+k)- V y (x + £ f , y) 


k 

h 


V kvT = c f 


or 


v 2 t 


£ 53 ! 
k at 


^ It (Diffusion 
equation). 


where K is the thermal diffusivity constant. 

In the case of the small oscillations ofauniform mem¬ 
brane, if w(x, y) is defined to be deflection, then the 
potential energy stored by deformation of the membrane 
is a function of Vw while the local kinetic energy of 
motion due to the mass of the membrane is a function 
0 w 

of TT-. Thus the time-varying transfer of energy from 
ot 

one state to the other is described by 


2 1 0w 

V w = —2 — 2~ (wave equation) 

c at 


where c is the wave propagation velocity. 

Curl: 

If we think of the vector-valued function, V, as a ve¬ 
locity field, we can interpret curl V in a convenient 
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way. Consider the tangential velocity around the sides righthand rule, the shaft then points in the direction 
of a small rectangle. of VxV and its angular velocity is proportional to the 

magnitude of VxV. 


- Vx (x+£,y+k) 

(x+h, y+ k) 

Vy (x + h,y+£) 

Vx(x + £'y) 


-Vy(x,y+£') 


( v 



The (spin) around the boundary is given by 


Vy( x + h> y + £) - V y (x, y + £') 


V x (x + | T , y) - V x (x + | , y + k) 


If VxV = 0, V is said to be an irrotational field. 

If V = VxA, V is said to be a solenoidal field. 

If V = V* (pthen V is irrotational. (VxV- <f> = 0). 
If V = VxAthenVhas zero divergence (v- VxA = 0). 

VECTOR TRANSFORMATION THEOREMS 


Theor 1: (Divergence or Gauss theorem) 

Let R be a region over which a vector valued function 
v is defined, then (differentiability assumed) 



where n is a unit exterior normal to aR (boundary of R). 
Theor 2: (Stokes' theorem) 


In the limit as h, k—►O we have 


circulation, ' _ d _\ VxV=curlv 

unit area ^ puv ax ay vx v curi v ‘ 


To visualize the result in three-dimensions, we think 
of an infinitesimal paddle wheel located at the point 
(x, y, z). This paddle wheel is angularly oriented so 
that its shaft has the greatest rate of spin. Using the 


Let R be a 3-dimensional region over which a vector 
valued function V is defined, and let v be a 2-dimen¬ 
sional subregion of R, then (differentiability assumed) 

f N • (VxV) = y V • ds 

V 8(T 


where N is a unit positive (sense to be determined by 9 ct) 
normal to a and ds is a direction vector which lies in da. 
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CHAPTER 14 

PARTIAL DIFFERENTIAL EQUATIONS - PART 1 


This section treats the class of partial differential 
equations which can be represented 

V- KV<t> = E<fc <t> = <t>(x v x n , t) 


where the operation, E, is defined by one of 


a<f> = 0 

(Laplace equation) 

£ 

II 

(Poisson equation) 

E0 = k0 

(Helmholtz equation) 

* ■ <=<*>!? 

(Diffusion equation) 

*</> = M(X)^ 

(Wave equation) 


The analog computer is used to integrate E0, in other 

words to generate E *. This will be discussed before 
considering the generation of V • K V0. The block pro¬ 
grams for these cases are 


Laplace: 



The high-gain amplifier forces E0 = 0 when connected 
to the remainder of the program. The capacitor is used 
to reduce high-frequency gain which may result in 
instability. 

Poisson: 


<*> 



The high-gain amplifier ensures E0 - k = 0. 


Helmholtz: 



The capacitor is added to prevent algebraic loop insta¬ 
bility. 

Diffusion: (N is the time scale factor) 


4 > 



Wave: (N is the time scale factor) 





With the blocks for E0 available, all that remains is 
to find how to represent V • KV0; the Laplacian of 0. 
The form of V will depend on the coordinate system 
(see Chapter 13). To illustrate the concept we will 
derive the Laplacian for one dimension first and then 
develop it for both two-dimensional Euclidean and Rei- 
mannian manifolds. The extension to three dimensions 
will be obvious. For simplicity it will be assumed that 
K is not a function of time. For the case where K is a 
function of time, the pot set to K is replaced by a multi¬ 
plier. 
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ONE DIMENSION 


V- KV0 = 


ax 


[r( X ) . 


Using the same approximation 


[v • KV^>], 


[KV0] x x -[KV0] x 1 
n+n 


n-r 


The derivatives will be replaced by finite differences: 


x 1 - x 1 
n+g n-g 


= S 6 


Since 


a_ 

ax 


f(x,t) = lim —- +AX?t | x ~ f - j — 
Ax —0 ^ 


Therefore, by substituting the different equations pre¬ 
viously obtained, we have 


then, for small Ax 


a_ . Ax .» f(x+ Ax, t) - f(x, t) 

ax 1 2 ’ ' ax 


Adopting the notations 



+ 1 


+ x 


and 


K ».l [V 1 - *n] K n-i [*n-Vl] 

x „.r x „ ' 


= /x 1 - X 1 \ S 4 
{ n+ 2 n ~2/ 

The program which generates Etf> n is 


f(x n , t) = f n for any n 


we have, for finite differences 

AX 


ffXg, t) - f(x r t) 


x 2 ‘ X 1 


For any n 

-(^4^) 4m-‘) 


AX 


Af 1 
n+ 2 
AX 


f n + l - f n 
x n+l ' x n 


An 


f(x n + r 4) - f(x n> l) 


X . - X 

n +1 n 


Thus 


0(x n+1 t) - <p(x. t) 

[KV^] 1 = K(x 4 ) • .— ’ - 

L x n4 n+ 2 x n+l- x n 


= K 1 

n+ I 


Tl+l 


- 


X - - X 

n +1 n 


Similarly 


[KV*] 




0 n ' Vl 


X - X 1 
n n -1 


-</>„-! +<#>0 n+1 



If x n+ ^ - x r -* Ax, for all n, then the program reduces to 
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ONE - DIMENSIONAL BOUNDARY CONDITIONS 


Boundary conditions of two types are commonly en¬ 
countered: 


4*. t) = f(t) [Type 1] 


Kfa) Vtf> (a, t) = f(t) [Type 2] 

Sometimes both types are specified simultaneously as 
a mixed boundary condition. These are approximated 
by 


f(t), x n = a 


K-J. v* n 1 - lit), x 1 - « 


n+« 


The programs are shown below for the Wave equation, 
but are typical for all S <f>. 


Type 1: 


-*n 
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Type 2: 


TWO-DIMENSIONAL EUCLIDEAN MANIFOLD 



Ax 


Mixed boundary conditions are represented by 



f 2 (t) 


AX 


This discussion can be extended to higher-dimension 
cases by treating each of the boundary points in the 
same fashion. 


A two-dimensional Euclidean manifold is commonly 
called an x, y-coordinate system. Here 

v - Kv * = (fc %) • ( K(x > y) H K(x ’ y) If) 

and 


V ' KV ^= fe K ff + 


K — = 

ay ay 




As for the one-dimensional case, finite difference 
approximations are made for the derivatives: we will 
adopt the notation 

^n, m = ^ X n ? 




H Vl * x n = Ax ’ y m + l - y m = Ay ’ for a11 n ’ m > 

then the equation reduces to 


Vim 


V,. 


2 ’ / \ 2 ’ / \ 

Ax yn+1, m" ^n, mj + Ax yn-1, m ^n, mj 


K 1 
n, m +2 

Ay 


K 1 
n, m-o 


rn, m+1 n, ml Ay I n, m-1 n,mj 


Ax AyE# 


n, m. 
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The program is 



TWO-DIMENSIONAL RIEMANN MANIFOLD Note that 


Let the metric for the manifold be given by (ds * arc 
length) 

ds^ = gij dx*dx^ . 

In what follows we will use summation notation (i. e., we 
sum on repeated indices). 

g = |gij| = g^G^ (no sum on i) 

so that G 1J is the cofactor of g.. in the determinant I g.. I. 

ij 



It can be shown, that if <p is a scalar function, 


v 2 * <v£V J -2t) 

Vg Vox / ax 


Thus, we have a general expression for V 0 in non- 
orthogonal curvilinear coordinates. We could treat the 
non-orthogonal case. However, that treatment is much 
more complicated than for the orthogonal coordinates. 
Since the latter is of greater practical interest, we will 
restrict our attention to 


We define 



= °> 1 * i> 

g lj = 0, i * j 
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Thus, in three dimensions, 


When approximated with finite differences 


. 2 l x 2 ,, 2,2 (A 3v2 

ds = g^dx ) + g 22 (dx ) + g 33 ( dx ) 


8 ~ g ll g 22 g 33 


g** = — (no sum on i) 
g ii 

We denote 

0 (x*, x^, x k ); x a = A^,a=i, j,k, are constant 

by *u,*; 

Consider the diffusion equation 
V • K V 0 = 2J0 

We will need to replace the space derivatives with 
finite differences. 


_ 9 _ / K V^2^3 dA 

ax 1 y VgJJ 9xy 


^33 d<t>\ ^ 

1I 2 3 

dx / Ax ax Ax 


2 3 

Ax Ax 


Ax 


/Kvg^v£r 3 \ 

V v^i / 


i+1/2, k, j J 


- 

r 

Vl, j,k"^i, j, k 

+ 


2 3 

AX AX 


AX 


/ Kv ^22 Vg 33 N \ 

V viTi ) 


2 > j) ^ 


^i-1, j, k~^i, j, k 


V . K V<*> =—— (K Vg g 1 ^ ) 

vg" 3x ax J 


= E<fc K=K(x\ X 2 , x 3 ) 


and similarly for the other terms. Let us consider 
the one-dimensional projection to determine the sig¬ 
nificance of these terms. From the definition of the 
metric, clearly 


Now 


ax' 


r (K Vg g ij ) = Vg H0 


g lj = 0, i*j 

11 Vg Vg 22 g 33 


V^g 


g n Vg^ 


Vgg' 


22 Vg ll g 33 


Ax v/g n 


Ax 2 vi^ 


AX Vg^ 

12 3 

are lengths in the x , x , x directions respectively. 

12 3 12 

We denote these by 1 ,1 ,1 respectively. Ax Ax 

Ax^Vg” is the volume V = l 1 !^!^. 


and 


Vgg' 




33 Vg ll g 22 


/ Kv ^ 22 v/ %3 a<A 

v “7 


In general, the value of the lengths and the volume for 
12 3 

increments in x , x , x depends on location. The one- 
dimensional projection of the problem is 


[Hr) )] 


. _3_ / 

( KVf^VgJg 

aA 

, I7iVk\ 1 

H'l’j’k'Hi.k)] 

ax 2 ( 

V ^2 

ax 2 / 

LI iW 1-j.j.k ' 

+ 1 

ax ' 

^ ^^33 


V i, j,k * ^i, j, k 
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<f> DEFINED HERE 



Thus, for example, we see that, in the Diffusion 
equation 



and 


i _ 2> i> k 


i 2 i 3 k\ 

l 1 / i+1/2, j,k 


are lumped conductivities defined at appropriate loca¬ 
tions and V. . , S is a lumped capacity, also defined at 
if ^ 

an appropriate location. The extension to the non-pro- 
jected problem in three-dimens ions is obvious. It is 
interesting to note that the theoretical result agrees 
with our intuition. 

ASYMMETRIC LATTICE 

The set of space points over which the operator 


be found such that the perpendiculars in each polygon 
enclosed by a branch loop intersect at only one point. 



V • KV<p 

is defined is called a lattice. In arbitrary regions one 
cannot find a metric in any simple way. Thus, the ap¬ 
proach developed for the Riemann manifold is npt ap¬ 
plicable. If a Euclidean system is used, then the bound¬ 
ary must be composed of a collection of orthogonal line 
segments. This has two disadvantages: the lattice 
density must be high on the boundary in order to repre¬ 
sent the boundary. The lattice density on the boundary 
determines the lattice density in the interior of the 
region, leading to an exorbitant equipment requirement 
for the program. The asymmetric lattice representa¬ 
tion of V • KVp will allow one to avoid both of these 
difficulties. 

Consider a lattice whose only restriction is that a set 
of perpendiculars to the branches (one for each) can 


Referring to the figure above, 

1. Let a point function p be defined over the 

lattice a 

2. Consider one point of the lattice where we 
arbitrarily choose <£q and denote by <f>& the s 

at branch connected neighbors 

3. Define f. to be the length of the branch between 
lattice points associated with <£q, 

4. Define r. to be the length of the _L to i . 

r. 

5. Define the symmetry ratio Y. = — 
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Now consider the representation of V** 0 at the lattice 
point 0 given by 

°<»o> - E xr <*, - +0> 

a 0 

Where a is the index set of the neighbors. We will 
choose the Y. and A such as to minimize error. Now 
for the contiguous medium (where we denote the loca¬ 
tion of lattice point i by S.)each 0(S.) is given in terms 

of 0 (S Q ) and the derivatives of 0 at S q by 


Define 


^x 2 £ Y i x i ’ ^y 2^ Y i^i 9 ^xy S Y i x i y i 

a a a 


Then 0 can be expressed 


0(d) ) = ft 0 + R 0 +6 0 

xv o p x ^xxo p y yyo p xy xyo 


+ higher order terms. 


♦@i> ■ * s „> * £ JT <*i I? -nlr’" 

n=l 

We have assigned a Cartesian metric to the domain so 
that f. is the pair (x., y.). Thus the representation 

O( 0 Q ) is given by 


We want to reduce this to canonical form. The result 
above is independent of an orthogonal transformation 
on (x, y). Note that there exists an orthogonal trans¬ 
formation such that p =0. Henceforth, without loss 
xy 

of generality, we will assume ft = 0 . 

xy 

Thus, 


£y £ — (x — + v —) n 

“ n! K i ax y i dy 1 

a n=l 


0(0 ) = b 0 + 6 0 + higher order terms. 

p x ^xxo p y yyo & 


2 ir «*i k * 

n=l a 


0(d) ) = ((d> + 4> ) 

vv o ; v 2 ' xxo yyo 


p x +p v 

+ ( 2 ) (0 xxo ' <?> yyo ) + higher order terms. 


0(0 ) = Y* Y.x. 0 + y Y.y. 

VY o ^11 ^xo ^ ri 


*. T ‘TEVi 1 ♦, 


+ 23 Y^.y. 0 xyQ + higher order terms 


We see that if p x is significantly different in magnitude 

from py, then 0 is a poor approximation of . Thus 

we require the polygon formed by the perpendiculars 
to have relatively good symmetry properties. The 
higher order terms can be reduced by increasing the 
lattice density. One can get an intuitive feeling for the 
symmetry properties of the polygon by observing that 
the system has a physical interpretation in terms of 
moments of inertia as shown below. 


Now we require 

£Vi = °’ £ Y i y i = 0 (1) 

a a 

( 1 ), above, is the x and y projection of £ Y-* • • 

a 

EV, ■ E r i 

a a: 


projected on x, y is zero since the polygon formed. 

a 

by the r^ is closed. Thus, (1) is automatically satis¬ 
fied. Evidently our choice of Y^ was suitable. 
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Interpretation: 


All that remains is to find the normalizing constant 
A . We choose 

0 ^ i ij is a branch from 0. to </>, then the lumped flux 

a _ + p y from polygon i to polygon j corresponding to Vtf> is 

o 2 r ii r 1 r * 

f K ’ S0 Y ii = 7^ is a lumped admit - 

12212 ijLJ ij 

A q - + y i ) = tance. A. is the area associated with lattice point j. 

a a 


MacNeal, R.H., tT An Asymmetrical Finite Differ¬ 
ence Network", Quarterly Appl Math,October 1953) 

Geometrically, 


V = jyr.t 
o 4^ i l 
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CHAPTER 15 

BASIC ITERATIVE PROGRAMMING 


REPETITIVE OPERATION intervals can be set approximately by the controls at 

the console. Usually the C interval is adjusted so that 
In Repetitive Operation (RO) the entire problem on the the problem solution time is long enough for one to 

computer is solved repetitively. This is useful for observe the significant part of the solution. The R 

determining the influence of parameters on the solu- interval is set so that all integrators return to correct 

tion. Potentiometers and switches are operated man- initial values. The RO controls are calibrated only 

ually while the computer solves the problem repeatedly. roughly. For a more precise setting of the IC inter- 

The RO mode cycle is R, C, R, C, etc. The length vals, the program below can be used with the RO 

of the R (initial condition or reset) and C (compute) selector switch at the console set to EXTERNAL: 



RESET 
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The output from the F-RT' terminal resets the divider, 
counters, and clock, because the logical function 
(GATE 1) + F-RT T is connected to the OV reset ter¬ 
minal. Note that RH=F-RT\ The clock is off as long 
as the F-RT T terminal is a logical 1. The F-RT f ter¬ 
minal is a logical 1 when the main machine mode is R. 
The FF is reset whenever 

(GATE 1) + F-RT T = 1 

When the FF is reset all the integrators will be in the 
R (IC) mode since the R output of the FF is connected: 
through a driver to the R terminal. Each time the FF 
goes to the S state, the divider and counters will be 
reset momentarily due to the connection of the in¬ 
verted output of Gate 2 to the +V terminal. The FF 
will be set when Gate 2 = 1. The logical equations 
for the gates are, 

G1 = m 1 • m 2 - m 3 
G2 = n 2 n 3* 

Thus, the reset and compute intervals are, 

R = 100m 3 + 10m 2 + m^ 

C = 100n 3 + 10n 2 + 


The mode duty cycle for integrators in IO is 



The insertion of the hold interval avoids timing prob¬ 
lems between mode switches for normal and com¬ 
plementary integrators. It also has another important 
purpose as shown below. Suppose a complementary 
integrator is tracking a problem variable, x, 



times the period of the divider output. 

Note: The connection to patchboard C terminal turns 
on mode switch compute light during C portion of cycle. 

ITERATIVE OPERATION 


and x is rapidly changing. Then if the C-integrator 
goes to compute when the N-integrator goes to initial 
condition, the C- integrator will not store the final 
value of the N-integrator due to the time constant of 
the IC-circuit: 


In iterative operation (IO), results obtained during or 
at the end of one solution of the problem are used to 
change parameters or the circuit configuration (switch¬ 
ing) for the next solution. The term IO has an obvious 
derivation from mathematics. In IO, integrators and 
DAC’s are used for memory. Thus, in particular, there 
is a need for integrators with opposite or complemen¬ 
tary logic. These are called complementary integra¬ 
tors. They must be in the R(IC)-mode when the pro¬ 
blem is being solved and in the Hold-mode while the 
normal integrator is in the R-mode. The C-mode can 
be used instead of the H-mode if the integrator has no 
compute input, but only an initial condition input. 



The chart below shows the relationship between inte¬ 
grator modes 

Normal Complementary 

Integrator Integrator 


IC 

H 

C 


C 

H 

IC 



N-OUTPUT 

C-OUTPUT 

N-MODE 


This problem can be solved if the C-mode command to 
the C-integrator is delayed: 



N-MODE 


Without the insertion of the hold interval this cannot 
be accomplished. In some problems, both the N-, and 
C-integrators are used for computation, in sort of a 
bang-bang fashion. For example, if there are two sub¬ 
simulations, each of which depends on the final values 
of the other, then both types of integrators are used 
for equation solution, not merely memory. In this 
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case delayed commands for both types are necessary. 
If all the variables are sufficiently slow, then the need 
for delayed commands is obviated. The implementation 
of mode duty cycles and delayed commands will be dis¬ 
cussed presently. 

The programmer needs a switch at the console to start 
10. This is the RO mode switch when the RO selector 
switch is in external. It follows that there must be a 
"safe" mode when not in IO. This is the standby or 
S-mode and is activated by the reset switch for itera¬ 
tive integrators. Consider the circuit 



It is plain that the N-integrator must be in IC in the 
S-mode. The C-integrator cannot also be in IC in the 
S-mode because the two integrators would form a posi¬ 
tive feedback system resulting in instability. The C- 
integrator cannot be in compute since it may have an 
initial compute input and consequently would integrate 
to overload. Thus, the C-integrator must be in hold 
during the S-mode. Finally, there may be non-IO 
integrators in the program and these must remain in 


Compute or Hold during the iterative cycle. The fore¬ 
going can be summarized by 



The patching diagram for 10 integrator control is 


F 



RESET 


HOLD 


NON-IO WITH HOLD 


F 

F 

/j RT 

^IRT 

R R<Kjl RO F 

RO F 

o~ — -Nb-—o 

RESET 

RESET 

R F 

R _F 

o-o—-o 

o-<£-— 

HOLD 

HOLD 

NON-IO WITHOUT 

NORMAL 10 

HOLD 




R I F 
^ RESET 



COMPLEMENTARY 10 


The F-RT terminal is energized when the machine 
mode is Reset, putting N-integrators in IC and C-inte- 
grators inHold. The F-Hold terminal will be energized 
when the Hold button is pushed, putting both N-integra¬ 
tors and C-integrators in Hold. There is no conflict 
with the RO C-, R-terminals since these terminals are 
disconnected from the patchboard in any mode except 
RO. When the machine mode is RO, the problem Hold 
or H-terminal on the patchboard will be connected to 


the F-Hold terminal. Consequently, to put both the 
N- and C-integrators in Hold, it is necessary that 
the H-terminal be energized and that neither of the 
RO C-, R-terminals be energized by C, R inputs 
from the patchboard respectively. The diagram above 
shows patching for various integrators with mode 
relays. If electronic mode control is used, then the 
relays are replaced with electronic switches as shown 
in Chapter 6. 
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Fig. 1. Iterative Operation Logic Control 


Assume FS1 (Function Switch No. 1) is up and FS2 is 
down. When the machine is in the Reset mode, the out¬ 
put of the F-RT T terminal will be a logical 1. This out¬ 
put is used to reset the divider and counters and to 
place both FF1 and FF2 in the R-state. The clock is 
off. FF1 and FF2 form atwo-bit binary counter. (The 
IO cycle H, C, H, IC corresponds to the counter out¬ 
put 00, 01, 10,11.) Gate H is ready. When the RO 
mode switch is activated the decade counters will start 
and both N, C-integrators will be in Hold. At 10H 2 + 


the coincidence at Gate H will reset the decade counters 
and advance the binary counter to 01. (The use of the 
output of Gate DH will be discussed later.) The output 
of the binary counter will ready Gate C, thereby putting 
the N-integrators in Compute, and the C-integrators 
in IC by energizing the C-terminal through its driver. 
At 10C 2 + C^ the coincidence at Gate C will reset the 

decade counter and advance the binary counter to 10. 
The output of the binary counter will ready gate 1 and 
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put all integrators in Hold. At 10H 2 + H^ the coinci¬ 
dence at Gate H will reset the decade counters and 
advance the binary counter to 11. The output of the 
binary counter will ready Gate R, put the C-integrators 
in Compute, and put the N-integrators in Reset by 
energizing the R-terminal through its driver. At 
10R 2 + R^ the coincidence at Gate R will reset the 

decade counters and advance the binary counter to 00. 
The circuit is now ready to go through the same cycle 
again. 


If FS2 is up, then the IO cycle will start with C (i. e., it 
will skip the first and only the first H interval). 


Delayed Hold commands are provided for the delayed 
complementary integrators by gate DH and the associ¬ 
ated FF. The logical equation for gate DH is, 

DH = H • DH X • DH 2 . 


Note that 


10DH 2 + DH X < 10H 2 + H x 

so that the FF will be toggled before the end of each 
Hold interval. This causes the Reset period of the de¬ 
layed complementary integrators to bracket the Com¬ 
pute interval of the normal integrators. 

Delayed Hold can be obtained by other means. Delayed 
Hold can be generated for both N-, and C-integrators 
by delaying the deactivation of the C, R inputs, res¬ 
pectively. The circuits for obtaining these delays are 
shown in the following diagrams. The inverters are 
those associated with the preceding gates. 


A 

o—#■ 





O.S. 




B 

o- 


-W- 



O.S. 




If FS1 is down, the logic switch found on the counter 
module can be used to single-step the computer through 
the iterative program. 

If the machine Hold switch is depressed, all integra¬ 
tors, IO and non-IO, will go into Hold. This happens 
because the R, C inputs at the patchboard are discon¬ 
nected from the F-RO and R-RO terminals in the inte¬ 
grator modules and the main machine mode takes pre¬ 
cedence. The clock input to the divider is also stopped 
because the RH terminal output is a logical 0. This 
saves the state of the divider, counters, and two-bit 
binary counter. Consequently, the iterative program 
can be put in Hold at any time while preserving its 
status. To return to the iterative program after being 
in Hold, it is only necessary to activate the RO mode 
switch. 

Sometimes it is desirable to control either the C or the 
IC-interval or both by an event in the program. This 
is accomplished by using comparator outputs instead 
of decade counter outputs as inputs to the appropriate 
gates. 


A simpler but less flexible IO logical control program 
is 


COUNTER 





frt' >-4 




H 

C 

R 

ov 
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This program will permit only an IO cycle composed 
of 10 time intervals. As before, the RH f output is 
used to gate the clock. 

MODE CONTROL NOTATION 

When delayed commands are to be used for an integra¬ 
tor it will be denoted by: 



depending on its type 

Integrators can be controlled separately by pairs, not 
under general IO control. 

In this case the following notation is used: 


ACCUMULATOR 



K 

This program will compute the sum of n successive 
values of f: 


The pot is for scaling. The output, e Q , as a function 
of time is 


Q> 

R = LOGICAL EXPRESSION (1C RELAY) 

H= LOGICAL EXPRESSION (HOLD RELAY) 


For example, R = L^ • L 2 * The drivers are still used 

as logical and power buffers. If a resistor rather than 
a capacitor is used for the feedback element of an inte¬ 
grator amplifier, then the mode control relays can 
still be used. The notation is 




R = LOGICAL EXPRESSION 
H = LOGICAL EXPRESSION 


Thus, a summer can be made to act as a sampler in a 
sampled data system by controlling the hold relay. 

MEMORY CHAIN 

A memory chain program is 



This circuit will store a sequence of values of f (t). The 
first integrator requires a delayed command to obtain 
the final value of f when the program which generates 
f goes into Hold. The chain can be any length desired. 



Note that the complementary integrator has two IC in¬ 
puts. The second is obtained by patching a resistor 
between the output of the normal integrator and the U 
terminal of the C-integrator. 

AUTOMATIC RESCALING 

Sometimes it is necessary to rescale a problem during 
the computer solution in order to meet the accuracy 
requirement imposed on the simulation. This occurs 
because of the wide range of variables involved. One 
good example of this situation can be found in the sim¬ 
ulation of a nuclear reactor start-up, where the neu¬ 
tron level may rise through six decades of power. A 
technique is presented here for automatic rescaling. 
Although only a simple example is used for illustra¬ 
tion, the technique is applicable to the reactor start¬ 
up simulation and many other similar problems. 

Suppose it is necessary to generate a function, f(t), 
where 


f(t) = K n e at 

a > o 

K = amplitude scale factor, 
n 

and as the function exponentially increases and reaches 
some limit, f (t) , it is required to perform ampli- 
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tilde rescaling and proceed with the simulation. The 
following program automatically fulfills these require¬ 
ments . 



COMPUTE INTERVAL 
GATE OF LOGICAL 
CONTROL PROGRAM 


The integrator required for the generation of f(t) is a 
normal integrator. The compute interval is controlled 
by a comparator. The function switch is in the Up 
position only at the start of the simulation and the ini¬ 
tial IC of the integrator is entered from potentiometer 
1. The switch could be replaced by a function relay 
energized by the FR bus. The IO logic program (Fig. 1 
or equivalent) goes through its cycle and waits in the 
Compute mode until an external signal commands it to 
go into the Hold mode. Note that the comparator out¬ 
put is used instead of the counter output to activate 
Gate C (Fig. 1). This signal occurs when the output 
of the integrator reaches a limit set on potentiometer 
2 and the comparator is triggered. A complementary 
integrator with a delayed command learns this value 
and feeds it to the U terminal of the normal integrator 
through potentiometer 3 and a resistor R. 

As the logic goes to the IC mode, a new initial condi¬ 
tion is established and the entire process is repeated. 


We see that time can be used as a dummy integration 
variable. Suppose g(x) is the solution of 

M(d) g(x) = b(t). 


That is 


g(x) = 1 M -1 (D)b(t)dt. 

Then f ) g(x)} can also be generated, using time as a 
dummy variable, from the simultaneous solution of 


.g(x) , 

f(g(x)) = j L 1 (D)c (t)dt 
f x -1 

g(x) =J M (D)b(t)dt. 


The amplitude scale factor K n is given by the follow- The P r °g ram for the generation of f(g(x)) is, 
ing equation n 


K 


n 


r (0.056) n 
K l 


n — 0, 1 , 2, 3- 


where R is in megohms, <5 is the setting of potentio¬ 
meter 3 and n is the number of scale changes. 

FUNCTIONS OF A DEPENDENT VARIABLE 


This technique is not highly accurate, but is useful if 
precision is not a requirement. Suppose we want to 
generate 


f(x), x = x(t) 

If f is the solution of an ordinary differential equation 
L(D) f(t) = c(t) 

Where L(D) is a differential operator, then 


f(x) 


r A -1 

J L 1 (D)c(t)dt. 



The comparator is energized when g (t) = x, so that 
t = g(x). Therefore, the output of the f circuit at this 
time will be 


f(g(x 2 ))- 
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Example: f(x) = Ke , -b 4 x ^ b 



The comparator will activate when t = x. The output 
of integrator N1 is Ke _t so that the output of integra¬ 
tor N2 is Ke' X . If these IO integrators compute at 
high speed, say with 0.001 fxfd capacitors and elect¬ 
ronic mode control, the output of N2 is a closed stair- 

i 

3 

Example: f(x) = x , x ^ 0 



"Xlt # 

case approximation to Ke for slowly varying x(t). 
This IO program can be used as a subroutine for a 
non-IO main program. 

The same general comments apply to the following 
examples. 









PERIODIC FUNCTION GENERATION 


Substituting (2) in (1) 


Suppose a non-10 program requires aperiodic func¬ 
tion. We can use an 10 integrator and a function gen¬ 
erator to provide the periodic function. The function 
sub-program is 


t-RATE 



Here the N-integrator is under regular IO control. The 
main program will see the output of the DFG as a per¬ 
iodic function because its non-IO integrators are in 
Hold during the time the N-integrator is inIC and Hold. 
In fact, the Hold can be eliminated if there is no other 
10 calculation to be performed. The non-IO integra¬ 
tor must have the logical control 


0> 

R=F-RT 
H 2 F-H 


so that it will be in Hold, not IC, when the IO integrator 
is in reset. 

NEWTON'S INTERPOLATION FORMULA 


f i = f < x i> 

f N= f o 


A f l‘ f O 


f (x) = f [x 0 ] + (x-x 0 )f [ V xJ 

+ (x-x o ) (x-x 1 ) f jx o , x v xj 

f(x) = f [x o ] + (x-x 0 )f [x o , xj + — 

+ (x-x 0 )'*'(x-x n l )f [x o ”«,x n ] + E(x) 


where 


[ -1 X W 

V*l] ■ I-TT 


E(x) = (x-x 0 )° * ° (x-x n )f x 0 ***,x n) x 


If f has polynomial form, then for sufficiently large n, 
E=0. A good approximation for f, in any event, is ob¬ 
tained by assuming E=0. The result is called Newton's 
interpolation formula. This formula can be used to 
provide a continuous representation of a function de¬ 
fined at n points. The intervals x., x.-l need not be 

equal. However, to simplify our derivation we wiU 
assume that these intervals are equal. It should be 
remembered that the continuous memory to be devel¬ 
oped can be used for unequal sample intervals with dif¬ 
ferent values for coefficients in the computer pro¬ 
gram. We let x =o and assume all x.-x. ^ = a. 


(E-l) n f 

f [V*" X n]=^^’ E V f n 


[V X 1> X 2] 


f r X l’ X 2] - f [V X l] 


x 2 ” X q 


and generally 


1 tVH .*n] ' 

'[w"' 1 .] - 1 [v•'vi] 


X - X 

n o 


We have 


f(x)-Ip o ] * (x-x o )f [x o ,x] -(1) 

'[Vl-'M ♦ <*- x i> , [v x i, x ]— <2) 




“ (E-l) m f 


. (x-x.) 

1=0 l' 


so that 


f §= f o + E 


n nul 

-— L-z n (2L.i) (3) 

^ m! i=o 'a 


For a specific number, k, of E a f Q a representation of 
f(^) can be obtained for (3). For n=3, we have 


and generally 


f [V • ‘ ,x n-l ,x ] f [V- • ’’ X n] + 

[V * *’ x n’ x ] 


(E-l) 2 f 

f $ = f 0 + ^-^o^ + — r -2 ^ f _ x) + 


(E-l) 3 f 
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collecting terms 


L< E - i)f o 


(E-l) 2 t 0 2(E-» 3 J 0 ] a 

2 + 6 J 'a' 


(E-l) 2 f o 3(E-1) 3 : 


— (~) ! 

J V 


(E - 1)Jf o /X«3 
6 V 


f « = f 0 + (i f 3-| f 2 + 3 f l-r f o )x 


+ (-| f 3 + 2f 2-| f l + f o )x2 


+ ^6 f 3 " 2 f 2 + 2 f l " 6 V x 



In the analog computer we can retain a continuous rep¬ 
resentation, f, of f by sampling J at n points (not 
necessarily uniformly spaced), f can be generated 
later on command with the interpolation formula. The 
computer program which implements this operation 
from four uniformly sampled values of f is shown 
below. 



Comparators can be used in conjunction with one ramp 
integrator to de-energize the R-relays of amplifiers 1 , 
2, 3, 4 in the proper time sequence. The integrators 
are switched from the initial condition mode to the 
compute mode when it is desired to M playback fT f(t). 
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If, instead of T(t), it is desired to playback f (T-t), the 
summers and switches are added to the integrator 
chain as shown below. 



The procedure is the same as before except that the 
switches are energized at t=T. 

Although only the circuit for four uniform samples of 
f is shown, the technique can be extended to the case 
of n non-uniform samples. The pot coefficients can be 
calculated previously with a digital computer or desk 
calculator. 


15-11/12 



CHAPTER 16 

SAMPLED-DATA SYSTEM SIMULATION 


Sampled-data systems result from the utilization of 
telemetry, digital control, multiplexed controllers 
and other devices of a similar nature in a control 
system. The analog computer is a useful tool for the 
study, evaluation and synthesis of both linear and non¬ 
linear sampled-data systems. 

THE SAMPLER 

The single-rate sampler is a switch whose periodic 
duty cycle is shown below. 



It f(t) is sampled by this switch, the output function 
is denoted f*(t). 



In most sampled-data systems, At is small com¬ 
pared to T. In classical analysis At is assumed to 
approach zero with finite area, so that f*(t) = f(o), 
f(T), ---. The introduction of this assumption may 
or may not significantly alter system behavior. When 
the analog computer is used for analysis, the sampler 
can be simulated exactly. The H-relays for the inte¬ 
grator amplifiers can be used to represent samplers . 
The duty cycle of these relays can be automatically 
implemented with the digital logic. 

The sampling operation produces sidebands. 


Lf* = F < s + i nw s ) > 

n= -°° 

(L indicates Laplace transform) 


2 t r 

where cj = In order that the sidebands do not 

s T 

overlap, we must have 

Ff - 0, H > fr . 

(The operator F is the Fourier transform.) 

If there is sideband overlap it is not possible (with 
techniques considered here) to recover f(t) from f*(t). 
Of course, in practice, it is usually not necessary to 
recover f(t) exactly. Consequently, judgment is re¬ 
quired to determine the extent of allowable sideband 
overlap. This is one question which is conveniently 
answered with the analog computer. 

FUNCTION RECOVERY 

An ideal filter for function recovery has the Fourier 
transform 

T(«) = A(w)e^^ 

A(w) = 1, M«a> c 
= 0, M>co c 

0(ct>) ~ -k Cl). 

Thus, if F(a>) is the Fourier spectrum of f(t), and 
F(a>) =0, |w| ^ a> 

then the Fourier transform of the recovered function, 
F'(«), is 

F'(w) = T(“)F*H = F(o>) e~ jk “ 

so that 

f(t) = i f (t - k). 

T(o>) represents an ideal low-pass filter and a "cardi- 
nal data hold" if k = 0. 

The impuse response of the ideal low-pass filter is 
« c Sin K <k - t)] 
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The filter is not physically realizable due to the pre¬ 
cursor (response for t < 0). Ideally < w c < ( w s “ ^ 

where co ^ 2 a;.. w should be sufficiently greater than 
s Is 

so that an adequate practical filter can be found. 


The foregoing applies to piecewise continuous func¬ 
tions. Special filter techniques such as clamping 
can be used for discontinuous functions Hold cir¬ 
cuits are useful approximations of the ideal low-pass 


filter. 


T(jw) = Te 


-jot>T 
2 


Sin 


CJrp 

T 


G(^)e 




GH 


Sin (^y~) 




Zero-order hold filter: 


</>H = 


coT 

2 


Let the zero-order filter satisfy the relationship 


f(nT+t) = f(nT), 0£t < T. 


That is, the hold circuit senses the value of f at 
t = nT and holds this value until a new value is 
sensed at t = nT + T. Suppose f(t) = u^ft) so that 

F(s) = j - 


Then 



A computer circuit which simulates a practical zero- 
order hold filter is shown below. Integrator modes 
are controlled by the IO logic. 



If the transfer function for the zero-order hold cir¬ 
cuit is T(s), then 


T(s)F*(s) = ~ 


- T(s) 


1 

sF*(s) 


1 



Thus 


T(j co) 


l-e- j " T 

j“> 


T(jw) 






1 M. Schwartz, "’Information Transmission,Modulation, 
and Noise, ” McGraw-Hill, 1959. 


If an ideal hold filter is to be simulated, then the main 
computation must be interrupted while the memory 
loads. 

First-order hold filter: 

For the first-order hold filter, we have the predictor 


f(nT +t) = f(nT) + £ [f(nT) - f(nT - T)] , (K t < T. 

If f(nT) = 1 for all n, then the output o(t) from 
the filter is 


o(t) = 1 + y > 0«t<T 

* 1 , t ^ T. 

Thus, if the transfer function for the filter is T(s), 
then 


O(s) = T(s) F*(s) 
1 


F*(s) - 


1 - e 


-sT 
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0(8) = (± + A ) (1 - e' sT ) 
S Ts Z 


T(s) = (| + A ) (1 - e‘ sT ) 
s Ts 2 


It can be shown that if T(ju>) is expressed 


T(joj) = G(w)e 


j <*>(«) 


introduction of fictitious samplers or the use of the 
modified Z-transform. The analog computer is an 
effective tool for evaluation of ripple since the con¬ 
tinuous time response of all elements of the system 
is available for observation. 

PULSE TRANSFER FUNCTION 

In what follows, it will be assumed that all samplers 
are synchronized and have the same period T. In 
order to find the pulse transfer function of an element 
of a sampled-data system, consider the simple sys¬ 
tem below 


then 


GO 


w) = \L 2 +i 


Sin ^ 


1 2 


o/T 

2 


<p(u) = - cj T + tan " 1 co 



<x) 


K-th order hold filter: 

The predictor for the k-th order hold filter is 
k 

f (nT + t) = ^ —7 A m f(nT) , o ^ t < 
m = 0 

where A is an operator such that 


nT 


r (t )- 


We have 


r*(t) 


G (s) 


C (t) 


■C*0) 


C(s) = G(s) R*(s) 

oo 

C*(s) = ^ C(s + jnw g ) 


n= - oo 


C*(s) = fjr G(s + jno a ) R* (s + jnoj s ) 


n= -oo 


Note that R*(s + jno> ) is periodic with period w . 

s s 

(Consider the sideband structure.) 

Thus 

R*(s + jnco g ) = R*(s), for all n 

i J 80 , 

C*(s) = R*(s) y ^ G(s + jnw g ) ... Equ. 1 


n — — oo 


C*(s) - G*(s) R*(s) 


Af(nT) = ^ f(nT - T). 


Practical low-pass filters: 

Recovery can also be accomplished with a variety of 
low-pass filters. The analog simulation of the sampled- 
data system usually will provide a means of deciding 
the filtering requirements. Various types of filters 
may be tried in the simulation. 

RIPPLE 

Various responses in the sampled-data system may 
contain components that are either completely or 
partially hidden by the sampling process. These are 
called ripple and may have a bearing on system 
stability. In classical analysis, the system response 
between sampling instants is examined either by the 


From the relationship between the Laplace and Z- 
transforms, it follows 


q( z \ - G (z) 
G(Z) ' R(z) 


Where G(z) is the pulse transfer function and is the 
Z-transform of g(t), the impulse response. For con¬ 
venience, we will use F(s), F(z) to denote the Laplace 
and Z-transforms of f(t), respectively, although F(s) 
and F(z) have different functional form in general. 

PULSE TRANSFER FUNCTION OF VARIOUS LOOPS 

For convenience, we define 


z{F(s)} = Z {L _1 F(s)) . 
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The system 


(see equation (1) above). Thus 



satisfies 

C(s) = G 1 (s)G 2 (s)R*(s) 

C*(s) = [g i (s)G 2 (s)R*(s)J * . 
C*(s) = [g^sJG^s)] *R*(s) 


= Z {G^sjG^s)} . 

It is convenient to denote Z {G^sjGgfe)} by G^Ggtz). 
It is important to note that G^GgCz) * GjCzOGgte). 

The computer block programs for the simulation of 
various simple sampled-data systems together with 
the closed-loop pulse transfer function are shown in 
the examples below. In each case the sampling func¬ 
tion is implemented with mode relays under IO logic 

1 TT 1 

control. In all cases G = - + a ’ n ~ ~s ' 

(The symbol / for an integrator means it is not under 
IO control. The normal IO logical program is used 
for the sampling amplifiers). 


Example: 


r (t) - 


C(t) 


—►(: (t) 
C(z) = GR(z) 


r (♦)—^ ^► 


C*(t) 


Example: 
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Example: 


Example: 











DESIGN TECHNIQUES 


Hoot-locus: 

Root-locus plotting techniques apply without modifi¬ 
cation in the z-plane. The region of stability is 
|z| <1. The interpretation of the locus is somewhat 
different. The family of circles 

(where G (s) is the compensating function) is given 

IZ | = cr = a T _ v C _, 

‘ ' by T(z). Then since 



represent constant damping, a, (not damping ratio) m/ C(z) 

for the response T ^ z ' = R( z ) = 1 + G c G(z) 


f* 



* 


we have 


The family of radii 


G c G(z) 


T(z) 

1 - T(z) * 


Arg z = 0 = oo T 

represent constant frequency for the above response. 
Thus, £,, the damping ratio is given by 



or equivalently 



This technique is useful when G c is realizable. 
Digital Compensation: 

A digital controller can be used instead of a com¬ 
pensation network above. This controller can be 
simulated with the analog computer if it is suffi¬ 
ciently simple. If not, then a hybrid system may be 
used with the analog computer simulating the remain¬ 
ing sampled-data system. 

MULTIRATE-ASYNCHRONOUS SYSTEMS 

Multirate and asynchronous sampled-data systems are 
at best difficult to analyze with classical techniques. 
The analog computer, with its asynchronous control 
structure finds powerful application for such problems. 


Pole-zero cancellation techniques: 


Suppose the desired pulse transfer function for the 
controller 
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CHAPTER 17 

PARTIAL DIFFERENTIAL EQUATIONS - PART 2 


PROPAGATION OPERATOR 

There is a class of equations which can be solved with 
the similarity transformation 


<*>(R) = (T _1 A X T) 0(9R) 


where T is a suitably defined integral transform over 
3R (boundary of the region, R) and propagation into R 

is alongthe x 1 coordinate of R. That is, the surfaces 

x = constant are’’concentric” to aR. (T~ A T) is called 
the propagation operator and A the propagation function. 

”1 x 

That T AT exists, can be shown if R is homogeneous 
in the boundary coordinates. That is, usually R can 

be non-homogeneous only in x 1 . We will use this ap¬ 
proach to find Green’s functions for the solution of dif¬ 
ferential equations where we know the required inte¬ 
gral transforms over aR. 

First we derive A for a few simple cases for the pur¬ 
pose of illustration. 

Consider the cases: 


u = 

-u 

(1) 

yy 

XX 

u , = 

u 

(2) 

yy 

X 

u 

u 

(3) 

yy 

XX 


We will use the Fourier transform and its inverse. We 
claim (for propagation along the x coordinate) 


u(x, y) = T 1 a x GH 


Then, 


yy 


XX 


= 4 ] A X G(co)e j " y du,. 


J -w 2 a x G(co)e jay do> 

- oO 

J in2A A X G(w)e^ doj 


j In A A X G(o;)e^ do; 


For (1) above 

1 f / 2 i 2, 

U yy +U xx = 2* [ + 111 A) ‘ 

A X e*^ G(w) do? = 0 

for arbitrary G. A sufficient condition for the vanish¬ 
ing of the integral is 

1 2 a 2 

In A = cj 


and 


A = e 

u(x, y) = t' 1 e" X G(co) 
u(0, y) = T' 1 G(o>) 


so that 


0(0)) = (T)u (0, y) . 


Consequently 


u(x, y) = (T _1 e‘ l“l x T) u(0, y) 

and(T"^e"l w l x T)is the desired propagation operator. 
For (2) we have 


jcjy 

j (-<*? - lnA)e G(^) da) = 0 


so that for arbitrary G 


A = e 


and the operator is 


For (3) we have 


T~ 1 e~ w X T 


2 1 2 A 

“OJ = In A 


In A = ± j o) 


i w 

\ = e ‘> *2 


-jw 
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Let 


so that 


u(x, y) = T' 1 ^ G x H + T* 1 G 2 (co) 


u + k u =0 

yyyy xx 


= T^e^M+T " 1 e- jC ° x G 2 (a,) 


Then 


u x (x, y) = T G^w) - T *j«e ^ uX G 2 («) 

u (0, y) = T" 1 G 1 (o.) + T" 1 G 2 («) 


u (0, y) = T j w G^(o)) - T jo> GgM 


G 1 H + G 2 (o>) = (T) u (0, y) 


G^co) - G 2 H = ~(T) u x (0,y) 


G l(“) = J [(T) u (0, y) + p (T) u x (0, y)] 

g 2 m = i [( T ) u (°>y) - ^(T)u x (o, y )] 

u(x, y) = ^T _1 (T)^ u( 0, y) 

+ ( T_ 1 lj^ (T) ) u x (0 ’ y) 


(t _ 1 u(0,y) 


- ( T_1 ^ (T) ) u x (0 ’ y) 

, y) = ^ T _1 coscox (T)^ u(0, y) 

♦ (t- 1 ^(T))u x ( 0,y). 


Let's examine one more case 
2 


£_ (EI } + m ®-« = o . 

9y ay ax 


Let El, m be constants and define 

m 


k = 


then 


4 2 

co + k In A ~ 0 


2 2 

o 

VT VT 


9 4 U) -CO 

ln 2 A = J —. — 




Vk” Vk” 


and we see from the last case 

2 


u(x, y) = ^T -1 cos(^^ )(T)^ u(0, y) 


2 

• / w x \ 

, sm( vr) , \ 

+1 T \ -Y-] (T) j u x (0, y) 

VF 


I 


SUMMARY FOR POLAR COORD'S (r, e). 


2ir 


(T) f (e) = / f(t) e' iH dt, = F(n) 


T' 1 F(n; = F(n)e jn0 = f(e) 


n=-° 


For v<f> = Owe have 

u(r, e) = (t _ 1 r‘l n l(T)) u(0, e) 
2 2 

For V (p + k <p = Owe have 


u(r, e) = (' 


1 J (kr) \ 

T [ j R (k) ] (T) / u(0 ’ e) 


El 


HIGHER DIMENSIONS 


We can derive similar results in higher dimension with 
the use of suitably defined transforms on aR. 
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GREEN’S FUNCTIONS 


We proceed to use the propagation operator to obtain 
Green’s functions. 


Case I: u + u =0. 

- xx yy 

u(x, y) = (t‘ 1 e" I" I X T ) u(0, y) 

oo 

J u(0, ;) e - ^ 4 dr, 


fi 7 u(0, 4) dt) 


j e‘H x + 


so that 


/ e ia,( y- 4) e” wx dc 


/ +cux io/y-r,) , 
e w dcj 


1 1 _ 2x 

x+j(£-y) + x-j(c-y) " x 2 + ( y . t ) 2 


_ l ( (x)u(o, U 


u(x, y) = - / 


7r ^ tit \£ 

X +(y-U 


Case II: u = u 
- yy x 


U(x,y) = 2 ^ / e W X e Ja,y daj 


/u(0, Ue'^dt, 


/e~“ x e^ y ’ 4 ^ = / e" w x cosw(y-U dw 


oo 2 

+ j / e w x sinco(y-£) da; 


OO _ 2 

2 /e" w X cosa;(y-^) da; 
0 


-(y-cr 

1 4x 

-e 


2 V?rx” 


so that 


oo ±<L 

»<x,y) - — / e ta iM-dt. 

4irV7T - 00 Vx 


Case ID: u = u 
- yy xx 


u(x, y) = * coseux u(0, y) 

* (t- 1 ^ t) „„<(>,,) 

Let f(x, y) = 1 cos^x u(0, y) 

g(x, y) - (t- 1 t) »„(0,y) 

a 

f = 4- / cosojx e^ wy do) 


/ u(0, 4)e" iw4 dt, 


■4 


/ 


jo;(y-£,) , 

cos<^x e J J do> 


= 4 jfu(0,l)dt 


°° . / v 00 

J*cosa;x e Ja * y ’^da; = / eosa;X COScj(y-C) da> 


/ e -“ x e ja,(y - 4) da, 


J’coswx e^ y_4 W> = n (x-y+£) + s (-x-yH)j 
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and 


Finally, 


o° 

f(x, y) = Y / u (°> £) [& (x-y+c) 

— 00 

+ 5(-x-y+t,)] d£, 

f(x, y) = j [u(0, y-x) + u(0, y+x)J 


Also 


, x 1 c sinojx ^ jwy 

g( x >y)=i; J—^~ e dco 


27T 


( 


Ju x (0, 1 ) e' jwt d4 


= i / u x < 0 ’^ [cosc^y-;) 


sin w(y-t,)] 


+ j sin co(v-U I dco 

for y-£>0: 

0 

' sincox C0Sc ^y_^ do) = 7T , y - £, < X 


f sincj: 

J w 


= 2ir , y - t, = x 
= 0, y - £, > x 


for y-£,<0: 


fsirKJX cos ( y .^) d co= n , £ - y < x 

J Cx) 

- 00 

■= 2ir , £, - y = X 
= 0 , 4 - y > X 

for the first case (y-f,>0) 

$- »,y>4>y-x 
for the second case (y-£<0) 

5 = » » y + x > t. > y 

so in order to meet both conditions 

J - »’y + x >^>y-x 


and 


1 J +x 

g(x, y) = g J u x (°> ^ 
y-x 


u(x, y) = i- ju(0, y-x) + u(0, y+x)J 

1 y+x 

+1 ( u x <°’ ^ • 

2 J y-x x 


Case IV: 


Consider u = u in the region 0^y^7r, O^x. 

yy x 


Tf(y) = jf(4) sin n£,d^ = F(n) 
0 

T -1 F(n) =- S, F(n) sin ny 
7r n=l 


u(x, y) = T _1 A X G(n) = j J^A^sinny) G(n) 

u =- 2 A X (-n 2 ) (sin ny) G(n) 
yy 7 r n= 1 


u = — 2 A x log A (sin ny) G(n) 
x 7 r n=l 


We have 


log A = -n2 


A = e 


1 2 

u(x, y) = T’ 1 e‘ n X G(n) 


u(0, y) = T -1 G(n) 
G(n) = (T) u (0, y) 


, 2 

u(x, y) = (T -1 e -n x T) u(0, y) . 


9 00 2 

u(x, y) = - e" n x sin ny 
n=l 


. / l u(0, 0 (sin ni) d£ 


a 

7T 

= - £u(0, 

( 


‘) 


[ £ e' (sin ny) sin 
n*l 


in nt,^ 


17-4 



or to write out a few terms 


EXAMPLES 


u(x, y) = — e“ x siny f u(0, £) (sin£)d£ + e~ 4x sin 2y 

n J o 


.1 f U(0, £,) (sin 20d£ 


u + u =0: 
xx yy 


Let t = £ - y 


U(x , y) -I 


— x.(y-s) 


+ e ^ X sin 3yf u(0, £) (sin 3c)d£- 
0 


All of these techniques can be implemented using re¬ 
petitive operation. The idea is to evaluate at high speed 
the equations for fixed values for the coordinate para¬ 
meters. Then by changing these parameters slowly 
the solutions can be found over the region, R, of def¬ 
inition for the equation. 


u(x, y) 


= 1 J xu(0, y+t) 

71 — °° X +t Z 


u(x y) fxu&y+t) 

( ,Y) ' J 0 x 2 + t 2 

. f *«(», y-t) 


The program is: 


u(o,y tt) 

-TdfgN 


xu (o,y±t) 


-im 

1 





r 

-CM 



H 

x • \ p) 

n 

1 

x 1 
100 

( 

-- 1 

H 



NJi-o u(x,y) 


R H Rn 


o I FUNCTION 

Kjy relay 

FF3 FR 


1 1 I 5 I 6 

COUNTER 
O 


FF MASTER RESET 




The FF T s act as an up-counter with successive trigger The first and second integrals are calculated during 

inputs to FF1. When the machine is in reset the counter the first and second compute intervals, respectively, 

output is zero and the integrators are in the standby The comparator is energized when y ± t reaches its 

mode. (The ROselector switch is in external). When maximum absolute value for each calculation. It then 

the RO pushbutton is activated, the N-integrators go terminates the compute interval. The final output from 

to compute. As the counter is incremented, the mode the accumulator is the function u(x, y) for the par tic - 

sequence is ular x, y chosen for the run. 

COUNTER MODE 

3 2 1 FF 

0 0 0 C 

0 0 1 H 

0 10 IC Let t = i - y. 


u(x, y) = --— J e 4x u(0, 0 dc, 

47rV7TX 


Note that (FF3) = 0 when the counter has reached 101. 
When this occurs there are no more trigger inputs to 
FF1 and the program stops. Before this, the trigger 
for FF1 is provided by the clock (L 2 ) when the machine 
is not in compute and by the comparator (Cj) when the 
machine is in compute. There are only two compute 
intervals. The function relay (FR) is energized (+) 
during the first compute interval and de-energized (-) 
during the second compute interval. 


The analog program is 


-lOO 


GATE I 

-lOO 


Both the logic program and the principle of operation 
of the complete program are the same as for the pre¬ 
vious example. 




Oil H 

10 0 C 

10 1 H 


1 f . (yl£) u ( 0 r) 

u(x ’y )= irJ e 4x d ; 
"°° 2 
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CHAPTER 18 

CORRELATION ANALYSIS 


PROBABILITY THEORY 


STATIONARY RANDOM PROCESSES 


Sample space : 

A sample space, S, is the collection of all outcomes of 
an experiment. 

Probability: 

P(S) is a function defined over S. P maps S into the 
interval 0 ^ P(S) 4 1. Pis not necessarily univalent. 
The integral (or sum) over the range of P is unity. 
In case S is finite, P is said to assign a probability to 
each point of S. In case S is infinite, P is said to 
define a probability density over S. Suppose S a is a 
nested collection of subsets of S such that S^ C S i + i 

for all i. The extremum of S a = S and the infemum of 
S^is empty. The measure of the range of P for the 
monotone sequence S a is the probability distribution 
for S with respect to the collection S a . 


EXPECTATION AND VARIANCE 


If the statistics of a process are time-independent, 
the process is said to be stationary. 

Suppose f(t) is a random function defined for 

-oo < t < °°. 

Let the domain of f be denoted by D. We define an 

ordered countable set {D } so that 
1 a r 


D. 

l 


x. < t < x. - 
l l+l 


X i +1 


x. > A . 
i 


UD. = D 

l 


If f. — f(D.), then the collection f is an ensemble. If 
11 a 

f is stationary and A is sufficiently large, the ensem¬ 
ble will be ergodic. (That is, the statistical properties 
of the members of the ensemble are identical.) Hence¬ 
forth, all ensembles considered will be assumed er¬ 
godic . 


The expectation of S, E(S), is the measure of SP(S). 
The variance, V(S); is the measure of 

S 2 P(S) - E 2 (S) . 

Example: 

Suppose S is given by -» < x < °° and 

1 

P(x) = (27r) 2 exp 

Then 


oo 



• oo 


V(x) = | x 2 P(x)dx-/x 2 = <r 2 

- OO 

Let the S be defined by -°° < x < x , then 

ot a ’ 


FIRST AND SECOND PROBABILITY DENSITIES 

Let S be the sample space consisting of all possible 
values for the ensemble jf.(t. + A)| where A< A . 

P(S) is said to be the first probability density. P(S) 
is independent of A, since f is stationary. The ensem¬ 
ble is ergodic so that 

6 

E (S) = £5. 27 

-6 

5 2 

V(S) = 27 5 f2 ft)* - E < s ) 

-5 

Let S T be the sample space consisting of all possible 
pairs of values (f.(t. + a), L (L + A + T)) for the en¬ 
semble where A + T < X. P T (S T ) is said to be the 
second probability density. P^S^) depends only on 
T, not on A. The ensemble is ergodic so that 


fr.-* ) 2 

7JL 




| P(x)dx 


6 

E T ( V = o L i m ~ 27 J f ft) f ft + T)dt. 

-<5 


is the distribution function. 


(It can be shown that this equation is valid for T > X .) 
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AUTOCORRELATION FUNCTION 

E t (S t ), above, is the autocorrelation function for f 
and is denoted by <^(T). That is, 


%< T > = a- 


u 

27 \ + T)dt 


-a 


The autocorrelation function has the properties: 

1) 0 ff (T) = <*> ff (-T) 

2) (o) is the average power of f. 

’ff 


3) 
Proof 


%(T) « <P H (o) 


|% (T) | h j l f(t) l ' l f(t+T) l dt 


Lim 

‘ a— 00 


h \ i * 

-a 

|l(t*T)| 2 - [|l(t)| - |I(t*T)|] 2 jeft 


<%<») - k 5 D«*)i - 

-6 

|f(t +T)jj 2 dt 

4) f € P^ -—► e P„ (is periodic with period v) 

Proof: f(t + T) = f(t + T +v) 


u 

where <#> = Lim ^ \ A* Sin (» n t + 4 > n ) 

5^00 

• ^Sin (c» n t + T + 0 n ) dtj 


Proof: 


Lim 

6-°o 


+ v)- 


A Sin(o) t + T + <#> )dt = 0, m*n 
mm m 


9) If f, f' e bd, then (T) = - 4> ff (T) . 
Proof: 

6 

<?> fT (T) = Lim ^ f f’(t)f'(t + T)dt 

6-* 00 * 

- o 


= Lim 2a 
6~+ 00 


[f(t)f’(t+T)] - 
-6 


u 

Jf (t)f”(t + 

-a 


T)dt 


d T im J- 

'32 Lim 2a 

dT 6— °o 


| f 1 (t)f 2 (t + T)dt 


-6 


CROSS-CORRELATION FUNCTION 


5) contains no phase information about f. 

6) For every f there exists a unique the con¬ 
verse is not true. 

7) If f & P for any v , then Litri (T) « 0 

v x -*► n 

8) If f can be represented by 

f = £ A n Sin(a 'n t + V ’ 
n = -°° 

then 

oo 

% = S 0 nn 

n = - 00 


Suppose f, g are two stationary random functions, f, 
g can be decomposed (as above) into two ergodic en¬ 
sembles f , g . The union of these ensembles is a 

or a 

composite ensemble of pairs (f, g)^. Let S,p be the 

sample space consisting of all possible pairs of values 
(f (t. + a), g. (t. + A + T)) for the composite ensemble, 

where A + T <A {\ is the ensemble duration). Since 
the ensemble is ergodic 

6 

E (S T ) =Um ^ | f(t)g(t ■+ T)dt 

a— °° -a 


E^, is the cross-correlation function for the pair (f, g) 
and is denoted by V T >. 
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" oo 

Lim \ g 2 (t+T) dt 

. <5 ^°° -l 

Example 2 (general random function): 

AUTOCORRELATION FOR A SUM OF FUNCTIONS 


Henceforth we will denote 

Suppose f is defined by 

6 

a) f(n6 +t) = f(nfi), o ^ t < 6 

Lim j f(t)dt 

6~*°° J 

- 0 

b) P [f(n<5) = x] = P(x) 

c) E [f(n6)J = o 

by Mean {f(t)} , or just M{f} wherever convenient. Let 
f(t) =fj(t) +f 2 (t). 

then 

%fT) = Mean {[^(tj+f^t)] [fjft+Tj+f^t + T)]} 

°0 

E [f(4)f(£. + T)] = j* [xP(x)] [x(l- -j dx + 

— oo 

= M (f^tJf^t + T)} + M |f 1 (t)f 2 (t + T)| 

+ M {f 2 (t)f 1 (t + T)| +M {f 2 (t)f 2 (t+T)| 

oo 

jj [xP(x)] [yP(y>] dxdy, |T| < 5 

- OO 

Thus, 

oo 

= [xP(x)] [yP(y)J dxdy, |t| > 6 . 

II l^ l^ Vl X Y2 

- oo 

OO 

c c 

The general rule for a linear combination can be de¬ 
duced easily. 

But \\ xyP(x)P(y) = o, so that 

- OO 
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Example 3 (binary noise): 



S 

- <x 


x 2 P(x)dx = V(x) (1- 



= o, T > <5 


^ff 



Binary noise is a special case of (2), above, where 

p« 4 u i ( x - v + i u i < x+x 0 ) • 

That is f takes on values of x , -x with probabilities 

of ~2 ' 

Then V(x) = X 2 and 
o 

<b t . = X 2 
r ff O 

Example 4 (periodic binary noise): 

Suppose f is defined as in (3), above, except that 
f(C) = f(C + 4 for large v. Then <f>^ is periodic. 



4>ff 



POWER SPECTRAL DENSITY 

The operator F denotes the Fourier transform, so that 
Fg * G(w). 
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# (m) is said to be the power spectral density, 
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WHITE NOISE 


White noise has constant power spectral density of ar¬ 
bitrary magnitude (say A). Thus, the autocorrelation 
function for white noise is 

* g g = Au i (T) ’ 

where u^ is the unit impulse. 
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g(t) e^ Wt dt 


Lim 

<5-*°° 


6 

j g^e - ^ dt, 


PERIODIC WHITE NOISE 


If f satisfies the requirements for white noise for some 
interval 6, where <5 is large, and if f (t+6.) =f(t) r then 


tf> ff (T) = u 1 (T-nd). 

n =-«o 
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FINITE CORRELATION FUNCTION 


If (1) is to be implemented directly, the program 
is 


If f, g are defined only for o ^t < 00 , then 
1 6 

4 > = Lim — (f(t) g(t+T)dt. 


% 


can be approximated by a sum: 


N 

0 . « Mean f(mAt) g(mAt + T) (1) 

^fg m=i 


In this case At must be small enough that the high 
freauency components of f, g are adequately represent¬ 
ed (see Chapter 16). N must be large enough that the 
mean of the sum is a good approximation to 


CORRELATION PROGRAMS 


Simple computer programs for (1) require that 
At > T . (2 ) 


To obtain the best approximation, At = T + e , where 
e is the minimum time required for sampling in the 
computer program. ( 2 ) imposes a limitation on the 
magnitude of T. That is if T is large, then so is At. 
Thus, for large T, relatively few samples are made 
of f, g. This restriction is not unduly severe, because 
0 ^ cannot be determined accurately for large T to¬ 
gether with finite histories for f, g. 



Integrator 1 samples f and then after an interval, T, 
integrator 2 samples g. Integrator 3 is in hold until 
both samples are made and then it samples the product 
f(t.) g(L + T). The output of integrator 3 is averaged 

by aMean Program (Chapter 3) which in turn generates 
0 ^. If T is to be specified accurately, the integrator 

capacitors must be as small as possible to minimize 
the effect of the IC time constant. 

The total sampling time can be reduced by replacing 
integrator 1 with a quantizer whose input is f together 
with dither (see reference), f will be represented by 
a two-bit word (output of the quantizer). 

The table below defines the quantizing function 


OUTPUT 

BINARY 

OUTPUT 

DECIMAL 

OUTPUT 

-2 ^ x <-l 

00 

-2 

-1 ^x < 0 

01 

-1 

0 ^ x < 1 

10 

+1 

CM 

VI 

X 

VI 

t-H 

11 

+2 


18-5 



Except for scaling the program is 



MEAN 

PROGRAM 


RTS 


DITHER 

PROGRAM 


DIVIDER - COUNTER 

|0 |T,|l| 


f-rt': 




ES MEANS 
ELECTRONIC SWITCH 


COUNTER 

TOT 


COUNTER 
l T 3|13 


The Dither Program is found in Chapter 10 and the Mean Program in Chapter 3. 


The correlation parameter, T, and the reset interval 
for integrator 1 are given by (10Tg + lOTg + T^) P 

and (101 3 + lOIg + I^P, respectively, where P is 

the period of the output of the divider. 


Reference: Korn, G. A., "Hybrid-Computer Tech¬ 
niques for Measuring Statistics from 
Quantized Data", Simulation, April 1965. 







CHAPTER 19 

ON-LINE DATA ANALYSIS PROGRAMS 


An analog computer can be used to advantage for on¬ 
line data analysis. Programs are given in this section 
which correspond to the more frequent operations on 
incoming data. The analog computer is used for both 
processing (analysis) and editing. The principles de¬ 
veloped in other chapters can be used for processing. 
The programs here are for simple editing. 


The program which measures error frequency is 


ERROR FREQUENCY 

Suppose an error is defined to occur whenever 
I f(t) - g(t) I > e • 



ERROR PEAK VALUE 

Suppose, with the above definition of error, it is re¬ 
quired to find a sequence of peak values of error, 

Max j f - g | 

Max | f - g | . 

The program is 



+ 100 
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All the programs indicated will be found in various 
chapters of this manual. The integrators in this pro¬ 
gram are under iterative control from the C, R logical 
control inputs. Pot 1 is set to a small value to provide 
positive zero derivative detection. 

TIME BETWEEN EVENTS 

Suppose two events, Ep Eg are defined by 


Ej: f(t) - g(t) = €l 

f(t) - g(t) ^ 0 

E 2 : h(t) - m(t) = 6 2 

h(t) - m(t) s 0 

A program which will determine a sequence of time 
lapses between Ep Eg (not Eg, Ej) is 
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CHAPTER 20 
SYSTEM OPTIMIZATION 


System optimization is a procedure for generating a 
response, y s , of a known experimental system so that 
it matches, as closely as possible, some desired res¬ 
ponse y d . y^ may come from tape or may be the res¬ 
ponse of a known but more complex system. In what 
follows we will assume that either an external device 
generates y d or else we know the computer program 
for it. A measure of the accuracy of fit of y s needs to 
be established. This is called the Criterion Function 
(CF). The CFhas different forms depending on whether 
dynamic or static optimization is required. The exact 
form of the CF in either case is a matter of technical 
judgment and will depend on the particular situation. 
The CF for the dynamic problem will need to take ac¬ 
count of the transient behavior of both y s and y^ over 
some interval of time. Typical choices for the CF are 

STATIC 

F = I K - y/ 


DYNAMIC 

f = So ' y d (t) ] 2dt - 

(Note that the CF is defined so that it is non-negative) 

The optimization problem amounts to finding maxima 
or mininr/a of the C F. Since y s is a function of its para¬ 
meters, say C^, ..., C n , so is the CF. Thus the CF 
can be viewed as a surface over the parameter plane. 
In general the C F will have several maxima and minima. 
We will restrict our attention to finding the nearest min¬ 
imum from some initial point in the parameter plane. 
Once techniques are developed for the minimum prob¬ 
lem, they can be applied with modifications, which will 
be obvious, to the maximum problem. 

We will develop two methods for finding a minimum: 
steepest descent and relaxation. The static and dyna¬ 
mic optimization problems will be treated separately. 

STATIC PROBLEM 

The family of curves f = constant on the surface for f, 
the CF, are constant elevation contour lines. Differ¬ 
entiating 

df = £ -H: dC. = Vf • dC = 0 

where C is the parameter vector — i. e., 

C = [Ci, C 2 , ..C n ] (Remember that the CF is a 
function of the parameters of y s .) 


We see that vf is perpendicular to dC and dC is a dif¬ 
ferential vector which is tangential to the space curve 
f = constant. The projection of the curve on the para¬ 
meter plane is identical geometrically to the space 
curve because it is at constant elevation over the plane. 
The projection ofvf on the parameter plane is perpen¬ 
dicular to the curve. It is evident that Vf corresponds 
to the slope of the surface f and Vf = 0 at the minimum. 
To get to the minimum from some point on the surface, 
then, all we need to do is follow - vf. (If we wanted to 
find the maximum we would follow + Vf.) The pro¬ 
jection of the path on the parameter plane is 



If the system is static (algebraic) we can change the 
parameters of y s while following - Vf without intro¬ 
ducing transients in the system. Now if 



then 

V f = (y g - y d )vy s . 

If we let e = y s - y^, then C can be determined by 
(since we follow - vf) 

c = - e {vy g dt. 

or in scalar form 

C dY s 

C 1 _ ." e fdCj dt> 

. r dy 

C n = ‘ e fer dt ' 
n 
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In practice, it is usually advisable to use gain factors, 
K^, so that 


C 


1 


-K x e 


7c 


dt 


ac P 


ay s 

* c 


d y 
_ 

2 ac„ 


= -y„ 


dt‘ 


d , 9y s (o) « _ 0 9y s (o) _ 0 
dt (_ aoT ) " °> ~dCT ' °' 


• r 

C = -K e \- 77 r- dt . 
n n J dC 

n 


K is chosen so that the procedure converges but at the 
same time is stable. When vf = 0, the integrals will 
be stationary and C will be the parameter vector (point 
on the parameter plane) corresponding to the minimum 
of f. Vy s is generated by differentiating the generating 
equation for y g with respect to the parameters C.. 

DYNAMIC PROBLEM 

Except under special circumstances, vf cannot be gen¬ 
erated on an instantaneous basis. Even if it could, the 
parameters could not be changed during the determina¬ 
tion of f. Since the calculation of f is carried out over 
a finite time interval, varying the parameters would 
lead to a result which is not the CF. 

Some definitions for f will permit Vfto be computed as 
often as is f . For example, if 


then 

vf = \ T 0 (y s - y d> v y s dt - 

As soon as we know the components of y s , we can use 
a steepest-descent method. After each solution a new 
position in the parameter plane is chosen by taking a 
finite step along -Vf. Now Vy s can be generated by the 
parameter- influence equation method. The differential 
equation for each component of Vy g is obtained by dif¬ 
ferentiation of the generating equation for y g . For 
example, if 


These parameter-influence equations can be solved 
with an analog program. Their non-homogeneous terms 
are obtained from outputs of the generating program 
for y s . It will be noticed that they have identical form 
except for the non-homogeneous terms. Due to equip¬ 
ment limitations, it is usually necessary to solve them 
one at a time, switching in the correct non-homogeneous 
inputs from the y^ program in a suitable way. This 
leads to the use 01 the relaxation method. 

With the relaxation method we change only one para¬ 
meter, C i? until the component, (vf)^ corresponding 
to the parameter is zero or nearly so. The pro¬ 
cedure is applied to all the parameters on a cyclic 
basis. An example of the resulting path in a two 
dimensional parameter plane is 


C 2 



meter-influence equations. Nothing needs to be changed 
in the influence equation program except the non-homo- 
geneous term. 


C 1 y s + C 2 y s = °’ y s (o) = °’ y s (o) = A ’ 
ay s ay s 

the equations for ^ can be obtained by dif- 

1 2 

ferentiating first with respect to and then C 

d h " d y c 


d 2 dy s 

- * C ( - ) + c 

aci • c i dt 2 ( a<V c 2 acj 


- y s’ 


d 9y s(0K _ 0 ^ 

dt K dC 1 } ’ dC 1 


= 0 


In the event it is not convenient to use the parameter- 
influence method, successive values of f can be cal¬ 
culated for a sequence of values of Cp This is con¬ 
tinued until two values of f bracket (vf)j = 0. Then 
the procedure is applied varying (L +1 and so on. This 
will also lead to the minimum for f. In either case the 
subscript indexing and switching of non-homogeneous 
terms is done with a digital logic program. In the 
latter case, successive values of f are retained by a 
two or three word memory chain and comparisons made 
by comparators. 

With either method the step size must be determined 
by a gain factor which ensures proper convergence. 
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CHAPTER 21 

MEDICAL APPLICATIONS 

Originally published by: Beckman Instruments, Inc. 
Author: Hiroshi Hara 


ANALOG SIMULATION OF VENTRICULAR PUMP 
ACTION 

Due to the complexity and non-linearity involved, the 
building of an analog model of blood circulation is one 
of the challenging problems in biomedical science. 
This example is by no means intended to cover the 
subject, but illustrates how iterative programming is 
conveniently used in simulating ventricular pump 
action. 


The ventricle has an inlet valve (tricuspid valve) and 
an outlet valve (pulmonary valve). A cardiac cycle 
consists of two periods called the diastole and the sys¬ 
tole, as shown in Figure 1. During the diastole, the 
heart muscle relaxes and blood flows into the ventri¬ 
cle as the inlet valve opens. At the completion of the 
diastole, the systole begins. During this period the 
heart muscle contracts and as the ventricular pressure 
exceeds the arterial pressure, the outlet valve opens, 
expelling blood from the ventricle. An electric analog 
of the ventricle is shown in Figure 2. The analogous 
quantities are indicated below. 

Blood Pressure ~ Voltage 

Blood Flow ~ Current 

Blood Volume ~ Charge 


nth CYCLE , (n+l) th CYCLE 


~ m 

DIASTOLE 

SYSTOLE 

DIASTOLE 

SYSTOLE 

- 1 *- 









FIGURE 1 


The contraction and relaxation of the muscular wall of 
the ventricle is simulated by decrease and increase, 
respectively, of capacitance c^ in some prescribed 

manner. If we let E = ^ = elastance of the ventricle, 

the pump action is assumed to be described by the 
following equations: 


[ E d]n- A - (A - E n-l le7 + B M>, 
l E s]n- E J*V' 


where 

|e = Elastance during nth diastole 
|e J n ” Elastance during nth systole 

E^ = Elastance at the end of n th diastole 
n 

E s ! = Elastance at the end of (n-l)th systole 
n-1 

fadln = Blood volume during nth diastole 


* n = Slope of |E s J n to be determined as a 

function of q^, the blood volume at the 
end of n th diastole. 

A, B, and r are positive constants. 


Since the static work done by the heart at a particular 
value of q« is known experimentally (Starling 1 s Curve), 
must be found as a function of q^ such that experi¬ 
mental results are duplicated. The work done during 
the systole is 


Tsystole 

w s= k / Vl2 dt > 

o 

where V£ = arterial pressure 

i 12 = blood flow out of ventricle 
k = constant 

d 

and consequently for a certain q n , must be such 
that |w g J n = [W s ] Starling. In other words, itera¬ 
tion with respect to is required. 

The passive circuit shown in Figure 2, equations (1) 
and (2), and an iteration circuit are implemented as 
shown in Figure 3. Note that two iterative cycles are 
required to simulate one cardiac cycle. The distinc¬ 
tion between the diastole and systole is made by a func¬ 
tion relay (Relay KO) which in turn is controlled by 
one of the memory integrators in a chain of four. The 
action of this circuitry is shown by the timing diagram 
in Figure 4. Since the compute intervals for the dias¬ 
tole and systole differ, the comparator which controls 
the compute interval of the iterative control program 
must, necessarily, require two different reference 
voltages which are switched alternately by relay KO. 
The function of the M2 - Ml memory chain is to make 
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ELECTRIC ANALOG OF A VENTRICLE 


r-1 

v o 


I_ 

VENTRICLE 

V 0 : VENOUS PRESSURE 
V, i VENTRICULAR PRESSURE 
V 2 ! ARTERIAL PRESSURE 
R 0 , . INLET VALVE RESISTANCE 
R| 2 : OUTLET VALVE RESISTANCE 
C| : VENTRICLE CAPACITANCE = -g-.WHERE E = ELASTANCE 
q .BLOOD VOLUME 

W s ! WORK DONE BY THE VENTRICLE DURING SYSTOLE * 
k / v 2 L \2 dt > WHERE k = CONSTANT 





R OI 

■A/W 


i 01 




1 I2 


AA/V 


C,(t) 
CHARGE, q 


i 12 


V 2 


C 2 


FIGURE 2 


the final value of q (blood volume) at the end of the nth 
diastole available throughout the nth systole. (M de¬ 
notes complementary memory.) Note that relay KO 
prevents this memory chain from learning the end 
systolic blood volume. 


In order to determine the required functional relation¬ 
ships between q^ and * n , a simple iteration circuit is 


provided. With some * n at the output of M3, systolic 

work is computed, and this is compared to the output 
of FG 2 which follows Starling 1 s curve. The error is 
repeatedly used to up-date * n until the desired * n is 

obtained. The above procedure is repeated for various 
values of v 0 until enough information is obtained to set 
up FG 1. With FG 1 in the circuit, a ventricular model 
is ready to be used as a part of a circulation model 
which will not be discussed here. 
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FIGURE 3 
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PROGRAM FOR ANALOG SIMULATION OF VENTRICULAR 
PUMP ACTION. 
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CHAPTER 22 

A PRACTICAL APPROACH TO ADAPTIVE CONTROL 


Originally published in: CONTROL ENGINEERING, May 1966 
Authors: J. W. Banham Jr. and W. L. Smith 


Controllers are linear elements that are often re¬ 
quired to operate in nonlinear systems. Thus they 
cannot be expected to provide optimum performance 
over a wide range of system operating conditions. 
However, through a linear representation of the non¬ 
linear system a controller can be designedwith adap¬ 
tive features that does provide optimum compensa¬ 
tion for transient system requirements. The follow¬ 
ing text describes the design of a controller and its 
conversion to the "adaptive" type (one whose dynamic 
parameters are continuously established as a function 
of transients in the system). The approach described 
has been breadboarded and tested for a shipboard sys¬ 
tem but the design techniques are applicable to other 
control systems. 

The maneuverability of a naval ship depends to a large 
degree upon the transient response characteristics of 
the steam generating and machinery systems. Most 
marine steam generators use forced draft; air for 
combustion is supplied to the fuel oil burners by vari¬ 
able speed blowers. The forced draft blowers in the 
machinery plant of a high speed combat ship are gen¬ 
erally powered by direct drive steam turbines and 
regulated by controlling the steam flow to the turbines. 
The rotational inertial characteristics of these ma¬ 
chines, over a wide range of operating speeds, contri¬ 
bute materially to the "sluggishness" of process re¬ 
sponse to changes in applied load. Designers interested 
in improving the system response must optimize the 
controller parameters of those system loops or sub¬ 
loops with the longer time constants. 

SYSTEM DESIGN STUDY 

In a conventional automatic combustion control sys¬ 
tem, controller settings are usually determined em¬ 
pirically at the load condition of maximum open loop 
sensitivity. Because controllers are linear elements 
and the forced draft blowers are nonlinear, the situa¬ 
tion is considerably less than optimum over a wide 
range of operating loads. 

a) Nonlinearities in the blower 

For any type of rotating machinery, the angular ac¬ 
celeration of the rotor assembly is proportional to the 
applied differential torque and inversely proportional 
to the polar moment of inertia. This linear relation¬ 
ship may be expressed by the simple differential equa¬ 
tion: 


» - 7 < T a- T r) 


where 

2 

w = angular acceleration, rad/sec 

2 

J = polar moment of inertia, lb-ft-sec 

T = driving torque developed by the prime 
a mover, lb-ft 

T^ = load torque imposed by the blower, lb-ft 

Several factors influence the applied torques. The 
torque developed by a steam turbine is a function of 
the turbine steam rate and its rotational speed. A 
typical turbine performance graph illustrating the func¬ 
tional relationship is shown in Figure 1. The para¬ 
bolic curve on this map represents the design load re¬ 
sistance for the turbine-fan unit. For this particular 
set of performance data, it was found practical to rep¬ 
resent the driving torque by the algebraic expression 

smax / \ w max / 


where 

G /G = turbine steam flow, fraction of 

s smax rating 

T = maximum locked rotor torque, lb-ft 

max 

m = a constant, lb-ft 

cj/w^ax = angular velocity, fraction of rating 

The steam flow to the turbine nozzles varies with the 
steam supply pressure and temperature and the steam 
valve flow coefficient and port area. The forced draft 
blowers considered in this study were equipped with 
a V-ported regulating throttle valve for which the 
characteristic flow-lift curve is as illustrated in Fig¬ 
ure 2. The supply steam state is reasonably constant 
throughout the load range. The load torque varies di¬ 
rectly with the blower system resistance and with the 
square of the rotational speed. 

T r “ P (comax) 
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LINEAR REPRESENTATION OF THE 
NONLINEAR SYSTEM 


All of the elements of the combustion air flow control torque terms may be reduced to linear functions ex¬ 
loop of amarine steam generating system are charac- pressed by the following equations: 

terized by some nonlinearity. The following analysis 
illustrates one method of representing the forced draft 

blower nonlinearities by a set of linear functions. AT = AG (T G mG„ w A ) 

as max so so o 

In the basic acceleration equation 



AT * 2Bco Aw 
r o 


the torque terms are functions of the turbine steam 

flow, G s and the turbine-fan speed, w, as shown in 

the text of this article. The accelerating torque, T , 

a 

was found to be a function of both steam flow and speed; 
the graphical relationship appears in Figure 1. In this 
study, the analytical relation between these variables 
was established by curve-fitting to manufacturers 
performance curves, but can also be developed by a 
detailed derivation of momentum equations applied to 
the turbine wheel. For a fixed blower air system re¬ 
sistance, the load torque varies with the square of the 
rotational speed. The constant of proportionality, /3 , 
is established by this resistance. 

Steam flow to the turbine is regulated by a throttling 
valve characterized by the arbitrary function 

G s = F(L) 


Gs is a function of valve lift as shown in Figure 2; a 

tangent constructed at the point (L , G ) has aslope 
. * o so 

given by 


AG„ 



L 


o 


AL 


Substitution of these linearized functions into the ac¬ 
celeration equation and applying a Laplace operator 
yields 

Aa> = Js [dL F ^ II = L L ^ T max G so 

L o 

- mG„ 0) ) - 2/3W ACJ 
so o' o J 


where L represents valve lift. 

Sufficient information is now available to permit the 
development of a set of transfer functions which de¬ 
scribes the open loop frequency response of the tur¬ 
bine-blower unit to small perturbations about a fixed 
operating level. The general perturbation response 
is given by 

A« = 7 / (AT a - AT r )dt 

If a small disturbance is applied at a particular steady 

state initial condition (L , G , oj ), the nonlinear 

o so o 


For any particular steady state operating level the 
terms d/dLF(L) Il=L q , T max G so , mG^ and 2p v Q 

are constants; the above equation may thus be expressed 
as a transfer function: 

“r- F(L)rT G - mG w 1 
Ao> _ dL v L max so _ so oj 

AL (b) “ 1 + [J/2/lcuJS 

This transfer function defines the response of shaft 
speed to valve lift for small perturbations given initial 
conditions defined by L q , G gQ , and 
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oj/qjxk ox Speed, fraction of rating 


Figure 1. Forced draft blower turbine performance 
curves. 



L/L m0 x Valve lift, fraction of full stroke 


Figure 2. Turbine valve steam regulating charac¬ 
teristics. The arbitrary function G = F (L) used in 

s 

Figures 3 and 7 is illustrated. Construction of a tan¬ 
gent at (G /G = 0.30, L/L = 0.53) illustrates 
s S max 

max 

the graphical approach to the linear representation of 
a nonlinear system. The tangent slope represents a 
valve gain of 1.12 percent of steam flow per percent 
of valve lift at a point corresponding to 53 percent of 
lift. 


where p is a constant for a parti cular system resist¬ 
ance, lb-ft. In the design study, the coefficient p 
was permitted to vary with the effective fan discharge 
area (established by the number of burner air regis¬ 
ters in use in the steam generator). 

This information led to the model of the system shown 
in Figure 3. The problem at this point was to devise 
a control law which would provide optimum compensa¬ 
tion for the nonlinear system over a wide range of 
operating speeds. For this design, the system is in¬ 
itially described in the frequency domain as a contin¬ 
uous series of linear systems, each similar in form 
to the others but having different parameters. 

b) Linear representation of the nonlinear system 

One approach to the linear representation of a non¬ 
linear system is to expand the nonlinear elements of 
the system by means of the Taylor series, then dis¬ 
regard all but the linear terms. This produces a 
linear system which approximates the real system re¬ 
sponse to small perturbations about a fixed operating 
point. Because the steam admission valve character¬ 
istic is an essentially arbitrary function, a graphical 
approach was used. In this method, all of the non¬ 
linear terms are represented as linear functions by 
constructing tangents at selected operating points along 
the curves which represent the functions of the alge¬ 
braic variables. The slope of each tangent curve is 
taken as the gain of the element under consideration, 
Figure 2. The resulting linear model is illustrated in 
Figure 4. The development of the values of the pa¬ 
rameters associated with this model is detailed in the 
box on page 22-2. n 

A linear model like the one in Figure 4 leads to a set 
of transfer functions instead of a single function. 
Each function in the set represents the open loop fre¬ 
quency response of the system to a small signal si¬ 
nusoidal disturbance superimposed on a particular 
steady state initial condition. The validity of the 
analytical model constructed in this fashion was con¬ 
firmed by comparison with actual frequency response 
data obtained from the operating units subjected to 
pulse test analysis at two load conditions, (Ref. 1). 
The family of transfer functions thus obtained was 
used as the basis for a series of control system de¬ 
signs /based on classical frequency response tech¬ 
niques applicable to linear systems. 

DESIGN OF THE CONTROL SYSTEM 

The following conditions were used as criteria for the 
dynamic performance of the closed-loop system: 

a. There must be no offset error in the steady 
state* 

b. No overshoot in the transient response to a 
"ramp" input disturbance is to be permitted. 

c. The time integral of the system absolute error 
is to be minimized. 

The conventional control loop employs a proportional- 
plus-integral controller which receives as its feed¬ 
back the output signal from a combustion air flow 
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DEVELOPMENT OF CONTROLLER DESIGN 


Conventional controllers used for regulating forced 
draft blowers are simple proportional-plus-integral 
devices characterized by fixed gain and reset rate 
over the entire operating range. The output signal 
of a linear controller of this type is given by 


e 


y e 



/ 


edt 


where e = error between input and feedback signals 
y - controller proportional gain 

— = controller reset rate 
r 

An adaptive feature was incorporated in such a con¬ 
troller by servoing the proportional gain and reset 
rate (integral crossover frequency) settings as re¬ 
quired to maintain an optimum combination of closed 
loop system response and stability over the designed 
range of forced draft blower operation. These con¬ 
troller parameters are shown plotted against various 
loop operating levels in Figure 6. Curve fitting of 
these graphical functions produced the following alge¬ 
braic expressions: 


y, : 


0 )/ 0 ) 


max 


where y' = variable controller gain 


y = optimum controller gain at = a> 


max 


where — = variable controller reset rate, rad/sec 

T* 


= optimum controller reset rate 

at«/ w = 0 
max 


= optimum controller reset rate 


max 


at w/w max = 1 


The output of the adaptive air flow controller may 
therefore be expressed by the equation 

"'■ 7 (^r) e *t/ edt 

* W *)/“* 


It is to be noted that the controller parameters have 
been expressed as functions of the control feedback 
signal (in this particular system, the blower rota¬ 
tional speed). To produce maximum loop gain during 
any transient disturbance, the controller gain and 
reset rate are automatically varied as a function of 
either the command signal $ or the feedback signal 
whichever of the two is least. An arrangement which 
produces this effect is illustrated in the computer sche¬ 
matic, Figure 7. Using the notation indicated on this 
diagram, the controller equation becomes 


w/w max = steady state loop operating level 

1 = 1 _ _1 

T T T 

o max 
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Figure 3. Forced draft blower block diagram. 



Figure 4. Block diagram of forced draft blower perturbation model at (Lq , ojq)* 


transmitter. For the "adaptive" design, it was elected 
to utilize the same basic controller concept but to vary 
the proportional gain and integral cross-over frequency 
to maintain optimum performance. 

The conventional design was initiated by selecting con¬ 
troller parameters which would give the entire control 
loop a gain margin of six db and a phase margin of 35 
deg; these margin goals were derived from previous 
experience. This conventional procedure was repeated 
for each of the forced draft blower transfer functions 
at the various operating load levels. The result, de¬ 
tailed in the box on page 22-4, wasasetof empirically 
optimum controller settings which could be correlated 


with the loop input or output to develop analytic func¬ 
tions that would specify the adaptive nature of the 
design. 

The air flow controller thus obtained was simulated 
together with the forced draft blowers and air flow 
system on an analog computer. The transient response 
to input disturbances comparable to the effects of ship 
maneuvering were then studied. Typical results are 
shown in Figure 5. Note that a small amount of over¬ 
shoot occurred which violates one of the system de¬ 
sign criteria. This observation led to a repetition 
of the design procedure, this time using gain and phase 
margins of 9 db and 45 deg, respectively. Figure 6 
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Figure 5. System transient response with initial 
optimization criteria. 




Figure 6. Controller parameters (servoed controller 
settings) as a function of loop operating level in the 
steady state. Open loop gain and phase margins are 
9 db and 45 deg. 


illustrates the relationship between the controller 
parameters and the loop operating level for perform¬ 
ance with the revised open-loop margins. 

A revised analog computer program was then pre¬ 
pared; a simplified version of the basic program is 
shown in Figure 7. Transient response of this simu¬ 
lated system appears in the chart recordings of Fig¬ 
ure 8A; the corresponding response of the simulated 
adaptive control system appears for comparative pur¬ 
poses in Figure 8B. It is evident from these record¬ 
ings that the adaptive system is considerably superior 
to the conventional system. Inspection of the error 
recordings discloses a reduction of five to one in the 
time integral of absolute error. At the same time 
there is neither overshoot nor permanent offset. Thus 
the design criteria have been met. 


FROM SIMULATION OF ACTUAL TESTS 


Results obtained from the analog computer studies 
supported a decision to breadboard a pneumatic con¬ 
trol system from standard hardware. The basic com¬ 
ponents used in the adaptive controller were pneu¬ 
matic proportional and integral amplifiers with servo- 
driven gain settings. Two of these units provided the 
necessary proportional-plus-reset control mode. The 
gains of these two devices were automatically varied 
as a function of either the input to the closed air flow 
control loop or its output, whichever produced the 
greatest loop gain during any transient disturbance. 
In the steady state, these variables are equal and 
produce the same gain in the open loop. 


Tests of the loop were conducted by imposing "ramp” 
inputs which corresponded to those used in the com¬ 
puter simulation, thus aiding in an evaluation of con¬ 
trol system performance. The recorded results of 
one of these transients are shown in Figure 8C. The 
transient performance of the closed loop system 
closely paralleled that of the simulated system, see 
Figure 8B. 


The air flow control loop comprises one of the open- 
loop elements of the steam generator’s automatic 
combustion control system. The principal benefit of 
the adaptive air flow control system when used as a 
part of the overall steam generator control system 
is its ability to reduce the deviation of controlled 
steam pressure and to eliminate the possibility of 
incomplete combustion (and accompanying smoke) 
during maneuvering. Although hardware requirements 
for systems of this type are presently more com¬ 
plex, and less reliable than for conventional systems, 
adaptive design techniques become attractive when 
applied to systems constructed from encapsulated or 
modular design control circuitry. The U. S. Navy is 
currently investigating adaptive systems of the solid 
state pneumatic (pure fluid) and solid state electronic 
types for use in the control of shipboard main pro¬ 
pulsion and auxiliary machinery. 
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Figure 7. Simplified analog computer program used to study conventional 
and adaptive controller responses. 
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Figure 8. Transient response recordings. A-simu¬ 
lated conventional controller. B-simulated adaptive 
controller. C-breadboard adaptive controller. 
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Glossary of Abbreviations 


C Capacitor 

C Compute 

C Clock 

CK Check 

D or DN Down 
DCU Decimal Counting Unit 
DFG Diode Function Generator 
DVM Digital Volt Meter 
ES Electronic Switch 

Fb Feedback 

FF Flip-flop 

F-RT Forward-Reset Real Time 
FR Function (Operational) Relay 

FS Function Switch 

H Hold 

I or ic i nitial Condition 

IJ Initial Junction 

IN Input 


10 Iterative Operation 

J Junction 

K Relay 

L Logical Variable 

M Meter 

NC N ormaI ly Closed 

NO Normal ly Open 

0 Output 

P or POT Coef. Potentiometer 
R Reset 

R Resistor 

R-H Reverse-Hold 

RO Repetitive Operation 

RTS R eset-T oggle-Set 

S Set 

SJ Su mming Junction 

T Toggle 

TL Trunk Line 


INV Inverter 


u 


Up 




APPENDIX II 


The symbology used throughout this handbook is in conformity with the objectives 
of Simulation Councils, Inc., namely, to encourage the use of uniform analog com¬ 
puter graphics. This article was prepared by SCI, and appeared in the March 1966 
issue of SIMULATION. Permission has been granted to reprint this article. 


Presenting... 

Uniform graphics for SIMULATION 


The symbols and the methods of laying out analog 
and hybrid computer diagrams presented here are 
advocated to alleviate the confusion caused by the 
uncoordinated invention of new symbols and dia¬ 
gramming practices. The increasing use of hybrid 
techniques and equipment has aggravated an already 
bad situation to the point that it is often no longer 
possible for one worker in our field to read another's 
diagram. Usually this is because symbols are devised 
and diagrams are drawn to include details peculiar 
to a particular kind of equipment. Such a wiring, or 
"patching," diagram is of course necessary for set¬ 
ting up and checking out an actual simulation, but 
hardware-peculiar details are only confusing to those 
with other kinds of equipment. With few exceptions, 
the use of a simplified signal-flow diagram to illus¬ 
trate technical articles is much more effective. 

With the foregoing in mind an SCi committee 
composed of 

GEORGE BURGIN 
JOE HUSSEY 
HANS JORGENSEN 
GRANINO KORN 
JOHN McLEOD 

selected the symbols and offers the following sugges¬ 
tions for their use. Primary considerations were cur¬ 
rent usage, clarity, and simplicity. We devised no new 
symbols and, unless there were overriding indica¬ 
tions to the contrary, we adopted those already in 
widest use. Clarity and simplicity, we believe, will be 
enhanced by the choice of unique shapes to repre¬ 
sent different components, and the elimination of all 
unnecessary details in diagrams. 

There was no intent on the part of our committee 
to set up standards for the industry. However, all 
diagrams appearing in SIMULATION will be pre¬ 
pared according to the committee's recommenda¬ 
tions (as they may be modified from time to time), 
and we hope that these recommendations will prove 
attractive to others. Suggestions for modifications 
and additions are solicited. 


General rules 

The following methods and symbols are recom¬ 
mended for the illustration of technical articles pre¬ 
pared for publication. Unless the purpose of the 
article is to describe the use of a particular kind of 
equipment, and the hardware details are pertinent 
to the subject, such illustrations should not be "hard¬ 
ware-peculiar." In other words, the objective should 
be to show signal flow, rather than "patching" 
details. 

The primary, or overall, system diagram should 
show only the essentia! signal flow. Where it is neces¬ 
sary to show details, separate diagrams should be 
made and referenced to the primary diagram by en¬ 
closing the detailed area of the primary diagram in 
dotted lines with suitable notation. 

The direction of signal flow should be indicated 
by arrowheads except where the shape of the sym¬ 
bols makes the direction of flow obvious. Primary 
signal flow (with the exception of feedback loops) 
should be from left to right, and, if practical, each 
"line" of symbols should be made to read like the 
mathematical relation it represents. 

The choice of whether to end a line and label it 
(preferably with the symbol of the variable that the 
signal represents) when it reaches the right-hand 
side of the diagram, and then indicate its continua¬ 
tion with the same label as it enters again at the 
left-hand side, instead of drawing in the connection, 
should be made on the basis of clarity; if a line re¬ 
turning the signal from right to left will cross many 
other lines and be hard to follow, it should not be 
drawn in. 

If a diagram involves a number of identical cir¬ 
cuits, only one should be shown in detail, while the 
others should be indicated by boxes with appropriate 
notation. 

Components should not be numbered unless they 
are referred to by number in the text. 

All amplifier gains should be shown just outside 
the amplifier at the point where the input enters. 
Unity gains should not be labeled. 


Always apply the test of clarity and simplicity. Ask 
yourself: "Is this the most understandable way to 
diagram this for those unfamiliar with the hardware , 
and less familiar with the subject , than IV 
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ELEMENT 


SYMBOL 


FUNCTION 


NOTES 



Distinction as to kind of signal is made only at 
beginning and end of line, and then only where 
there are two or more kinds of signals to be shown 
on one diagram. Arrowheads should be omitted 
only if direction of signal flow is obvious. 


Analog 


Mechanical 


No connection 


Connections 


No arrows; shape clearly indicates direction of 
flow. 

No label on unity-gain input. 

"N" indicates where number of amplifier should 
be placed if, and only if, referred to in text. 


Amplifier 


Summers 


Label high-gain input "SJ" or "G" only if perti¬ 
nent; the point at which c enters (at top or bot¬ 
tom of symbol) indicates a gain of 


Summing amplifier (with high-gain input) 


Note that the numbers at the inputs are now only 
relative gains; /*, an unspecified large number 
which is a function of the particular amplifier, is 
the overall gain. 


High-gain amplifier 


When amplifiers are used only to change signal 


sign they should be smaller. Orientation should 
be that which makes for greatest graphical sim¬ 
plicity, i.e., they may be shown in vertical, or in 


Inverting amplifier 


feedback, paths. 


Standard integrator mode is that of basic problem. 
The point at which a enters (at top or bottom of 
symbol) indicates that it sets the initial value of x, 
so it is redundant to label it 1C. 


Standard 


Integrators 


•a —/(10b —5c —j— d)dt 


Arrows indicating mode control should be shown 
only when two or more integrators are operating 
in different modes simultaneously. The mode- 
control input should then be labeled or coded to 
indicate the controlling mode. 


RESET when R=1 
COMPUTE when R = 0 


Three-mode 


RESET when R = 1 

or R = 5=1 
HOLD when 5=1, R = 0 
COMPUTE when R=S=0 


When used in this way the size of the integrating 
capacitor may be important, as it affects the track¬ 
ing lag. In such cases the value should be indi¬ 
cated as shown. The quantity tracked or held is 
—(a -f- b + c). 


TRACK when R = 1 
HOLD when R = 0 
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ELEMENT 


SYMBOL 


FUNCTION 


NOTES 


Integrators 

(continued) 


Potentiometers 


Function generators 


Memory pair 



Integrator 1 in 
Track for 5=1 

Integrator 2 in 
Track for R = 1 


Integrator 1 in TRACK for 5=1, integrator 2 is 
TRACK for R = 1. in certain iterative problems y 
is required to have some value for the first itera¬ 
tion; thus an "initial" initial condition must be 
furnished as shown. 


Two-terminal 




Three- a 
terminal 


x 



x=fca 

(load connected) 

a = 1,6=0 
(load connected) 


The mathematical relationship shown for the 
three-terminal pot is given only to define k. It 
will not usually hold during problem solution be¬ 
cause if 6 were always zero a two-terminal pot 
could be used. If it is desirable to define k in any 
other way, its meaning should be clearly indicated. 


Servo-set 
a_ 


- 0 - 


x=ka 


This should be used only if the fact that N is servo- 
set is significant to the solution of the problem, as, 
for instance, when it must be automatically reset 
after each run of an iterative problem. 


Digital 




x=k*a 


As shown, the symbol indicates a digital pot with 
integral (committed) inverting amplifier. The sym¬ 
bol for non-inverting digital pots should not in¬ 
clude the amplifier. 

"Indicates quantized value. 


Arbitary 


a 

f 

hS x 


V 


x=f(a) 


The number N should be used only when the 
function generator is referred to by number in 
the text, or when it is desirable to explain what 
the function f is, either by a footnote, or by a 
separate graph. 


Mathematical functions (typical) 


a 

, i — 


V 

Multiplier 

a 


> 


-fa 


+b 


Divider 


\ >> 




=vs~ 


-x, x=- 


ab 


aa 

IT 


a =reference voltage 


In cases where the function generated can be rep¬ 
resented by a standard mathematical symbol, the 
t should be replaced by the symbol. 


These symbols should be used for all analog multi¬ 
pliers and dividers; if the kind is significant to the 
solution of the problem, it should be so stated 
in the text and/or noted on the diagram. Because 
equipment differs, the sign of the output for spe¬ 
cified signs of both inputs should be given, other¬ 
wise inversion will be assumed. 

The number N should be used only if the com¬ 
ponent is referenced in the text, or if more than 
one channel of a multichannel device is used. In 
the latter case the subscript, in this case /, identi¬ 
fies the channel. In many cases, it will be found 
that a drawing can be simplified and clarified by 
drawing the same multichannel device in more 
than one signal flow path. In this case the number 
would be the same, but the subscript would be 
different for each channel. 


Resolvers 




X x_ 

L. 


N jl 


x= R cos 0 


Polar to rectangular Rectangular to polar y—R sin Q 


$ = arc tan y/x 
R = Vx 2 -f y 2 


If a resolver has additional outputs they should be 
shown only if used. 


Comparators 

With binary output With true and false U=1 (a < 6) 

binary outputs U=0(a>6) 



u 


a 

- -> 

6 

N 

j 


The symbol for a relay comparator can, of course, 
be made by combining either of these symbols 
with that of a relay. 
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ELEMENT 

Relay 


Digital Inverter 

Digital/analog switch 

Limiters 


Converters 

Diode 

Diode gates 


Flip-flop 

Black box 



SYMBOL 


FUNCTION 


NOTES 



L4=x (U=0) 
U m =x (U= 1) 


If the diagram is drawn as shown, the designations 
of arms and contacts are superfluous. However, 
it is often desirable, for clarity, to show the relay 
coil in one place on a diagram and the contacts 
in signal-flow paths elsewhere. In this case, 
would designate the normally Open (U = 0) con¬ 
tact of channel a, relay N; the Arm of channel a 
of relay N, etc. 


t 


U a 


U a =a (U= 1) 


Feed-back 



x——a (/<a<u) 

x——/ (a < I) 

x ——u (a > u) 


x==a (/ < a < u) 
x— / (a < /) 

x—u (a > u) 


Analog-to-digital Digital-to-analog 

. x* a’ 



Solid state 




x==a (a>x) 
x=0 (a < x) 


AND 


NAND 





AND 

NAND 

OR 

NOR 


x=a-6-c 

x^a-b-c 

x==a -f- b -j- c 
x^a -f-6-f c 


The tilde ('-') is used instead of a bar to indicate 
negation. This is recommended because of the 
many other meanings of a bar over a variable. 

ON when binary signal, U, is digital "one." 


u— larger value at which output is limited 
/ = smaller value at which output is limited 



•Indicates quantized value. 


Output is high (logic 1) if and only if all inputs 
are high. 

Output is low (logic 0) if and only if all inputs 
are high. 

I Output is high (logic 1) if and only if at least one 
input is high [i.e., any or all of the inputs are high). 

Output is low (logic 0) if and only if at least one 
input is high (i.e., any or all of the inputs are high). 

As there are many kinds of flip-flops, symbolic 
representation of the operation of any but the 
simplest is not advised; if the operation is not 
obvious, an explanation in the text or a footnote 
with the diagram is recommended. 


BB N 


FT 


This is Black Box number N, to be used if no sym¬ 
bol is given here, or in case of doubt\ It can have 
any number of inputs and outputs, but the nature 
of the signals should be indicated and all signal- 
flow directions should be designated by arrow¬ 
heads. Here there are three analog inputs giving 
rise to one analog and one binary output. 


The function of a "black box" can often be 
made obvious by properly labeling the inputs 
and outputs. However, if it is not obvious, the 
function should be explained in the text or by a 
footnote with the diagram. If internal details are 
of interest, they should be shown in an appro¬ 
priately labeled auxiliary diagram. 
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